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MOSQUITOES AND MALARIA. 


A STUDY OF THE RELATION BETWEEN THE NUMBER OF 
MOSQUITOES IN A LOCALITY AND THE MALARIA RATE, 


By H. WAITE, M.A., B.Sc. 


PROFESSOR Ross, in his Report on the Prevention of Malaria in Mauritius*, 
has dealt briefly, from a mathematical point of view, with the influence of 
mosquitoes on malaria and he suggested in a letter to Professor K. Pearson that a 
further treatment from the more mathematical side would be of value. Such is 
the origin of the present investigation and it is to the above-mentioned Report 


that I am entirely indebted for the numerical data upon which this enquiry is 
based. 


It has been proved conclusively that malaria is caused by the bite of 
a mosquito which has already bitten a previously infected person; and that the 
infection is probably conveyed only by the sub-family of mosquitoes called 
Anophelina. Hence, in order that the disease may spread in a neighbourhood, 
there must be present a source of infection in the shape of persons suffering from 
malaria, and a sufficient number of anophelines to act as carriers. Moreover, 
since recoveries are constantly taking place, it is clear that the increase or 
decrease in the number of cases will depend largely on the relative frequency of 
anophelines and of human beings. 


Malaria Rate. The ratio of the number of persons affected with malaria 
to the total population of a district at a given time is called the Malaria Rate of 
the district at that time. In general, the rate is continually changing owing to 
(a) new infections, (b) recoveries, (c) emigration and immigration, (d) the birth 
and death rates, and (e) the extent to which cases are isolated, as well as owing to 
changes in the mosquito population. 


As emigration and immigration vary considerably in different localities, and in 
the same locality at different times, their influence on the malaria rate cannot be 


* Report on the Prevention of Malaria in Mauritius, by Ronald Ross, D.P.H., F.R.C.S., D.Sc., 
LL.D., F.R.S., C.B., published by Waterlow and Sons, London, 1908. 
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satisfactorily dealt with except in particular cases where the necessary statistics 
are available; neither would results in general terms be of much practical use. 


It is probable that the risk of infection may depend partly on the age and 
past history of an individual; for example, a person who has completely recovered 
from one or more attacks may become partially or entirely immune. So far, 
however, I have not found any statistics which deal with this phase of the 
question. 

The combined effect of the birth and death rates will, as a rule, be small 
compared with the effects under (a) and (b); eg., the birth rate for Mauritius 
in 1906 was 2°8 per thousand per month and the death rate 3°3; of the deaths, 
1:2 per thousand were due to malaria. These figures do not differ appreciably 
from the average taken over a long period of years. 


I propose, therefore, in the present paper, to confine my attention to the effects 
on the malaria rate of the first two of the above causes and to leave for a later 
consideration the minor contributory causes. 


I shall also assume in this preliminary survey that both the human population 
and the anophelines are uniformly distributed over the area in question; that 
persons of all ages, whether they have previously suffered from malaria or not, 
are equally liable to infection; that the number of anophelines per unit of the 
population remains approximately constant during the period under review; and 
further that infected and non-infected persons are equally exposed to the risk of 
being bitten. 


It is obvious that the results so obtained will require modification in practical 
application and ought to be taken into account when the theory is developed to 
a higher approximation ; they will be, at the best, tentative and suggestive, but 
it is hoped that they will tend to throw light on the problems which arise when 
the specialist has to consider remedies for the disease. 

Infection Rate. The number of new infections in a locality during any 
particular month has been called the Monthly Infection Rate. This depends 
chiefly on the number of infected anophelines in the neighbourhood, that is, the 
number which have bitten previously infected persons. 

I shall in general follow the notation adopted by Professor Ross (vide Report, 
p. 31, et seq.). 


Let p denote the average population of the locality during the enquiry, m the 
number of infected persons*, h the number of healthy persons at the beginning of 
the first month and im the average number of the infected persons whose blood 
contains enough of the malaria parasites to infect anophelines; evidently 7 is a 
fraction. Next, let a denote the average number of anophelines in the locality 
capable of carrying malaria and ba the average number of these which succeed in 
biting a single person during the month; here b must be a very small fraction, 


* This differs from Professor Ross’ use of m. 
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since probably many mosquitoes never succeed in biting men at all, and of those 
that do succeed, only a few will have bitten a particular individual. Hence the 
average number of anophelines which have bitten the suitably infected persons 
during the month will be baim. 


Now if any of these insects are, in their turn, to infect human beings, they 
must survive for at least a week or ten days, in order to give time for the 
parasites to mature within them; and by no means all of them will survive so 
long. Let s be the average proportion of anophelines which can survive for that 
period; then sbaim will be the number of infected mosquitoes which have 
survived long enough to infect men in their turn, where s is also a fraction. But 
not all of these will find opportunities to bite human beings again, though they 
have survived long enough to do so.. Let f be the proportion which succeed in 
biting. Then fsbaim will be the average number of infecting anophelines which 
succeed in biting men; also, if each of these bite a separate person and only one 
person, the same expression will denote the average number of persons bitten by 
infecting mosquitoes during the month. Call this number n. The most probable 
distribution of the n bites between malarial and non-malarial persons will be in 
the ratio m:h; so that the number of new infections during the month will be 
fsbaim h/p, or nh/p*. 


Now each time a healthy person is bitten by an infecting mosquito the number 
of malarial patients is increased by unity, while the bites inflicted on those who 
are already suffering from the disease make no difference in the number of cases. 
This is equivalent to increasing m and decreasing h by h/p for each bite. 

[By taking mean values in the above I am assuming that if 

z= F(a, &, ... @y); 
then Z may be taken = F'(%, %,,...%,), and I have been led to make this 
assumption in a first investigation by finding that the modal values of the 
distribution approximate very closely to the means in a number of particular cases 
fully worked out.] 


Recovery Rate. The ratio of the number of recoveries to the total number of 
malaria cases during any period is called the Recovery Rate for that period. 
Professor Ross has estimated that of a given number of infected persons 50°/, are 
ill after three months, that is, the recovery rate is ‘5 per three months. 

If m, be the number of cases at the beginning of the period, m,, m., ms, at the 
ends of the first, second and third months respectively, and r the recovery rate 
per month, 

Mm, =M— mr =m (1-71), 


m,=m,(1— 7) 


m (1—ry, 
m,=m,(1—-—r)=m,(1—ry. 
* Professor Ross assumes that to a first approximation m will be the number of new cases in the 
month. 


54—2 
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But Ms = m,(1 — °5), 
“. (l—rf= 55, 
log (1 —r) =} log °5 


1 

3 
= 1:899657, 
1—r= 7937, 


and . 9” ='2063 per month. 





In the numerical examples it will be necessary to know the recovery rate for 
the average period between two consecutive bites. Suppose there are n bites in 
a month and r is the recovery rate for the period between two consecutive bites, 
then, as before, 

(l-r) 


I 


1 — ‘2063 


= ‘7937, 


1 
log d-r)= ra log ‘7937 





. ne 
=- x 1899657 
nt 
"100343 
—— > 
it 
whence the value of 7 is readily obtained. 
I will now consider a formula for giving the number of cases at the end of | 
a month in terms of the number at the beginning, the total population, the f 
number of anophelines and the recovery rate. | 


Let m, be the number of cases at the beginning of the month out of a popula- 
tion p; m the number after the ¢” bite and r the recovery rate when there are 
n infecting bites in the month; then | 


P= Ti3.5 
Me= M_, 4 I —7Tm-, 
P 
l 
or mMm=mM-,|1—---r)4+1. 
». 2 


] , ; ’ 
Put #& for (1 “—" r) and this equation may be written 
} 

m=mR+1. 


Now r remains constant for the month under consideration and therefore 
also R. Hence 
M4, = Me» R T E 


m,=m,R +1, 
m, = Mp R+1. 


From these equations 


1— Rt 
1-R’ 


m, =m, Rt + 
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and therefore, for the end of the month 


1 

My = MR” + i- Ro" dens deivakeunessiaaaecs (A). 
Now, for any particular month n = fsbaim where m is the number of cases at 
the beginning of the month, (or the average number for the first half of the 
month and the second half of the previous month); and Professor Ross has 


estimated the approximate values of the constants f, s, b, ¢ to be i: >< and 4 
o =p 
respectively ; hence 
eS 
"Stes m 
_ am 
~ 192p’ 
also r ='2063 when n } 1, 
: 34: 
and log (1 —r)=- oe = when n ¢ 1. 
U 


We can thus by successive applications of the formula (A) obtain a series of 
values of m at the ends of consecutive months when we know the value at the 
beginning of the period, the total population and the average number of 
anophelines in the neighbourhood. Although it is impossible to estimate this 
last quantity with any degree of accuracy, yet by comparing the results of 
numerical examples, using different values of a, we get an idea of the effect on the 
malaria rate of reducing the anophelines in any given proportion. 

Condition for a Stable Population. Suppose that a certain locality has a 
population p and m, cases of malaria. This latter will remain stationary if the 
number of new cases is equal to the number of recoveries. With a given value 
of m, and p the number of new infections depends on the number of anophe- 
lines in the neighbourhood. Hence with the anophelines present in a certain 
proportion to the human population, the malaria rate is stationary and if they are 
present in a greater or less proportion there will be a corresponding increase or 
decrease in the number of cases. 


Let m, be the number of cases at the end of the month ; then, by Equation (A) 


1 — R* 
m, = mR” + ae 
The value of m is stationary if m, = m =m, Le. 
1 es R 
ee ed 
1 \1-—R)=0 
or (m ey (1 — &£”")=90, 
either | ee ad) | eee nank Guigbla dale oesea ele 
1 
or m—— = () 
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If 1—R"=0, n=0 for R=1 -; —r and is therefore fractional; also if n=0 


the case reduces to that in which there is an entire absence of mosquitoes: hence 
in this enquiry 1— R” #0 and we must proceed with Equation (ii). This may 
be written 


Now when n>1 


log (1 - r)=- 100868 ; 


n 








Hence by giving a series of values to m we may find the number of anophelines 
corresponding to any stable population. 
Taking p= 1000, m=5, we get 
r= ‘2 — ‘001 = °199, 
log (1 — 7) = — 096367, 
100343 _ 096367, 
n 
and n = 104125. 
But <—— ; 
192p 
= 192 x 1000 x 1:04125 
5 


= 39984. 


Similarly we get the table of corresponding values of m and a for stable 
g p g 
populations given below. 


TABLE A. 


Table of Corresponding Values of m and a for a Stable Population 
of a Thousand. 

















| 
m a m a m a 
as) eee Salt : 
5 39984 150 | 52032 600 110914 
8 41886 175 53648 650 126720 
10 42553 200 | 55339 700 147892 
15 | 43539 250 | 59062 750 177525 
20 | 44148 300 | 63308 800 221752 
25 | 44605 350 68193 850 | 295897 
| 50 46251 400 73883 900 | 444118 
75 47662 450 80612 950 | 885583 
100 | 49069 500 88680 | 
125 | 50520 550 98534 | 














H. WartE 427 





Again, suppose m = 5 and that a has such a value as will make n = 1, ice. 


da Ai ne 
192000 1 or a= 38400. 
We now have r=-2063 and rm=1°0315; hence rm>n; or the number of 


recoveries exceeds the number of new cases and the value of m decreases 
continuously. 


From the curve connecting these quantities (Diagrams I and II), we see at 

a glance the value of m corresponding to any given value of a and we gather that 

any appreciable change in the number of anophelines will involve a corresponding 
gradual change in the number of cases; or, in other words, that, given any 

definite number of anophelines per thousand of the population there is a particular 


malaria rate which will ultimately be reached by a gradual change in the 
| value of m. 











Number of Anophelines in 1000's. 
\ 








38 [ a a a a Y Ee ee ee ey —— 
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 3 
Number of Infected Persons. 


Diacram I.* A. Stable Population Curve, Part 1. B. Curve of Equation (ii) p. 429. 


| * See p. 429, for the construction of these curves. 
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) 560 58! 


1 4 4 ri — ee ttle ential 
0 600 620 640 660 680 700 720 740 760 780 800 820 840 860 880 900 


Number of Infected Persons. 


Dracram II.* 


When n <1, we have (p. 424) 


(1 —r)" = 7937, 


or 
Now for 


also 


a 


stable 


log, (1 —r) 


population 


See p. 429 for the construction of this curve. 


Stable Population Curve, Part 2. 
direct relation between m and a (these symbols denoting 
values for a stable population) may be investigated as follows: 


= log, ‘7937 
“ie: 


23104 
n° 
“25104 
l-r=e ” 
] 1 
r= _—-, 
m Pp 
am 
n= 
192p’ 
1 1 ‘23104 x 192p 
— = @é@ am 
p m 
1436p 
= e am 





corresponding 
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Put c for 44°36 and the equation is* 





429 
:.4 = J : 
oo a SS  nccetcisnnrer eoccceceees sowense 
pm 
Expanding the right-hand side we get, to a first approximation, 
1+ 1 — 1 =]— a ; 
p m am 
which reduces to 
pre 

a = —*— 

P —m 
A table 


of values of a and m found from Equation (ii) is given in Table B 
below; the i 


graph is the upper curve in Diagram I, the portion of the graph 
common to this equation and to Table A, is given in Diagram II. 


The graph of Equation (ii) fits very closely that of Table A when m exceeds 
100, and accurately when m approaches p; for small values of m however the two 





* The general form of this curve is given below. 


* : 0) on the m axis 
pt+l1 

(b) in the second quadrant with the axes of coordinates as asymptotes; (c) in the fourth quadrant, the 

axis of m and the line m=p being asymptotes, That part of the branch (a) which is above the hyperbola 

am=192p, alone has any direct application to the present problem, 


It has three branches :—(a) in the first quadrant, 
the line m=p being an asymptote ; this branch ends in a point d’arrét at ( 














Diacram III. 








SSIS 





General Form of Curve of Equation (i). 


The dotted line is the hyperbola which limits 
the practical solution of the part of the curve in the quadrant AOM., 
Biometrika vu 
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curves differ considerably and it is necessary to proceed with Equation (i) to 


a higher degree of approximation. 


TABLE B. (Equation (ii).) 





m a 
—_ — 

0 44360 
10 | 44808 
50 46695 

100 49288 
150 52188 
200 55450 
250 59146 
300 63371 











Writing Equation (i) in the form 


ke 1_ ~) 
am loge (1 + Pp m 


and expanding the right-hand side, we have 


pe _ (1 ~\ = 5 1_ily,! 3. sy 
“Gas m 2\p = 3 (5 =) m 


cee 1 —my? 1 = 
or os = — + 5 (p te += (P = a a 
am pm a paw 2) pw 


: y—m : : 
Putting z for — , this may be written: 
pm 





pe —_" o eS « . , 
ae + 5 23+ 9 2 approx., 
ss, am _ 1 (1 ,! raise i 
po 2 : = ) 
7 ae t 3 
‘si Ss SS + (5045: ‘>> 
Ste Ot a, AG 
= 4 1 _ a= z 2 
z | 2 * 4 sa ) 
= ee See 
= 55 — Te APProx., 
+. q-—P’c__ pe _¢(p—m) 


p-m wm 12m? 


The above relations cease to hold good when n<1, i.e. when am <192p; 
hence the limit of applicability of the Stable Population Curve, obtained from 


Equation (iii), is given by its intersection with the hyperbola 


am = 192p. 
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Part of this hyperbola is shown on the left of Diagram I as a broken line 
from which it appears that in a population of a thousand, there can be no stable 
condition with a smaller average-number of infected persons than about 45. 


The numbers in Table A are all found for a population of a thousand, but the 
general form of Equation (iii) is applicable to a population of any magnitude so 
long as the number of infecting bites does not fall below one per month. For 
example, when p = 20000, and m = 200, Equation (iii) gives a = 893941, while the 
formulae on p. 426 give a = 893966, the difference being less than ‘003 per cent. 

It will be of interest to compare the values of a found from Equation (iii) with 
those of Table A in a few. typical cases : 


m a, by Equation(iii) a, from Table A 
10 42554 42553 
50 46249 46251 
| 100 49067 49069 | 
500 88676 88680 





c(p—m) 
12 


In the last three cases the term ,— is negligible, while in the last case 
m 


me. : 
the term — is only about ‘01 per cent. of the whole. 
m 
Numerical Examples. In the numerical illustrations which follow, the crease 


or decrease in malaria, month by month, has been obtained by successive applica- 
tions of Equation (A), p. 425, in the form 








Ms = Ms_,R,"% + ; ~ as 
where Ns = isap? 
R,=1 -5-t 
log (1 —r;)=—- _— [n>1], 
r, = *2063 | [mn > 1}. 


Seven-figure tables have been used throughout in the calculations but the 
results are given to one place of decimals only. 


(Note. It is assumed, in all the examples, that Professor Ross’ recovery rate 
holds. This would rapidly reduce a large number of cases to the vanishing point 
in the absence of new infections, whereas it is well known that a few persons 
suffer from relapses after living many years in countries where there is no 
possibility of re-infection. (See Report, p. 16.)] 

55—2 
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Example 1, Case (i). p=1000, m,=50, a=64000. (The number of anophe- 
lines for a stable population is about 46000.) 


TABLE C. 


(Giving the value of m for consecutive months.) 


| | 
Months | 1st Year le Year | 3rd Year | 4th Year | 5th Year | 6th Year | 


| 


| aan | 


1 539 | 1180 | 1988 | 2585 | 286 | 297-7 
2 580 | 1246 | 205°0 | 261°7 | 287-9 | 298-2 
3 | 624 | 1313 | 2110 | 2649 | 289-2 | 298-5 
5 670 | 1380 | 2168 | 267°3 | 290°3 | 299-0 
5 71:7 | 1449 | 292-4 | 270% | g91-4 | 299-4 
6 76-7 | 151°8 | 227-7 | 2731 | 292-4 | 2996 





7 82-0 158°6 233°0 275° 293°3 300°1 
8 87°5 165°5 2380 277°4 294°3 300°4 
9 93°2 172°6 242°6 279°5 295-0 300°6 

10 | 99°1 179°5 247-0 281°5 295°7 301-0 

11 | 105°2 186-1 251°2 283°3 296°5 

















12 | 111°6 192°6 254°9 284°9 297°1 








Example 1, Case (ii). p=1000, m=50, a= 24000. 





TABLE D. 
Months Ist Year | 2nd Year | 3rd Year | 4th Year | 
1 45°] 13°8 5:8 4:98 | 
2 40°7 12°6 5°7 | 4:96 
3 36°8 11°5 5°5 | 4°94 
j 33°2 10°5 5:4 4°92 
5 30°1 9°6 53 | «(4:91 
6 27°2 8°8 53 | 4:90 
7 24°7 8:0 52 | 4:89 
8 22°4 74 5°1 | 4°88 
9 20°3 6°9 5°] 4:87 
| 10 18:3 6°5 51 487 | 
| 11 16°6 6°2 50 | 4:86 
12 15°2 6°0 50: | 4386 





It will be seen from this example that with a comparatively small number of 
infected persons and a quantity of anophelines somewhat in excess of the number 
which would keep the population stable, we get a steady rise in the malaria rate 
until nearly one-third of the whole population are infected; on the other hand, 
the destruction of rather more than half the mosquitoes will produce a fall in the 
rate until the cases are so few that they can be easily and effectively dealt with. 


For the graphs see Diagram IV. 


Example 2. p=1000, m,=500, a =(i) 96000, (ii) 48000, (iii) 24000, (iv) 0. 
(The number of anophelines corresponding to a stable population would be about 
89600.) 





—— 
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Number of Infected Persons. 
& 
~“ 











2 <6 8 10 12 14-16 18 20 29 ga 96 08-30 BP 3a Se a2 4s 46 a8 G0 02 8a G6 88 60 6? 64 ce 
Months. 
Dracram IV. Curves illustrating Example 1. A. Case (i). B. Case (ii). 
In Case (i) we get an increase and in each of the other cases a decrease in the 
value of m according to the following Tables : 
TABLE E. TABLE F. 
Case (i), a=96000. Case (ii), a = 48000. 











l 
} | Months | lst Year } Ist Year 2nd Year, 3rd Year 4th Year | 5th Year 
| 
| ae aa : i a 
| 1 | 507°9 455°3 | 235-0 169°7 1389 | 121°3 
2 | 513-9 118°7 | 227-0 166°3 137°0 | 120°2 
3 518°7 388°3 219°6 163°1 135°3. | 11971 
| 4 | 5292-7 362°4 | 212°8 160°2 133°6 | 118-0 
| 5 926°3 340°2 206°5 | 157°3 132°0 | 117°0 
| ¢ 528°8 321:0 | 200°7 154°6 1305 | 116-0 
7 530°7 304°1 | 195°3 152°0 129°0 | 115-1 
| 8 532°3 289°3 190°3 149°6 127°6 | 114-2 
9 | 583-5 276°0 185°6 147°2 126°2 113°3 
10 934-4 264°2 181°3 145°0 1249 | 1125 
11 535-2 253°5 177°2 142°9 123-6 | 111°7 | 
12 | 535°8 243°9 173°3 140°8 122-4 | 1109 
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TABLE G. 
Case (iii), a= 24000. 


TABLE H. 
Case (iv), a= 0. 














| ' 

Months | Ist Year | 2nd Year 3rd Year | 4th Year| 1st Year | 2nd Year | 3rd Year | 

a) | | Bo Se: See! Pe ioe 

| 

| 1 427°0 | 94:8 | 27° L 8°8 396°9 24°8 16 | 

2 369°8 | 85:0 24°7 8°1 315°0 19°7 
= 320°8 | 76:2 224 | 2500 | 156 
| 4 279°8 | 68°5 20°3 | 198°4 12°4 
ze . 245'°3 | 61°6 185 | 157°5 9°8 
a. 215°8 | 55°5 168 | — 125°0 78 

y 190°6 | 50°0 15°3 | 99°2 6:2 

8 168-8 | 45°2 13°9 a 78-7 4°9 

9 149°9 40°8 127 | — 625 | 39 | 

10 132°9 | 36°8 ne | - — 496 | 31 | 

11 | 118% 33°3 05 | — 39°4 | 2°5 

2 105°9 | 30°1 ~~) = 313 | 19 | 

| | 


This example serves 





to show the effect (i) of leaving the mosquitoes alone, 


(ii) of destroying half, (iii) of destroying three-quarters of the mosquitoes, 
(iv) of exterminating them (but see Note on p. 431); the comparison is best seen 


by means of the 


540;-—r 


curves in Diagram V. 





r + — r 
—_——— 


520 


Ze 


Number of Infected Persons. 
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Diacram V. Curves illustrating Example 2. 





A. Case (i). 


B. Case (ii). 


C. Case (iii). 


D, Case (iv). 











H. Watrre 435 


The foregoing results point to a stable population for any given number of 
mosquitoes, but it is perfectly clear that this number is itself subject to constant 
change, often of considerable magnitude, and that, in consequence, a stable 
population is rarely, if ever, attained in practice. 


The minor contributory causes mentioned on p. 421, will also affect the 
malaria rate to an extent which cannot be neglected, while another important 
item to consider is the distribution both of the human population and of the 
mosquitoes; also whether the most densely populated districts are at the same 
time those most infected with mosquitoes. These and many other points will 
require careful investigation with the aid to be derived from the most recent 
and complete observations before the subject can be said to be exhausted, and 
I am quite well aware that the present paper merely traverses the outskirts of an 
important and extensive field of research. 


But still it seems to me to establish two very important points, i.e., (i) Given 
a number of mosquitoes greater than that corresponding to the “stable population 
value,” the number of malaria cases will tend to increase until a stable population 
value is reached, (ii) Given a number of mosquitoes less than that corresponding 
to the stable population value, the number of malaria cases will tend to decrease 
until that stable population is established. 


In both cases the number of malaria cases tends to asymptote to the stable 
value. The amount of malaria does not increase or decrease indefinitely, but 
tends to attain a definite prevalence. Where the “stable value” means a large 
number of malaria cases, the right step seems to be the reduction of the number 
of mosquitoes; on the other hand, where it means a small number of cases it 
should be possible to segregate and isolate these cases. 


In conclusion I wish to express my grateful thanks to Professor Karl Pearson 
for much valuable assistance in the preparation of this paper. 


[ Postscript. Since completing this paper, I have, through the kindness of 
Professor Ross, been able to see proofs of part of his forthcoming treatise on 
Malaria. From a comparison of the two accounts it will be seen that there is 
complete agreement between us on the essential and fundamental points, but some 
difference in numerical details. This difference is due to an attempt on my part 
to treat the question by fuller and more rigid mathematical methods which would 
be out of place in a treatise written from the medical standpoint, while the 
simpler methods employed by Professor Ross give results which are sufficiently 
accurate for practical purposes. 


The principal points of agreement are : (a) for a given number of anophelines 
per unit of the population the number of malaria cases will gradually rise or fall 
to a fixed value at which it will remain stationary, and (b) when the anophelines are 
less than a certain number (about forty per unit of the population) there can be 
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no stable condition and the malaria cases will gradually decrease and finally 
disappear. 


The differences are: (c) the number of cases in the “stable population,” 
corresponding to a given number of anophelines, is higher in Professor Ross’ 
treatment than in mine, and (d) the changes in the malaria rate are more rapid ; 
e.g., in a particular example, where the population is 1000, the number of 
anophelines »cr person 100 and the initial number of cases 500, we have: 


Months...... 0 1 2 3 4... finally 
Cases (Ross) 500 525 544 560 571 600 
» (Waite) 500 512 521 529 535 556 


The divergence seems to be chiefly due to the difference in the time units 
employed in the two methods of treatment. Professor Ross has used the month 
throughout and has taken the value of m constant during each month, while 
I have used the average time between two consecutive infecting bites as my unit. 
The fact that m is increased by unity each time a healthy person is bitten by an 
infected mosquito and is continually being diminished owing to recoveries, fully 
justifies, in my opinion, the adoption of this unit. 


Part of the difference (d) is easily accounted for as follows: Suppose m is 
increasing; then in finding m,,, from m,, the number of new cases obtained by 
Professor Ross is too great; for m, is the smallest value of m during the month, 
and therefore the chance of malarial sufferers being bitten again is greater and 
that of healthy persons being infected less than that found. Again, the number of 
recoveries found for the month is, for the same reason, too small, and thus both 
causes combine to make the value of m,,, too great. 


Similar reasoning shows that when the malaria is decreasing the value found 
for m,4, is too small; ie. in both cases the changes in the malaria rate are too 
rapid. It is evident, however, that the difference on this account becomes smaller 


as the stable value is approached. H. W.] 





























TABLES OF THE TETRACHORIC FUNCTIONS FOR 
FOURFOLD CORRELATION TABLES. 


By P. F. EVERITT, B.Sc. 


Explanation of the Tables. 


THE method of the fourfold table for determining correlation was described by 
Pearson in Phil. Trans. A. vol. 195, pp. 1—47. 











. b a | b a+b 
? k c | d c+d 
amen" wale 
} ate b+d N 








Following his nomenclature, the normal correlation surface is divided into four 
parts by two planes at right angles to the axes of « and y at distances h’, k’ from 
the origin and these distances h’ and k’, when measured in terms of the standard 
deviations of their respective variables, are called h and k. The volumes or 
frequencies in the four divisions are represented by a, 6, c, d in the manner shown 
in the plan and it will be seen, that 6+d and c +d, owing to the position given to 
the point of intersection of the traces of the dividing planes, cannot exceed }N. 
In the paper referred to, the correlation coefficient r is determined from the 
equation 
d b+dc+d 
ei mae 
where H, K are the ordinates of the normal curve of area N corresponding to the 
abscissae h and & and consequently dividing the curve into areas, of which the pro- 


© /pn caw ae 
1 \n! ‘ 


; b+d c+d : a oo ‘aa . 
portions to the whole are W and V respectively, while 0 and @ are given by 
Un “ hn = (n on; 1) Vp-2, Wy = kWn == 1) Wn —2) 

%=1, %4=h, wm=1, WM=k. 
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Now let 
Ht, , KW, 
T™=—— and fr, =, 
"Nn Vn! 
then the equation becomes 
d b+dc+d. 2% ; 
a ae + S(tat 7), 


: : ; c+d ; 
It is clear from the above, that 7,’ is the same function of a as T, is of 


b+d ; ; : ; ‘ 
y and that one table of these functions will serve for both, if we enter the 


d 


os OF os ee i 
table with Vv as argument for the latter, and with 
p 





l : 
nt for the former. It 
v) 
.. b+d c+d ' ; ‘ 
should be noted that these quantities 7 , are identical with 4(1-— a), 
4 4¥ 

where 4 (1+) and a are used as arguments by Sheppard* in his published tables 
and to avoid ambiguity have been tabulated under that heading. 


In the present tables the values of the first six + functions, henceforth to be 
termed tetrachoric functions, have been computed for values of 4(1—a) from 
‘001 to ‘500 by successive increments of ‘001; the last column contains the values 
of h+ (or k) corresponding to the value of $(1— a) given in the first column, 
and required for computing the functions of higher order than the sixth as 
well as for the probable error of +. 

In the auxiliary table are given the values of p, and q, required to compute the 
functions of the seventh to twelfth orders by means of the difference formula 


Tr= hprTr- — QnTn-2- 


Illustrations of the Use of the Tables. 
(a) With interpolation. 


Consider the hypothetical table given below : 


1668 131 1799 





137 64 201 








1805 | 195 2000 





| 
| 


bid | 
Here *¢_.09750, °**—-10050, % =-032000. 
N N 


* W. F. Sheppard, ‘‘ New Tables of the Probability Integral,” Biometrika, Vol. 1. p. 174. 
+ x of Sheppard’s Tables. 
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Now enter the table twice using interpolation, once with 4(1—a) equal to 
‘09750 and a second time with it equal to ‘10050, and we have 


4 (1-2) 7 T2 T3 ™%] 75 ™%; 
09750 +°17228 +°15787 +4°04779 -—-06018 —-06693 +-00854 
10050 +°17614 +°15926 +4°04567 -—-06275 -—-06652 +4-01111 


Multiplying the numbers in each column together, we find the equation 


‘032000 = 009799 + °030345r + 0251427? + -002183r° 
+ °003776r4 + °0044527° + -0000957° ; 
whence, solving by Newton’s method, 


r= 3501. 


An illustration of the calculation of the probable error of r when found as 
above is given in Pearson’s paper on p. 36. 


(6) Without interpolation. 
Using the same table as in example (a), we have as before 


b+d_ 


09750, —— =:'10050, @ _ 939000. 
N 


N N 
Entering the table with $(1 — a) equal to ‘098 and -101 we have 
4 (1-a) Ty T2 T3 T% T5 7; 
098 +°17292 4+°15811 +4°04744 -—-06061 —-06687 +-00897 
‘101 +°17678 +°15948 +4°04531 --06317 -—-06644 +4°01153 
Multiplying out as before 
032000 = 009898 + -030569r + °0252157* + °002150r* 
+ 00382974 + 0044437 + -0001037° ; 
whence, solving by Newton’s method, 
r = °498., 


This particular illustration, chosen so that the true values of $(1—a) fall 
midway between tabulated values, shows the maximum error caused by working 
without interpolation for values of 4(1 —4).of ‘100, and in practically every case 
this error will be negligible when compared with the probable error of r. 


A study of the differences of the functions shows that for values of } (1 — a) 
greater than ‘100 the error introduced by working without interpolation will not 
be appreciably greater than in the example given and may quite easily be less; 
for values of 4 (1 — @) less than *100 it will be desirable to use interpolation if the 
greatest accuracy attainable is desired, but even in very unfavourable cases such 
errors will rarely become as large as the probable error of the result. 


56—2 
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The tetrachoric functions of higher order thu the sixth (easily found by 
means of the auxiliary table) will be chiefly required when r is large and the series 
consequently converges less rapidly. In connection with such cases, it is well to 
recollect that in the series on the right-hand side of the equation successive 


terms need not necessarily converge but that alternate terms are absolutely 
convergent *. 


Another Use for the Tables. 


The tables can also be used with advantage in the work of computing the 
frequencies in the divisions of any frequency table, when the correlation coefficient 
r is known and the frequency is normal. Such problems can be classified into 
three cases according to the fineness of the grouping and a brief outline of the 
methods to be employed is given. 


Case I. Fourfold table. According to the purpose to be served, the 
dividing planes are fixed either by h and & or b+d and c+d being known; 
the value of d is then calculated from 


d b+dc+d % ote 
, Ss ee 


when d is known, b, c and a are also known from the total column and row. 


Case II. Ninefold table. It is required to divide a normal frequency surface 
into nine parts A, B,...£as shown in the diagram. Then it will always be possible 














| A B C | 
p| gz F 
| 
S ewe. L | 
| 
| | | y | 
| | 


to consider one or more of A, C, G, Z in turn as the d of a standard fourfold 
table and find its value from the equation given in Case I. The column and 
row totals are either known or obtained directly from values of h and k. From 
the known divisions and the column and row totals the values of other divisions 
and groups of divisions are next obtained; the table is then divided afresh, if 
necessary, according to the particular circumstances of the case, the values of 
other groups of divisions obtained, and so ultimately the value of each separate 
division not already known may be found by differences. 


* Pearson, loc. cit. p. 10. 
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Case III. More than three columns or rows. 


Find the value of as many of the corner divisions as may be treated as the 
d of a fourfold table by the method of Case I; then find the column and row 
totals, if not already known, and from thein the values of other divisions and 
groups of divisions. The table is then divided afresh and the values of other 
groups of divisions obtained and this process repeated until the values of the 
remaining separate divisions may be found by differences. This case is really an 
extension of Case II. 


In conclusion I desire to express my indebtedness to Prof. Pearson for 
suggesting the construction of the tables and for much valuable advice as to 
methods of calculation and arrangement. I also wish to express my thanks to 
Mr Sheppard for allowing me to use a hitherto unpublished table of the 
probability integral, from which the values of corresponding to given values of 
4(1—a) have been taken, the published tables not providing sufficient decimal 
places without troublesome interpolation. 


Ausiliary Table for calculating Tetrachoric Functions 


beyond those Tabled. 





] 
| n Pu | In | 
| 
| 
| 7 | -37796 | -7152 
} 8 *SOL78 
9 82496 | 

10 "84327 | 
11 30151 "85812 
12 28868 ‘87039 


Ta ANT 1— YnTn-2- 








442 


$(1=a) 


‘001 
“002 
‘003 
O04 
"005 


‘006 
007 
‘008 
“009 
‘O10 


‘Oll 
“012 
‘O13 
“O14 
O15 


‘016 
O17 
‘018 
| O19 
| °020 


‘021 
| *Q22 
"023 
"O24 
"025 


"026 
“O27 
‘028 
"029 
"030 


| 031 
| "032 
| *033 
| O34 
| 


O35 


‘036 
"O37 
‘038 
“089 
“040 


“O41 
‘O42 
‘O43 
Ohh 


O45 


O46 
O47 
048 
O49 
| "050 


| 


Tetrachoric Functions for Fourfold Correlation Tables 


TABLE OF THE TETRACHORIC 


vi 


00337 
‘00634 
“00915 
“01185 
“01446 


01700 
01949 
“02192 
02431 
02665 


“02896 
"03123 
03348 
“03569 
‘03787 


“04003 
04216 
"04427 
04635 
"04842 


05046 
"05249 
“05449 
“05648 
"05845 


“06040 
06233 
“06425 
‘06615 
“06804 


06992 
07177 
07362 
07545 
07727 


07908 
“08087 
“08265 
“08442 
‘08617 


08792 
08965 
09137 
“09309 
‘09479 


09648 
“09816 
09983 
“10149 
“10314 


Te 


00736 
“01290 
01778 
02222 
02634 


-03020 
03386 
03734 
“04066 
04384 


“04690 
“04985 
"05270 
"05545 
“05811 


“06069 
“06320 
“06564 
“06801 
‘07031 


*07256 
‘07475 
‘O7688 
“07897 
“08100 


“08299 
“08493 
“08682 
“O8868 
“09049 


09227 
“09400 
*09570 
‘09737 
“09900 


“10060 
"10216 
*10370 
“10520 
*10668 


“10812 
"10954 
*11093 
*11229 
11363 


"11495 
“11623 
"11750 
“11874 
“11996 





++tet t4+¢4+ 


t+++++ 


EERE FEEEHE FHEEHHE FETE Ftttte Ftettt 


+++++ 


T3 


01175 
01885 
“02446 
02918 
03326 


‘03686 
“04008 
“04298 
“04561 
“04800 


“05020 
“05221 
“05406 
‘05577 
05735 


“05880 
‘O6015 
06139 
“0254 
‘06361 


‘06459 
“06549 
‘06633 
‘06709 
*06780 


‘06844 
*06903 
"06956 
‘07005 
07048 


07087 
07122 
07153 
‘07179 
07202 


07221 
*07237 
“07249 


°07258 


‘07264 


‘07268 
07268 
‘07266 
07261 
‘07253 


07243 
07231 
‘07217 
‘07200 
07181 


4 


+ °01391 
+ 01968 
*02335 
*02587 
“02764 


“02887 
“02970 
“03021 
‘03047 
03053 


“03041 
“03014 
“02975 
“02926 
‘02867 


“02801 
02727 
“02647 
“02562 
02472 


‘02378 
“02280 


HEHE FEEHHE FHEttH FHttte +++ 


‘01634 
“01520 
“01404 


+ 

oo 

+ 

“fb 

+ °01287 
+°01168 
+°01049 
+ 00929 
+ ‘00809 
+ ‘00688 
+ ‘00566 
+ 00444 
+ 00322 
+ *00200 
-f 


*00077 
00045 
00167 
“00290 
“00412 


— 00534 
— ‘00656 
— ‘00778 
— ‘00899 
— ‘01020 


I++++ t+4+4+4+4+ 


FUNCTIONS. 
7} 7; 
“01134 + 00415 
01269 + °00053 
“01228 ~ ‘00328 
‘OLLI — ‘00687 
“00952 — ‘01017 
00770 — 01318 
‘00575 - O1592 
‘0037 1 — ‘01841 
00164 | — 02067 
00044 | —:02271 
00253 | — 02457 
00460 | — ‘02625 
- 00664 | — 02777 
00866 | —-02914 
01064 | — :03037 
| 
| 

01259 | -°03147 
“01449 — 03246 
01636 — 03334 
01818 — ‘03411 
01996 — ‘03479 
‘02170 | —-:03538 
02340 | —-03589 
02505 | —-03632 
02666 | -— 03667 
02823 | — ‘03696 


‘02976 
03125 
“03270 
‘03411 
03547 


‘03680 
03810 
03935 
“04057 
‘04176 


04291 
“04402 
‘04510 
‘04615 
‘04717 


04815 
04910 
*05003 
“05092 
- 05178 


- ‘05261 


"05342 
“05420 


- "05495 


2 and h are essentially positive. 


*05567 


‘03718 
03734 
03744 
*03749 
*03749 


"03744 
‘03734 
‘03720 
‘03702 
‘03680 


03654 
“03624 
"03592 
"03556 
- ‘03517 


03475 
*03431 
“03384 
03335 
03283 


*03229 
03173 
03115 
“03055 
“02994 


| 
| 
| 


3°09023 
2°87816 
‘74778 
*65207 
*D7583 


bo bo bo 


“51214 
*45726 
“40892 
*36562 
*32635 


bo bo bo & bo 


*29037 
*25713 
°22621 
19729 
17009 


bo bo bo bb 


14441 
12007 
‘09693 
“07485 
‘05375 


bo bo bo bo bt 


03352 
“01409 
"99539 
‘97737 
“95996 


— — = bS bD 


"94313 
“92684 
91104 
*89570 
*88079 


— ee 


“86630 
"85218 
"83842 
*82501 
“81191 


“79912 
‘78661 
‘77438 
‘76241 
*75069 


— ee eet 


*73920 
°72793 
‘71689 
*70604 
“69540 


— et et 


"68494 
67466 
"66456 
*65463 
164485 


—— i 
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TABLE—(continued). 


| | | 














140-3 9 " — )- ee oo ee aoe “at 
a oa ae | | 
| 
f 051 | °10478 | -12115 +°07160 | —°01140 — *05637 — 02931 1°63523 | 
| 052 10641 | °12232 | +°07138 | —-01260 — °05704 — 02866 1°62576 | 
| 058 "10803. | +12347 +°07113 | — 01380 — 05769 — ‘02799 1°61644 
"054 “10964 "12460 | +°07087 | -°01499 — 05831 — *02732 1°60725 | 
055 “11124 | +12571 | +4+-°07058 | — -01618 — 05891 — 02662 159819 | 
| } | 
056 "11284 | 12680 | +°07028 | — ‘01736 — *05949 — 02592 1°58927 | 
057 11442 | -12787 | +°06997 | -—-°01854 — 06004 — °02520 1°58047 
058 11600 | -12892 | +°06964 | —-01971 — 06057 — 02447 157179 
059 11756 12995 | +°06929 | — °02087 — ‘06107 — °02373 1°56322 
‘060 11912 13096 | +-06893 — 02203 — 06155 — °02298 1°55477 | 
061 “12067 13196 | + °06855 — 02318 — 06202 — ‘02222 1°54643 
| 062 13293 | -+ 06816 — 02433 — 06246 — °02145 1°53820 
"063 "13389 | + 06775 — °02547 — *06288 — 02068 1°53007 | 
064 | °13483 | -+°06734 — 02660 — 06328 — 01989 1752204 | 
065 12679 | 13575 | +°06690 02773 — 06365 — ‘01910 1°51410 | 
066 12830 13666 | +4+°06646 | —-02884 —°06401 | —°01830 1°50626 
067 "12981 | -13754 | +°06601 | — ‘02996 — 06435 — ‘01749 1°49851 
068 13130 | +13842 | +°06554 — ‘03106 — 06467 | — 01668 1°49085 
| 069 13279 | :13927 | +°06506 | —-°03216 — ‘06498 — 01586 1°48328 
| ‘070 *13427 | #°14011 | +°06457 | --°03325 — 06526 | - 01504 1°47579 
| | 
| “O71 *13574 =| 14094 +°06407 | —°03433 — 06552 — "01421 1°46838 
O72 "13720 14175 +°06356 | —°03541 = 06577 | —°01337 1°46106 
073 "13866 14254 +°06304 | —-03648 — 06600 — 01253 1°45381 
OT 4 "14011 “14332 +°06251 | —-03754 — 06621 — 01169 1°44663 
O75 14156 14409 + ‘06197 | — 03859 | —-06641 | — 01085 1°43953 
O76 "14299 "14484 +°06142 | — 03963 | — ‘06659 “01000 1°43250 
O77 14442 | °14558 +°06086 | —-04067 —°06675 | —-00915 1°42554 
078 14584 | -14630 | +-06029 | --04170 | —-06690 | — -00829 1°41865 
O79 “14726 | °14701 + °05971 — 04272 — 06703 — 00743 1°41183 
‘080 14867 | - :14771 +°05913 | --043874 | —-06715 — 00658 1°40507 
"081 “15007 | 14839 | + °05854 | — 04474 — ‘06725 — 00572 1°39838 
082 “15146 | = -14906 + 05794 — 04574 — 06733 — 00485 1°39174 
‘083 "15285 ‘14971 + °05733 — 04673 — ‘06741 — ‘00399 1°38517 
O84 "15423 | =*15036 + °05671 — ‘04771 — ‘06746 — 00312 1-37866 
085 "15561 | -+15099 + °05609 — 04869 — ‘06751 — 00226 1°37220 
086 15698 | -+15160 +°05546 | — 04965 — 06753 —°00139 | 1°36581 
O87 "15834 | °15221 + 05483 = 05061 — 06755 — 00053 1°35946 
088 “15970 | *15280 + °05418 — 05156 — 06755 + °00034 1°35317 
‘089 ‘16105 | 15339 +°05353 | — °05250 — 06754 | +°00120 1°34694 | 
"090 "16239 | 15396 +°05288 | —°05344 — 06751 (es 00207 1°34076 | 
091 "16373 | °15451 +°05222 | — 05436 — ‘06748 +°00294 1°33462 
092 16506 | -15506 +.°05155 — 05528 — 06743 + "00380 1°32854 
"098 16639 | 15560 + °05088 — *05619 — 06736 + °00467 1°32251 
O94 "16770 | -15612 + *05020 — ‘05709 — 06729 + °00553 1°31652 
095 "16902 | *15663 + °04952 — ‘05798 — ‘06720 + °00639 1°31058 
| 
| 096 *17033 “15713 + 04883 — 05887 — ‘06710 + °00725 1°30469 
097 17163 15763 +°04813 — 05975 — 06699 +°00811 1:29884 
098 *17292 ‘15811 + 04744 — 06061 — 06687 + 00897 1°29303 
099 17421 *15858 + 04673 — 06148 — ‘06674 + 00982 1°28727 


100 ‘17550 | °15904 + 04602 — ‘06233 - ‘06660 + 01068 1°28155 
ne ee ees eee | j 
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*22473 
"22580 
“22686 


"22792 


TABLE—(continued). 


! 





"17025 
*17035 
“17044 


“O1356 
‘O1l275 


*O1194 


T1 | T2 T3 
‘17678 "15948 + °04531 
*17805 "15992 + 04459 
*17932 “16035 + °04387 
*18058 “16077 + 04315 
“18184 "16118 + 04242 
*18309 “16158 + °04169 
"18433 “16197 +°04095 
*18557 "16235 + 04021 
“18681 "16272 + °03947 
“18804 *16308 + °03872 
"18926 "16343 + °03797 
“19048 “16378 + ‘03721 
“19169 ‘16411 + °03646 
*19290 "16443 + 03570 
“19410 "16475 + °03493 
*19530 “16506 + °03417 
"19649 *16536 + 03340 
“19768 "16565 + 03263 
“19886 *16593 + ‘03186 
“20004 “16620 + °03108 
“20121 “16647 + °03030 
*20238 “16672 + °02952 
*20354 *16697 + ‘02874 
*20470 16721 + °02796 
"20585 "16745 +°02717 
*20700 "16767 + 02638 
*20814 "16789 + °02559 
‘20928 “16810 + *02480 
‘21042 *16830 + °02401 
*21155 “16849 + °02321 
*21267 "16868 + 02241 
‘21379 “16886 + °02162 
*21490 “16903 + ‘02082 
‘21601 "16919 + 02001 
*21712 “16935 +°01921 

21822 *16950 + 01841 
21932 “16964 + °01760 
22041 ‘16978 + °01680 
22149 “16990 + °O01599 
22258 17003 +°01518 
22365 17014 + 01437 

+ 

o 

a 

+ 


‘22898 
"23003 
‘23108 
*23212 


23316 


*17053 


‘17061 
*17069 
“17076 
‘17082 
*1L7O87 


‘O1113 


“01032 
“00950 
“OO869 
“OO788 
+ ‘00706 


++++ 








— ‘09307 _ 
- ‘09356 


— 09405 


“05078 
“05028 
- ‘04978 
- ‘09452 - 
— ‘09499 — 


“04928 








“04612 
‘04678 
“O4744 
*O4809 


“05374 
"04939 
04505 
‘04073 


T% TS | T%; | h 
— 06317 — ‘06644 +°01153 | 1°27587 
— 06401 — ‘06628 +°01238 | 1:27024 
— ‘06484 —- ‘06610 +°01322 | 1°26464 
— 06566 — 06592 +°01407 | 1°25908 
-- 06647 —*06572 | +°01491 | 1:25357 
} | 

— 06727 — ‘06551 | + °O1575 1°24808 
— 06807 —°06530 | +:°01659 1°24264 
— ‘06886 — *06507 +°01742 1°23723 
— 06964 — "06484 | +°01825 123186 
— ‘07041 — ‘06459 +°01908 | 1°22653 
— ‘07117 — ‘06434 + 01990 1°22123 
— ‘07193 — °06408 + 02072 1°21596 
— ‘07268 — 06381 + 02154 1°21073 
— ‘07342 — °06353 + °02235 1°20553 
—°07415 | — 06324 + °02316 1°20036 
— 07488 — 06294 +°02397 | 1°19522 
— ‘07559 — 06264 | +°02477 | 1°19012 
— ‘07630 — 06233 + °02557 1°18504 
— ‘07700 | —-06201 + °02636 1°18000 
— '07770 — °06168 + °02716 1°17499 
— ‘07838 — ‘06134 + 02794 1°17000 
— 07906 | —-06100 + °02873 1°16505 
— 07973 | -—-06065 | +-°02950 1°16012 
—°08039 | — -06029 | + 03028 1°15522 
— ‘08105 | — "05992 | +°03105 1°15035 
—°08169 | —-05955 | + ‘03181 1°14551 
— ‘08233 — ‘05917 + '03257 1°14069 
— ‘08297 — ‘05878 + °03333 1°13590 
— ‘08359 -- 05839 + °03408 1°13113 
— 08421 — ‘05799 + 03483 1°12639 
— ‘08482 | —-°05758 + °03557 1°12168 
— ‘08542 — ‘05717 + 03631 111699 
~ ‘08601 — ‘05675 + 03704 1°11232 
— ‘08660 — 05632 + 03777 1°10768 
-~ -O8718 — 05589 + ‘03850 1°10306 
— 08775 — 05546 + °03921 1°09847 
— ‘08831 — 05501 + °03993 1°09390 
— ‘08887 — 05456 + 04064 1°08935 
— ‘08942 — ‘05411 + 04134 1°08482 
— ‘08996 — *05365 + °04204 1°08032 
~ 09050 — ‘05318 + °04273 1°07584 
- ‘09103 — ‘05271 + *04342 1°07138 
— ‘09155 — 05224 + ‘04410 106694 
— ‘09206 — ‘05176 + ‘04478 1°06252 
- ‘09257 — 05127 + 04545 1°05812 

a 

+ 

a 

+ 

u 


‘O4877 





“OA87T4 


‘03643 























P. F. Everirr 445 





TABLE—(continued). 



































| | $(1—a) 77 T2 Ts 7, 7 | 7; | h 
| — : | | 

| *151 *23419 *17092 +°00625 | — 09546 | — 04825 + 04938 1°03215 
°152 *23522 “17097 + 00543 —°09592 | —-°04774 +°05002 | 1°02789 
*153 *23625 *17100 + °00462 — ‘09637 — 04721 + °05065 1°02365 
"154 *23727 ‘17103 + 00380 — ‘09681 — 04669 +°05127 | 1°01943 
°155 "23829 ‘17106 | +°00298 — ‘09725 — ‘04615 + °05189 | 1°01522 
156 *23930 *17108 + 00217 — ‘09768 — 04562 +°05250 | 1:°01103 
157 24031 17109 | +-00135 — 09810 — 04508 +°05311 | 1:00686 
+158 24131 17110 | +-°00053 — ‘09852 — 04454 +°05371 | 1:00271 
159 24232 ‘17110 | -- -00028 ~ ‘09892 — 04399 +°05431 | 99858 
160 *24331 17109 | —-00110 - 09933 04344 | +°05490 | 99446 
“161 "24430 *17108 | — ‘00191 — 09972 — "04288 | +:°05549 | “99036 
162 "24529 17107 ce °00273. | —-10011 — 04232 + °05607 *98627 
163 *24628 17104. | —-00355 — *10049 — 04176 +°05664 | -98220 
“164 *24726 17102 | —-00436 — *10087 — 04120 + °05721 | ‘97815 
"165 *24823 *17098 | — ‘00518 — 10124 — 04063 + °05778 | ‘97411 
166 *24921 *17094 — ‘00599 — ‘10160 — *04006 + °05834 ‘97009 
167 *25017 *17090 — 00681 — ‘10196 — 03948 +°05889 | 96609 
168 *25114 “17085 — 00762 — *10231 — ‘03890 +°05943 | -96210 
169 ‘25210 | -°17080 — 00844 — ‘10265 — 03832 +°05998 | 95812 
“170 °25305 | °17073 — 00925 — *10299 — 03774 + ‘06051 | *95417 
171 25401 | ‘17067 | —-01007 — *10332 —°03715 | +°06104 | -95022 
172 *25495 | 17060 | — 01088 — 10364 - 03656 + 06156 *94629 
173 *25590 "17052 | —-01169 — 10396 — *03597 + °06208 "94238 
"174 -25684 | -°17044 Ly 01251 — ‘10427 — 03537 + *06260 ‘93848 
*175 *25778 | *17035 — 01332 — °10458 — *03478 + °06310 "93459 
176 | -25871 | -17026 | —-o1413 | —-10487 | —-08417 | +-06360 | 93072 
ev 25964 | ‘17016 | —-01494 — °10517 — ‘03357 + °06410 “92686 
178 *26056 | *17006 — 01575 — ‘10545 — ‘03296 + °06459 “92301 
179 *26148 "16995 — 01656 —°10573 — 03236 + °06507 91918 
180 "26240 "16984 — ‘01737 — ‘10601 — 03175 + °06555 "91537 
"181 26331 | -16972 | —-01817 — ‘10627 — ‘03113 + 06603 ‘91156 
"182 26422 | -16960 | —-01898 — *10653 — 03052 + °06649 ‘90777 
183 26513 | °16948 | --01979 —+10679 — °02990 + *06695 -90399 
184 26603 16934 — 02059 — ‘10704 — 02928 + °06741 -90023 
"185 *26693 "16921 | —-02146 — ‘10728 — ‘02866 + °06786 *89647 
"186 26782 | -16907 | —-02220 — *10752 — 02803 + °06830 "89273 
187 "26871 | °16892 | —-02300 — °10775 — 02741 + °06874 “88901 
"188 *26960 *16877 — *02380 — ‘10798 — ‘02678 + °06917 "88529 
189 27049 | -16861 | —-02460 —+10819 — 02615 + °06960 "88159 
-190 27137 | -16845 | —-02540 — 10841 — 02552 +°07002 | -87790 

| | 
191 27224 | -16829 - 02620 — ‘10861 — 02489 + "07044 87422 
"192 ‘27311 | +16812 | —-02700 — *10882 — 02425 + 07085 “87055 
+193 27398 "16795 — 02779 - *10901 — 02362 + °07125 “86689 
194 27485 "16777 — 02859 — 10920 — *02298 + 07165 86325 
195 27571 16759 | —-02938 — ‘10939 -- 02234 +°07204 | 85962 
196 ‘27657 | +16740 | —-03018 — *10956 — ‘02170 + °07243 *85600 
197 °27742 | #*16721 | —-03097 — ‘10974 ‘02106 + °07281 "85239 
198 ‘27827 | °16701 | —°03176 “10990 02041 + 07319 *84879 
“199 ‘27912 | -16681 | — ‘03255 — ‘11007 — ‘01977 + °07356 "84520 
*200 27996 | ‘16661 | — 03334 — *11022 — 01912 + °07392 84162 | 
Biometrika vir 57 
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*201 ‘28080 “16640 
"202 28164 "16619 
"205 *28247 "16597 
"204 | °*28330 *16575 
*205 28413 *16553 


| +206 *28495 *16530 
"207 *28577 *16506 
*208 "28658 “16483 
209 | 28739 ‘16459 
“210 | *28820 *,6434 
211 | 28901 | -16409 
212 | -28981 | -16384 

213 | +29060 | -16358 
214 | 29140 | -16332 

| 215 | 29219 "16305 
"216 *29298 *16279 
"O17 *29376 *16251 
218 *29454 “16224 
‘219 | +29532 | -16196 
‘220 | +29609 ‘16167 

| 
*221 "29686 “16139 
“222 *29763 “16110 
223 *29840 “16080 
"224 ‘29916 "16050 
"225 *29991 *16020 
-29G *30067 *15990 
"227 30142 "15959 
‘228 *30216 *15927 
*229 *30291 *15896 
+230 30365 15864 
*231 “30439 "15832 
*932 *B0512 “15799 
+233 *B0585 *15766 
23) 30658 *15733 
"O35 *B0730 *15699 
"236 *30802 *15665 
+23} "30874 ‘15631 
"238 *30945 *15596 
+239 ‘31017 *15561 
‘240 ‘31087 "15526 
241 ‘31158 “15490 
"242 *31228 "15454 
243 *31298 *15418 
Dhl, ‘31367 "15382 
"24D "31436 *15345 
"246 *31505 "15308 

| "247 ‘31574 ‘15270 
"248 31642 *15232 
"249 “31710 *15194 
*250 “31778 “15156 


TABLE— (continued). 


— 03569 
— 03648 
— ‘03726 


— ‘03804 
— ‘03882 
— 03959 
— °04037 
— ‘04114 


— ‘04192 
— ‘04269 
— ‘04346 
— 04423 
— ‘04499 


— ‘04576 
— ‘04652 
— ‘04728 
— 04804 
— ‘04880 


— "04956 
— ‘05031 
~— ‘05107 
— *05182 


— ‘05257 


— *05332 
— ‘05406 
— ‘05481 
— "05555 
— *05629 


=— 05703 
- 05777 
- "O5851 

— 05924 
— ‘05997 


— ‘06070 
— ‘06143 
— 06215 
— ‘06288 
— ‘06360 


— 06432 
— 06504 
‘06576 
— ‘06647 
— ‘06718 


— ‘06789 
— ‘06860 
— ‘06931 
— 07001 
— ‘07071 





™%} 7% 7; 
*11037 — ‘01848 + °07428 
"11051 — ‘01783 + ‘07464 
*11065 — ‘01718 + 07498 

- 11079 — ‘01653 + 07532 
‘11091 — 01587 + *07566 
"11104 — ‘01522 + 07599 
“11115 — ‘01457 + °07632 
"11126 — ‘01391 + °07664 
‘11137 — 01326 + 07695 
"11147 — ‘01260 + ‘07726 
“11157 — 01194 + ‘07756 
“11166 — °01129 + °07786 
11174 — ‘01063 + °O7815 
11182 ~ 00997 407844 
“11189 — 00931 + ‘07872 
“11196 — ‘00865 + ‘07899 
*11203 — ‘00799 + 07926 
*11208 — ‘00733 + ‘07952 
11214 — ‘00667 + ‘07978 
“11218 — ‘00600 + °08004 
11223 — 00534 + °08028 
*11226 — ‘00468 + °08052 
*11230 — ‘00402 + ‘08076 

11233 — ‘00335 + ‘08099 

11235 — 00269 + 08122 

11237 — ‘00203 + 08144 

11238 — ‘00136 + 08165 

11239 — ‘00070 +°08186 

11239 - 00004 + °08207 

11239 + ‘00063 + *08226 

11238 + 00129 + *08246 

11237 + °00195 + ‘08265 

11235 + 00262 + °08283 

11233 + ‘00328 + °08301 

11230 + 00394 + 08318 
*11227 + ‘00461 + °08334 
*11224 + *00527 + 08351 
*11220 + °00593 + ‘08366 
"11215 + 00659 +°08381 
‘11210 + °00726 + °08396 
*11205 + ‘00792 + 08410 
“11199 + °00858 + °08423 

-°11192 + °00924 + 08436 
"11185 + *00990 + '08449 
‘11178 + °01056 + ‘08461 
‘11170 + O115? + ‘08472 
“11162 +°01188 + 08483 

—*11154 + °01253 + "08494 
“11145 +°O1319 + °08504 
"11135 +°01385 + °08513 


*83805 


83450 


*83095 
*82742 
*82389 
*82038 
*81687 
*81338 


“80990 
"80642 


*80296 
*79950 
*79606 
°79262 


‘78919 


*78577 


*78237 


*77897 


‘77557 
‘77219 


‘76882 
“76546 
‘76210 
“75875 


*75541 


“75208 
‘74876 





*70952 


‘70630 


*70309 
“69988 
‘69668 
“69349 
“69031 


68713 


68396 
“68080 


"67764 
67449 
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TABLE—(continued). 

















| | | | 
| $(1-a) T% T2 | T3 | T% | TS 7% h 
"251 31845 | °15117 — ‘07141 — °11125 + °01450 + °08522 *67135 
"252 *31912 *15078 — ‘07211 — "11115 + 01516 +°08530 | °66821 
D253 °31979 | -15039 —°07280 | -—-:11104 + 01581 + °08538 “66508 
"24 "32045 *14999 - 07350 — *11093 + ‘01647 + 08546 “66196 
"255 “32111 “14959 — ‘07419 —°11081 + °O1712 + °08553 “65884 
"256 “14919 — ‘07488 —°11069 + °O1777 + °08559 *65573 
257 | 14879 - 07557 | —-11056 + °01842 + °08565 "65262 
258 "14838 — 07625 | —-+11043 + 01907 + °08571 "64952 
259 "14797 — ‘07693 — °*11030 + 01972 + 08575 64643 | 
260 "14756 — ‘07761 — ‘11016 + °02037 + °08580 64335 | 
261 32501 14714 ~ 07829 —°11002 + °02102 + 08584 64027 | 
"262 *32565 “14672 — ‘07897 — *10987 + °02166 + *08587 *63719 
263 32628 *14630 — ‘07964 — ‘10972 + °02231 +°08590 “63412 | 
264 ‘32691 14588 | —-08031 — "10956 | +°02295 | +-08593 | -63106 | 
265 | 32754 "14545 —'08098 | —:10940 | +°02360 | +°08595 | 62801 
| | 
266 32817 "14502 — 08165 — 10924 + 02424 +°08596 | °62496 
°267 | °32879 "14459 — ‘08231 — *10907 + °02488 +°08597 | °62191 
"268 *32941 "14415 — "08298 — ‘10890 + 02552 +°08598 | °61887 
*269 *33003 *14372 — ‘08364 — °10873 + 02616 +°08598 | “61584 
‘270 | °33065 "14328 — “08429 — *10855 + *02680 +°08598 | “61281 | 
‘271 | *33126 "14283 — ‘08495 — °10837 + °02743 + °08597 | “60979 
“272 *33187 "14239 — ‘08560 — ‘10818 + °02807 + ‘08596 “60678 
; "273 *33247 "14194 — ‘08625 — 10799 + ‘02870 + 08594 “60376 
1 ‘274 =| «=°33307 “14149 — ‘08690 — ‘10779 + °02933 + 08591 “60076 
"275 *33367 “14104 — ‘O8755 “10759 + °02997 + 08589 “D9776 
‘276 *33427 "14058 “O8819 — °10739 + °03060 + ‘08586 “59477 
‘Q77 "33486 "14012 — ‘08883 —°10718 + 03122 + °08582 *d9178 
"278 33545 *13966 — ‘08947 — ‘10697 + °03185 + °08578 “58879 
“279 33604 °13920 — ‘09011 —°10676 + °03248 + °08573 "58581 
*280 *33662 *13873 — ‘09074 — 10654 + °03310 + ‘08568 “58284 
*281 *33720 "13826 — ‘09137 — °10632 + °03372 +°08563 | "57987 
“282 33778 13779 — °09200 — *10609 + °03434 + °08557 “57691 
*283 *33836 *13732 — 09263 — *10587 + °03496 + *08551 "37395 
"284 "33893 *13685 — 09325 — *10563 + °03558 + *08544 “57100 
"285 *33950 *13637 — ‘09388 — *10540 + °03620 +°08536 | °56805 
286 *34007 *13589 — ‘09450 ~-10516 | + °03681 + 08529 56511 
es; *34063 *13541 — ‘09511 — 10491 + °03743 +°08521 | 56217 
"288 “34119 *13492 — 09573 — *10466 + °03804 + 08512 "55924 
289 34175 "13443 — ‘09634 —*10441 + °03865 +°08503 | °55631 
-290 "34230 *13394 — ‘09695 — ‘10416 + °03926 +°08494 | *55338 
"291 "34286 *13345 — ‘09756 —'10390 | + °03987 + ‘08484 | *55047 
292 34341 13296 — ‘09816 — 10364 + *04047 + °08473 54755 
"293 "34395 *13246 — ‘09876 — *10337 +°04107 + °08463 54464 
294 "34449 "13196 — ‘09936 — *10310 + °04168 + ‘08451 54174 | 
*295 *34503 "13146 — ‘09996 — 10283 + 04228 + °08440 "53884 
*296 *B4557 *13096 — ‘10056 — *10256 +°04287 + 08428 53594 
*297 “34611 *13046 -~ "10115 — *10228 + °04347 +°08415 *53305 
298 “34664 "12995 —°10174 —*10199 + °04407 + 08402 53016 
299 *B4717 "12944 *10233 — 10171 + 04466 + °08389 “52728 
300 34769 "12893 — ‘10291 .% "10142 + "04525 + 08375 52440 


57—2 
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"301 *34822 
"3802 "34874 
“303 "34925 











“304 *34977 
"B05 *35028 
306 
307 
308 
309 
310 
‘311 
"312 
"313 
“B14 
"B15 
‘316 | *35572 
‘317 | °35620 
"318 *35667 
“319 | °35714 
320 | -35761 
‘321 | +35808 
"322 | +35854 
“322 | +35900 
“B24 35946 
"325 *BD991 
"326 *36037 
S27 | °36082 
+328 36126 
$29 | °36171 
*330 *36215 
“331 | °*36259 
"332 *36302 
333 *36346 
| "S84 36389 
} 335 *36431 
‘336 | 36474 
"337 "36516 
"338 *36558 
| "339 “36600 
| +340 "36641 
B41 “36682 
342 36723 
343 “36764 
“Bh “36804 
B45 


~ 


"B46 “36884 
"B47 *36923 
348 *36962 
B49 “37001 
350 "37040 


| 
| ‘36844 
| 
| 
| 





"12686 
"12634 


"12581 
*12529 


“11660 
"11604 
"11548 


"11492 
“11436 
"11379 
"11322 
*11265 


“11208 
“11151 
*11093 
“11036 
"10978 


“10920 
*10862 
"10804 
*10745 
‘10687 


‘10628 
"10569 
*10510 
*10451 
“10391 


"10332 
*10272 
"10212 
“10152 
“10092 


TABLE—(continued). 





Correlation Tables 





| 


T3 | 


— '10349 
— ‘10407 
— 10465 
— *10522 
— *10580 


— *10636 
— ‘10693 
— ‘10750 
— 10806 
- "10862 


—°10917 
— 10973 
— ‘11028 
— °11082 
—°11137 


— ‘11191 
— 11245 
— °11299 
— 11353 
— ‘11406 


— 11459 


11512 


— ‘11564 
— 11616 
— ‘11668 


— ‘11720 
—°11771 


- °11822 


— "11873 
— 11923 


— ‘11974 
— °12024 


— "12461 
— 12509 


— °12555 


"12602 


— 12649 


— 12695 
— 12741 
— 12786 
— "12831 
— 12876 


T%} 


‘10113 
“10083 
*10053 
“10023 
“09992 


09961 
09930 
“09889 
“09867 
‘09834 


*09802 
‘09769 
09736 
‘09703 
09669 


“09635 
09600 
‘09566 
09531 
09495 


“09460 
09424 
‘09388 
‘09351 
‘09315 


‘09278 
“09240 
“09203 
‘09165 
09127 


“09088 
09049 
‘09010 
‘0897 | 
‘08932 


“08892 
*O8852 
‘O8811 
‘08771 
‘08730 


“08689 
‘08647 
“08606 
*O8564 
08522 


08479 
08437 
‘08394 
“08351 
08307 











! 
| Ts | T% h 
| 
| a 4 
| + 04584 + 08361 52153 
+ 04643 + °08347 51866 
+ 04701 + 08332 51579 
+ 04759 + 08316 51293 
+°04818 + 08301 “51007 
+ 04876 + 08284 50722 
+ °04933 +°08268 | °50437 
+ 04991 +°08251 | -50153 
| +°05048 | + 08233 49869 
| +°05105 | + 08216 49585 
+ 05162 | + °08197 “49302 
+ °05219 + 08179 “49019 
| +-05276 | +-08160 ‘48736 
| + +05332 +°08140 | 48454 
| + -05388 + 08121 48173 
| + 05444 + 08101 47891 
+ °05500 + ‘08080 47610 
+ °05555 + ‘08059 "47330 
+ 05610 + °08038 47050 
| +°05665 | +-08016 “46770 
| 
+ °05720 +°07994 | °46490 
+°05775 +°O7972 | °46211 
+ 05829 +°07949 | -45933 
+°05883 | +:07926 | -45654 
+ °05937 + °07902 | "45376 
+ 05991 +°07878 | -45099 
+ 06044 +°07854 | +4482] 
+ ‘06097 + 07829 44544 
+ ‘06150 +°07804 | *44268 
+ °06203 +:07779 | -43991 
| 
+ 06255 +°07753 | °43715 
+ °06308 +°07727 | °43440 
+ °06360 + 07701 43164 
+ 06412 + 07674 | +42889 
+ 06463 +°07647 | -42615 
| 
+ 06514 +°07620 | 42340 
+ 06565 +°07592 | -42066 
+ 06616 +°07564 41793 
+ 06667 + °07535 41519 
+°06717 | +°07507 | -41246 
| | 
+°06767 | +:°07477 | -40974 
+ 06817 +°07448 40701 
+ ‘06867 +°07418 | -40429 
| +-06916 +°07388 | +40157 
| +°06965 | +-07358 39886 
| | 
| +-07014 | +-07397 | -39614 
| +:07062 | +-°07296 39343 
+ 07110 + ‘07264 39073 
| +°07158 + 07232 “38802 
| +:07206 + °07200 "38532 
| sa: LP) SER 
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TABLE—(continued). 



































3 (1 —a) v1 | T2 | ° 73 ™] 75 7 h 
ia pete : Cae Se sie 
351 | °37078 | -10032 | —-12921 | —-Ua264 +°07254 | +°07168 | *38262 
B52 “37116 09971 — °12966 — -08220 +°07301 +°07135 | -37993 
353 “37154 09911 —+13010 - 08176 + 07348 +°07102 | 37723 
B54 37192 09850 — *13054 — 08131 + °07395 +°07069 | 37454 
"355 “37229 09789 ~°13097 | —-08087 +°07441 + 07035 37186 
356 "37266 09728 — *13140 — 08042 + °07487 +°07002 | “36917 
357 ‘37303 09667 — 13183 — -07997 + 07533 +°06967 | °36649 
358 37340 09606 — *13226 — 07952 +°07579 | +°06933 | “36381 
"359 ‘37376 09544 | —-13269 — 07906 | +:07624 | + 06898 36113 
“360 37412 09483 ~ 13311 — 07861 | +°07669 | +°06863 "35846 
| ’ 
“361 ‘37447 09421 --°13353 ~-07815 | +-07714 | +°06827 ‘35579 
362 37483 ‘09359 — 13394 — 07768 +°07758 + 06792 35312 | 
"363 37518 -09297 — *13436 ~ 07722 + 07803 + 06756 "35045 | 
“364 "37553 09235 — 13477 — ‘07675 + 07847 | + 06719 34779 
365 37588 09173 | -—°13517 — ‘07629 + 07890 | + 06683 34513 | 
| 
"366 37622 | 09111 | -—°13558 | —-07582 +°07934 | +:°06646 34247 | 
“367 "37656 09048 | — +13598 — 07534 +°07977 | + °06609 33981 | 
368 | ‘37690 | 08985 | —-13638 — 07487 +°08020 | +°06571 ‘33715 
369 | ‘37724 | -08923 | —-13677 — 07439 +°08062 | +°06534 33450 
| 370 | *37757 08860 | — 13717 — 07391 | +-08105 | + 06496 33185 | 
| 371 | +37790 08797 | —°13756 — 07343 | +°08147 | +°06458 "32921 
‘372 | 37823 08734 | —-°13794 —-07295 | +°08188 | +°06419 32656 | 
3738 | °37855 08671 | — 13833 ~-07246 | +:08230 | +°06380 “32392 | 
‘374 | °37888 | -08607 —°13871 —°07198 | +°08271 | +°06341 | °32128 | 
BI5 | ‘37920 | 08544 — *13909 — 07149 +°08312 + 06302 | 31864 | 
| “376 | 37951 | -08480 | —-13946 | —-07100 | +-08352 | +-06262 | "31600 
‘377 | °37983 | -08416 —'13984 | —-07050 + 08392 +°06222 | °31337 
“378 | 38014 | -08353 —-14021 | —-07001 + 08432 +:06182 | °31074 
379 | °38045 | -08289 —'14057 | —-06951 + 08472 +°06142 ‘30811 
| 380 | 38076 | -08225 — "14094 | —-06901 + 08512 +°06101 “30548 
‘381 | °38106 | -08160 — 14130 | — 06851 + 08551 + 06061 “30286 
382 38136 | -08096 —-14166 | —-06801 +°08589 + 06019 30023 
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fk. INTRODUCTION. 


Volumes have been written concerning Natural Selection—its existence or non- 
existence, its potency or its inefficiency—but much to the discredit of biologists, 
thousands of pages of this literature may be explored without finding a single 
oasis of quantitative data to relieve the monotony of the polemic desert. 


Fortunately there are a few notable exceptions, among which memoirs by 
Weldon* +, Bumpus{, Beeton and” Pearson§, Weldon||, Di Cesnola™, Crampton **, 
Di Cesnolat+, Kellogg and Belltt and perhaps a few others may be mentioned 


as contributing to the direct quantitative evidence concerning the potency of 
natural selection. 


These researches, highly valuable in themselves, are only a slight beginning 
in the scientific investigation of a field of cardinal biological importance. Perhaps 
in the earlier days of Darwinism it was not possible to do much more than collect 
qualitative evidence and to reason from analogy; one great value of the researches 
just cited is that they show that the exceedingly complex problems of natural 
selection can be made the field of quantitative investigation, and now that this 
fact has been demonstrated I think that we should stigmatize as merely pseudo- 
scientific discussions of selection not based on quantitative data. 


Some years ago it occurred to me that it might be possible to ascertain some 
facts concerning selective elimination by comparing the characters of ovaries which 
develop into mature fruits with those which do not. That this is a distinct 
problem from that of the selective elimination of individuals I am quite aware. 
It is a consideration of the fitness of certain (morphological) types of organs 


* Weldon, W. F. R.: ‘‘An attempt to measure the death-rate due*to the selective elimination of 
Carcinus moenas with respect to a particular dimension.”’ Proc. Roy. Soc. Lond. Vol, tv11. pp. 360—379, 
1895. 

+ Weldon, W. F. R.: Presidential address, Zoological Section. Trans. Brit. Assoc. Bristol, 1898, 
pp. 887—902. 

+ Bumpus, H. C.: ‘The elimination of the unfit as illustrated by the introduced sparrow, Passer 
domesticus.”? Biol. Lect. Del. Mar. Biol. Lab. Wood’s Hell, 1897—98, pp. 209—226, 1899. 

§ Beeton, Mary, and K. Pearson: ‘‘ Data for the problem of evolution inman. II. A first study of 
the inheritance of longevity, and the selective death-rate in man.” Proc. Roy. Soc. Vol. txv. pp. 290—305, 
1899. 

Weldon, W. F. R.: ‘‘A first study of natural selection in Clausilia laminata.” Biometrika, 
Vol. 1. pp. 109—124, 1901. a 

§] Di Cesnola, A. P.: ‘Preliminary note on the protective value of colour in Mantis religiosa.” 
Biometrika, Vol. mm. pp. 58, 59, 1903. a 

** Crampton, H. E.: ‘“ Experimental and statistical studies upon Lepidoptera. I, Variation and 
elimination in Philosamia cynthia.” Biometrika, Vol, m1. pp. 113—180, 1904. See also Journ. Exp. 
Zool. Vol. 11. pp. 425—430, 1905. ; 

++ Di Cesnola, A. P.: ‘‘A first study of natural selection in Helix arbustorum.” Biometrika, Vol. v. 
pp. 387—399, 1907. ; : Dh 

tt Kellogg, V. L. and Ruby G. Bell: ‘Studies of variation in insects.” Proc. Wash. Acad. Sci. 
Vol. v1. pp. 324—329, 1905. 
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produced by individuals rather than of the fitness of the individuals with given 
peculiarities. All will admit, I believe, that from the morphological and physio- 
logical view-point this problem is well worthy of careful study. In comparison 
with the results which may be obtained in quantitative investigation of selective 
elimination of individuals, the results should have real interest to students of 
evolution. 


A first study was made of the ovaries of the leguminous plant Cercis*. Here 
the results were negative, but there were certain limitations to the material which 
I have never succeeded in removing. The resu®#s were published for two reasons. 
First, I wanted if possible to interest others in a problem which needs the results 
of several independent researches for its solution. Second, I had the present 
investigation nearly far enough along to see what the outcome would be, and 
I wanted to draw the conclusions from my first set of data quite independently 
and entirely unbiased by the evidence of a second study. 


II, MATERIAL. 


The fruit of the American Bladder Nut, Staphylea trifolia, is familiar to 
botanists as an inflated, membranaceous capsule of (generally) three locules, the 
dorsal sutures of which project to form three lobes, The carpels are (generally) 
united throughout the greater portion of the length of the fruit. From 4 to 





Fic. 1, Diagram showing the external appearance (A) and cross section (B) of fruit of Staphylea. 
Note the one seed and the arrangement of the aborted ovules in the front locule. The cross 


section of a young ovary (C) shows diagrammatically the disposition of the ovules of each locule in 
two rows along the inner angle. 


* Harris, J. Arthur: ‘Is there a Selective Elimination in the Fruiting of the Leguminosae ?”’ 
Am, Nat. Vol. xxii. pp. 556—559, 1909. 
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12 ovules are borne in two series along the inner angle of the locule. None, one 
or a few of the ovules in each cell or carpel develop into smooth, bony seeds. 


The flowers are borne in racemes or panicles. Only a small percentage of 
those formed in an inflorescence mature fruits. I am able to present no reliable 
statistics on this point, but I would say that in the individuals which I have 
studied not over 10 per cent. of the flowers of an inflorescence which ripens fruits 
at all mature their ovaries. What proportion of inflorescences fail entirely to ripen 
fruits I am quite unable to say. 


It has not seemed worth while to secure quantitative data on this point. The 
exact proportion varies from season to season. It is always so small that a very 
stringent elimination is obvious. 


This memoir is limited to two series of data, (1) that especially collected for the 
purpose in 1908, and (2) a small lot of ovaries in different stages of development 
taken for other purposes in the spring of 1906. 


1. The 1908 Collections. 


The collection of this material and the counting of the number of ovules per 
locule was carried out by my two assistants, Miss Eva N. Dixon and Miss Rose 
M. Pechmann, and my sister, Miss Nellie L. Harris. I regret not having been able 
to attend to the sampling personally, but the necessity for an investigation of the 
problem of selective elimination in Staphylea became apparent to me after I had 
left the Missouri Botanical Garden, and the necessity of devoting the spring 
season to experimental plantings at the Station for Experimental Evolution, 
precluded my visiting the Garden at the flowering season for Staphylea to make 
the collections myself. 


As for the countings, I cannot state too highly my indebtedness to these 
assistants for the faithful way in which the arduous task was carried out. In 
all over 21,000 locules had to be opened and the number of ovules counted. For 
the mature fruits the opening was simple, but the counting of the aborted ovules 
is not always an easy task. For the immature ovaries the work was much more 
trying to the patience and eyesight. It was quite out of the question to deal 
with the fresh material, and the difficulty of dissection and counting was increased 
by preservation in alcohol. First the ovaries had to be dissected out of the 
flowers and then each of the locules opened with a needle and the number of 
ovules counted. Possibly some of my fellow botanists would have had the deftness 
of fingers and the patience to carry through the countings, but I would have found 
the task impossible, and even Miss Pechmann and Miss Dixon occasionally nearly 
lost heart during the several months the work was under way. 


The collection of material to determine whether or no there is a selective 
elimination is a difficult problem in a form like Staphylea. The number of 
flowers produced is very: large, and it is quite impossible to take account of 
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even a small proportion of them. The inflorescence is long enough to allow the 
lower ovaries to be developing into young fruits while the upper ones are 
still in the unopened buds. A collection of the flowers before any elimination 
has taken place cannot be made for comparison with the ripened fruits, for 
before the more distal flowers are mature enough for the ovules to be counted 


with safety, the more proximal ones may be falling from the plant or developing 
into young fruits. 


It is desirable that the sample represent not only all parts of the plant, but 
also all regions of the inflorescence, since, should the ovaries be correlated with 
their position on the axis, the constants from material carelessly collected might 
not at all represent the population from which it is drawn. 


The labour of gathering ovaries with a record of their position on the inflor- 
escence seemed to be prohibitive. After careful consideration it became apparent 
that the collection of the spring material must be carried out in a rather rough 
and ready manner, The following was adopted. 


When a shrub was well in flower it was shaken gently to dislodge the flowers 
which had ceased to develop and were ready to fall. Sometimes flowers were on 
the ground which had been dislodged by wind and rain, but these were never 
taken for fear that they had been removed by too great violence. After the tree 
had been shaken and the fallen sample collected, a gathering as nearly random as 
possible of the opened flowers which remained was taken. 

In a few cases only a single collection of fallen and apparently developing 
flowers was made from an individual, but usually the shrub was visited again 
after a few days and a second gathering was made in the same manner as 
the first. 


Finally a series of matured fruits from these individuals was taken in the 
autumn. 

In the two spring series only ovaries dissected from well-opened flowers 
were examined. All buds were discarded for fear that the ovules might not yet 
all be laid down. By discarding any young fruits which fell we avoided one danger 
of overweighting the proximal region of the inflorescence. 

From the foregoing remarks it appears that the 1908 ovaries fall into three 
classes : 

A. Ovaries from opened flowers which were eliminated from the shrub. 

B. Ovaries from opened flowers which were apparently continuing their 
development at the time of sampling. 


C. Ovaries which had completed their development to mature fruits. 


Series B should represent the original population of ovaries more nearly 
than A, since it contains many which would be eliminated later. Because of 
the methods necessarily employed in collecting the samples, I do not look upon 
the comparison of Series A and Series B as very critical. Only a relatively small 
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proportion of the ovaries develop to maturity; Series A doubtless contains many 
which are non-selectively eliminated, while B includes many which a little later 
would have fallen from the tree and been classified as A. 

It seems worth while, however, to carry out the comparison between Series A 
and B. The critical tests are those between B and C and A and C, the most 
significant being the latter. 

In Staphylea the normal type of ovary is 3-merous. Occasionally individuals 
are found which produce a large proportion of 4-merous fruits. Such a one is 
shrub 36, which was omitted in all these comparisons because I wanted to deal 
with as nearly a “normal” population of 3-merous ovaries as possible. I feared 
that the 3-merous ovaries produced on a shrub which bore a very high percentage 
of 4-merous fruits might be in some characteristics—morphological or physio- 
logical—different from those formed by the growing points of an individual with 
almost exclusively trimerous tendencies. 

Shrub 28 bore no flower in 1908. 


In dealing with the raw data it has seemed desirable to treat the material 
for each individual separately, and for the following reason. Rather extensive 
experience in dealing with biometric constants calculated from the fruits of indi- 
viduals has shown that there may be very real differences between those of 
different plants. If each individual did not contribute the same quantity of 
material to the three collections studied, it is quite possible that in dealing with 
lumped samples differences in constants might arise quite independently of any 
influence of selective elimination. I think that with 28 individuals the danger 
of error from this source is not great, but conclusions based on the results from 
28 individual small samples will be stronger than those resting upon merely one 
large composite series. 


2. The 1906 Collections. 


As indicated above these collections were not made primarily for a study of 
selective elimination at all, but for a study of some questions of fertility and 
fecundity. Constants calculated from them indicated the necessity for the detailed 
study attempted in 1908. The results are treated here as supplementary to those 
from the more detailed investigations made later. 

The collection comprises 270 inflorescences* bearing partly developed fruits, 
none of which had reached a length greater than 20 mm., and was made primarily to 
determine the change in the correlation between the number of ovules formed and 
the number developing into seeds at different stages in the development of the 
fruit. The fruit increases in size very rapidly during its early stages of develop- 
ment, and only a rough measure of its size was considered worth while. Three 
classes, 6—10 mm., 11—15 mm. and 16—20 mm. in length, were recognized+. The 

* The correlation between the number of flowers formed and the number of fruits still continuing to 


develop at the time the inflorescences were examined is discussed in Biometrika, Vol. v1. pp. 440, 441, 
1909. 


+ These are exclusive of the short styles. 
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shape of the fruit varies considerably, and the length affords only a rough measure 
of the size, but it is ample for present purposes. 


Fruits which have reached a length of 11—15 mm. may be considered to be 
a selected class as compared with those 6—10 mm. in length—if selection occurs 
at all—and those of 16—20 mm, in length may be regarded as a selected class as 
compared with both of the lower classes. 


III. STATEMENT OF PROBLEMS AND METHODS OF INVESTIGATION. 
In this memoir five characters of the ovary are considered. They are: 
(1) Number of ovules per locule. 
(2) Total number of ovules per ovary. 


(3) Radial asymmetry of the ovary with respect to the number of ovules 
per locule. 


(4) Number of locules per ovary with an “odd” number of ovules. 


(5) Number of locules per ovary. 


To demonstrate that there is a selective elimination for any or all of these 
characters, it is necessary to show that there are differences between the constants 
for the ovaries which develop into fruits and those which do not, of such an order 
of magnitude that they cannot reasonably be referred to the probable errors of 
random sampling. In these pages I follow the rather common example of statis- 
ticians in regarding differences of at least 2°5 times their probable errors as 
significant. In addition I have demanded reasonably constant results from the 
series of individual plants. : 


The working hypothesis of the existence of a selective elimination for any 
character whatever may be tested by determining the difference between the 
constants for: 


(a) A sample from an original population before elimination and a sample 
of the eliminated individuals. 


(b) A sample from an original population before elimination and a sample 
of the individuals remaining after elimination. 


(c) A sample of eliminated individuals and a sample of individuals unaffected 
by the elimination. 


Here “individual,” “sample” and “population” are used in their commonly 
accepted and convenient statistical sense. By individual we understand a par- 
ticular ovary, by sample the collection of ovaries with which the limitations of 
time permitted us to deal, and by population the total ovaries produced by a 


given shrub, or the total produced by all the shrubs investigated, as the case 
may be. 
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A constant calculated from a sample does not necessarily describe precisely 
the population from which it was drawn, for it may be too large or too small by 
an amount depending upon the-chance errors of random sampling. As biologists 
we are interested in our samples only in so far as they indicate for us the charac- 
teristics of the population. The degree of significance to be attached to the 
constant calculated from a sample as a description of a population as a whole is 
determined by the probable error of the constant. The probable error of random 
sampling has nothing whatever to do with the care of the naturalist in doing his 
work, as some biologists seem to believe. Be the observer never so keen and 
conscientious, the significance to be attached to his conclusions is always to some 
extent dependent upon the number of observations before him on his work-table. 
This is merely a conclusion drawn by common sense from common observation. 
Fortunately in many fields of work the mathematician can give us formulae for 
determining the relative weight to be attached to constants calculated from 
different masses of data. These words of explanation are inserted for the benefit 
of fellow-biologists who may still be unacquainted with the true nature of the 
probable error of random sampling so much used by biometricians, for throughout 
these pages I find it necessary to return continually to this fundamental safe- 
guard in the interpretation of the significance, as descriptions of populations, of 
constants drawn from samples of the population. 


In the investigation of the problem of selective elimination on the Cercis 
material, I was necessarily limited to the comparison of the constants for the 
original population with those for the eliminated ovaries, since a severe frost 
prevented any of the ovaries maturing. Here we are able to compare eliminated 
ovaries, the original population after partial elimination, and the matured fruits. 

Further details concerning problems and methods will be taken up under the 
individual problems, to which we now pass. 


IV. Discussion oF PROBLEMS ON BAsiIs OF ACTUAL DATA FOR 1908. 


Problem 1. Is there a difference in the mean number of ovules per locule, or per 
fruit, in the ovaries which are eliminated and those which develop to maturity ? 


Seriations of the number of ovules per locule are given for these collections of 
material in Appendix, Tables XXVIII—XXX. The data for the total ovules 
per fruit are given for the three collections from the 28 individuals in Tables 
XXXI—XXXIII. The means, standard deviations and coefficients of variation 
for the number of ovules per locule, and the standard deviations and coefficients 
of variation of the total ovules per ovary are given for the 28 individuals of the 
three collections and their totals in Tables I—III. 


In comparing the means for the three series it is unnecessary to treat both 
ovules per locule and total ovules per fruit, for the latter is necessarily three times 
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TABLE I. 


Physical Constants for Series A. 











OvuLes per LocuLe OvuLes PER Ovary 
Number siekcbonte 0 TR OE On ee ee oe a eS 
of 
Shrub Average Standard Coefficient Standard Coefficient 
and | Deviation and of Deviation and of 
Probable Error | Probable Error| Variation Probable Error | Variation 
| 

11 8°117 + 022 574+ 016 7°07 + ‘20 1°143 + ‘054 4°70 + ‘22 
12 7:237 + 036 “920 + 025 12°72 + °36 2060 + ‘098 9°49 + *46 
13 7°693+°047 | *856 + °033 11°13 + °46 1°842 + 124 7°98 + °54 
14 6°440+°041 | °752+°029 | 11°69+-46 1°816 + +122 9°40 + °64 
15 7°893+°033 | °842+°023 | 10°67+-°30 1°737 + ‘083 7°34 + °35 q 
16 7°226+°040 | ‘965+:°028 | 13°35 + °39 2°327 +°117 10°74+ °55 
17 7°873+°035 | ‘908+ °025 | 11°53+ ‘32 2°144+ °102 9°08 + *44 
18 7°467 + ‘039 "708+ °028 | 9°49 + °37 1°414+ ‘095 6°31 + °43 
19 7°181 + -038 ‘936+ 027 | 13-°03+°38 2°007 + ‘101 9°31 +°47 
20 =| 6658+°045 | 1:025+-032 | 15°39+-48 2-091 +°111 10°47 + 56 
21 7453+ 043 | “7884-031 | 1058+ -42 1°775 + °120 794434 
22 75704038 | 7324-027 | 9°67 + °36 1°557 + °100 6°86 + "44 
23 =| «6060+ -045 818+ °032 | 13°50+°53 1°894 + 128 10°42 +71 
24 7128 + 045 *888+°032 | 12:464+°45 1:907 +°117 8'924+°55 | 
25 8'154+ ‘061 1256+ °043 | 15°40+-54 3°329 + °198 13°61 + ‘82 
26 69534046 | 8434-033 | 12°13+°18 1°844 +124 8°84+°60 | 
ay | 7067+ -043 ‘894 + ‘030 12°66 + °44 2°128 + +126 1004+ 60 | 
29 | 6663+ -030 781+ ‘021 11°72 + 33 1°763 + “084 8824-42 | 

| 30 | 5-917+4-029 "659 + 020 1114435 | 1:220+-005 6°87+°37 | 

| 31 | 7:237+-030 ‘783 + 022 10°82+°30 | 1°687+-080 7°77 +°37 
32 7713+ 035 723 + 025 9°38+°32 | 1°557 + -092 6°73 + °40 
33 6648+ °035 | °856+ 025 12°87+°38 | 1°946+-098 9°76 + °49 
34 7°053 + 046 *838 + 033 11°88+°47 | 2:082+°140 9°84 + ‘67 

| 35 7490 + 029 741 + ‘020 990427 | 1:513+-072 6°73 + °B2 
37 8°337 + ‘050 12924-03836 | 15°504-44 | 3-2544-155 13°01 + °63 
38 6020+ 025 | ‘6484-018 | 10°76+°30 1°489 +071 824+ °40 
389 =| 7933+ ‘041 668+°029 | 8:42+°37 1°470+ ‘111 6°17 +°47 
40 7118+ °043 *884 + ‘030 12°42 + *43 1°973+°117 9°24 + °55 


the former. The frequencies and constants for total ovules per ovary are con- 
veniently tabled here for a later problem. 


For purposes of comparison we may now take the differences between the three 
series. We use: 


B—A, or the constant for the developing ovules minus the constant for the 
eliminated ovaries. 


C—B, or the constant for the matured fruits minus the constant for those in 
early stages of development. 


C—A, or the constant for the matured fruits minus the constant for the elimi- 
nated ovaries. 


The three comparisons are made in Table IV. 
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TABLE II. 
Physical Constants for Series B. 
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| 7. Cine ewe en 
OvuLEs Per LocuLe OvuLES PER OvARY 
| Number Zo | 
of | | 
| Shrub | Average Standard Coefficient Standard Coefficient 
| | and Deviation and re) Deviation and | of 
| Probable Error | Probable Error Variation Probable Error Variation | 
| | | 
——— | | | 
11 8°157 + 026 672+ 018 8°24+°23 | 1°300+-062 5°31 + °25 
2 | 6°913 + :038 990+ °027 | 14°32+°40 | 2°427+-011 11°70+°57 
| 13 | 7-827+-024 | -6244-017 | 7974-22 | 1-109+4-053 4°72 +23 
L ag 6°573 + °050 912 + 035 | 13°87 + °55 2289 + °154 11°61+4°79 
| 45 8°087 + ‘031 ‘7864°022 | 9°73+°27 | 1:487+-071 6°13 4°29 
16 7141+ 043 1038+ °030 | 14°534+°43 | 2°6084+°131 | 12°17+°62 
17 8137 + 033 847+ °023 | 10°414°29 | 1°795+°086 | 7:354°35 
18 | 7:080+°035 “906 + 025 12°79+°36 | 2°1684+°103 | 10°21+°49 | 
19 | 7°246 + 043 -993+°031 | 13°704-43 | 2-257+-120 10°38 + °56 
20 7°107 + 037 915+°027 | 12°87+°38 | 2°049+-103 9°61+°49 | 
2 7°772+ 025 6044018 | 7°77+°23 1°366 + ‘071 5°86+°30 | 
22 | 7620+ 031 *640+°022 | 8:40+-29 1°252 + -074 5°48 + 32 
23 6°140 + ‘040 721 + 028 | 11°75 + °46 3°544+ -239 19°24 + °34 
24 | 7:026+:046 947+ °082 | 13°48+4°47 2°269+°134 | 10°77+°64 
25 8883 + :046 1°190 + 033 13°40 + °38 25174120 | 9454-45 | 
2% | 7:093+-049 "882 + °034 12°44+ 49 1:929+°130 | 9°07+°62 | 
2y =| «74824 -042 960 + °030 12°91+°40 2°135+°114 958+°52 | 
29 6°383 + ‘034 ‘8744024 | 13°694°38 2-007 + 096 10°48+°50 | 
30 6047 + :032 760 + 023 12°57 + °38 1°6744°087 | 9234-48 | 
31 | 6:860+°032 | -833+-023 1214+ °34 1:960+°093 | 9534-46 
32 7°677+°029 | 7384-020 9°62+°27 | 1°729+-082 | 7-51+°36 | 
33 6°357 + 036 864 + °026 13°59 +°41 1°896+°098 | 9°94+°52 
3 7°089 + 042 1:036 + 030 14°61 + °43 2°3224°117 | 10°92+4°55 
| 36 7°730 + 022 ‘557 + O15 7°21+°20 1:046+°050 | 4°51+°22 
| $7 8°763 + 041 1:046 + :029 11°94+°33 | 2°4424°116 | 9°294°-45 
3 5°833 + 026 668 +018 11°454+°32 | 1°4114+°067 | 8064-39 
39 8°240 + -027 ‘699+ °019 | 8494-24 | 1:429+-068 5°78 +28 
| 40 7°780+ "034 "874+ '024 11°24+ °31 2°006 + 096 8°59 + °41 





The means of the eliminated ovaries will be sensibly identical with those of the 
developing series if there is no selective elimination. If there is a tendency to the 
elimination of those with high numbers of ovules, all three of these differences will 
have the negative sign, while, if it is the ovaries with the lower number of ovules 
which are the more likely to be eliminated, the difference will carry the positive 
sign. These statements apply only to such differences as are greater than those 
liable to arise through the errors of random sampling. When differences are slight, 
the errors due to this source may be considerable enough to change the sign ef a 
constant from positive to negative. The probable errors of the differences are also 
given, and to facilitate comparison the ratios of the differences to their probable 
error, are also tabulated. 

Before considering the ratios of the individual differences to their probable 
errors we may examine the table as a whole, basing our conclusions concerning the 
presence or the absence of a selective elimination merely upon the signs of the 
differences in the mean number of ovules per locule. 
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If there be no selective elimination whatever, the means calculated from the 
three collections should be sensibly identical—that is to say, identical except for 
the differences due to the errors of random sampling. If our series of shrubs were 
very large we should expect half of the differences between the values of any 
two constants to be negative and half to be positive, if these differences were due 
merely to chance differences in the gathering of the fruits composing the samples, 
just as we would expect 50 per cent. heads and 50 per cent. tails in tossing a coin 
a few hundred times. We have only 28 throws of our coin, to use our figure, and 
so we cannot expect to get exactly 14 of each sign, but we ought not to get very 
wide divergences from these numbers unless there is some biological factor at work 
to modify the proportions. 


That there must be some such factor is quite evident from a casual examina- 
tion of the tables. Taking first the differences between the eliminated ovaries and 
those which are continuing their development, we find that the mean number of 


TABLE IIL 
Physical Constants for Series C. 








) See A ee See ee 

















| 
Ovutes per LocuLe | OvuLES PER OvARY 
Number |————_—____—_- i | — — 
| 

a. Average Standard Coefficient Standard | Coefficient 
and | Deviation and of Deviation and of 

| Probable Error | Probable Error | Variation Probable Error Variation 
scree Rema ousted Wace 
11 | 8'360+-030 | 760 + 021 9-09 + °25 1:804 + 086 719+°34 | 

12 | 7:480+°031 *793 + °022 10°61+°28 | 1°657+°079 7°38 + °35 

13 | 6°757 + 036 | 4 15+ °025 13°54+°38 | 1°891+-°090 9°33 + °45 

14 6°588 + 035 *746 + °025 11°32+°38 | 1°664+-096 8°42 + °49 

15 8°713 + 032 *836 + °023 9°59+°27 | 1°980+°094 7°57 + °36 

16 8°190 + ‘033 | *860 + °024 10°51 + °29 | 2°127+°101 8°66 + *42 

17 8°730 + ‘041 1°051 + ‘029 12°03 + °3 2°671+°127 10°20 + ‘49 

18 | 7°587+°029 | °754+ 021 9°94 + 28 L°861 + ‘089 8°18+ °39 
19 | 7°850+°'027 | ‘689+-019 8°77 + °24 1°639 + ‘078 6°964°33 | 
20 7°643 + ‘028 *709 + ‘019 9°28 + *26 1°538 + ‘073 6°71+°32 | 
21 8°067 + °025 654+ °018 8°11 + °22 1°544+ *074 6°38+°31 | 

22 7°903 + ‘020 *511+°014 6°47+°18 1°151 + °055 4°86 + °23 

23 6°261 + ‘039 | ‘777 + 028 12°41 + °45 1°644+ °101 8°75 + °54 

2s | 7°5387+°028 | -713+-020 9°46 + 26 1°378 + 066 6°09 + -29 

25 | 5974-077 1°990 + ‘055 20°74 *59 1°930 + 092 6°70 + 32 

26 7°618+°030 | ‘668+°021 8°77 + 28 1°429 + ‘078 6°25 + °34 

27 8°103+°027 | ‘688+°019 | 8:49+°24 1°347 + ‘064 5°54+°26 

29 | 6:917 +°031 | °802+°022 | 11°59+°32 1°763 + °084 8°49 + °41 

380 6°387+°026 | ‘666+°018 | 10°43+°29 | 1°:294+-062 6°75 + °32 

31 7°463 + ‘029 ‘750 + ‘021 10°05 + *28 1°679 + ‘080 7°50 + °36 

82 7°810+ ‘026 | ‘664+ ‘018 8°50 + *24 1°471+°070 6°28 + °30 

| 8&8 6°833+°030 | °*783+°022 11°45 + 32 1°775 + ‘085 8°66 + *42 

| 34 7°647 + 029 ‘740 + 020 9°68 + ‘27 1°541 + ‘073 6°72 + °32 
35 7830 + °022 ‘578+ ‘016 7°38 + °20 1°127 + ‘054 4°80 + '23 | 
37 8°533 + ‘040 1°021 + ‘028 11°97 + °33 2°112+°101 8:°25+°40 | 
88 | 7:003+°035 “908 + '025 12°96 + 36 2°343 + °112 11°15+°53 | 

| 89 3°290 + ‘029 *743 + ‘020 8°97 + °25 1°540 + ‘073 6°19 + °30 
| 40 | 8°233 + ‘028 *721 + 020 8°75 + 24 1°546 + ‘074 6°26 + °30 
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TABLE IV. 


Difference between Mean Number of Ovules per Locule for 
three Series of Ovaries. 











B—A CO—B C—A 
Number _—— tans l = 
ie Difference Diff. Difference Diff. Difference Diff 
and | ——— and | ——. and ; —— 
Probable Error | / diff. | Probable Error KE diff. | Probable Error | ~ ‘iff. 
-—- —.|—— meaunen taanseeee — A as 
11 *040 + 034 is | *203 + °039 5°14 "2434 °037 | 6°55 | 
12 — *323 + 053 615 | ‘567 + 049 11°47 | °2434°047 | 5:14 | 
13 133 + °053 2°51 | —1°070+°043 24°82 — ‘937 + 059 15°85 
14 133 + 065 2°05 015+ 061 24 1484+°054 | 2°73 
15 "193 + 045 4°31 | 627 + 046 13°71 | °820+°046 | 17°75 | 
16 — 085 + 058 1:46 | 1:°049+°045 23°21 | -964+-052 18°58 | 
17 263 + °048 5°44 593 + °052 11°30 | ‘857+ °054 15°83 
18 333 + °053 6°34 | — 213+ 046 465 | °120+°049 2°46 
19 064 + °058 1°11 604+ °051 11°89 | °669+°047 14°29 
20) -449+°058 | 7-70 5364-047 | 11°50 | -985+-052 | 18-83 | 
21 319 + ‘050 6°35 294+ 035 | 8:50 | -613+-050 | 12°19 
29 051 + 049 1-03 283 + °037 7°68 | 334+ °043 | 7°70 
23 “080 + -060 1°3: "121 + 056 oa "201+ °060 | 3°37 
2, ~ +102 + 064 1°60 511+ 054 9:53 | -409+-053 7°77 
25 *729 + ‘076 9°55 ‘713 + 090 790 | 1°443+-098 | 14°66 
2G *140+ ‘074 1°88 *525 + °057 9°21 | 665+ °055 | 12°03 | 
27 — | ~=S *367 + °060 6°10 “670 + °050 13°51 1°037+°051 | 20°40 | 
29 — *280 + °033 8°51 533 + 046 11°57 "2534+°044 | 5°81 | 
| 80 "130 + 043 3°03 340+ ‘041 8°24 | ‘470+°039 | 12°18 
| 31 —°377+°044 8°46 603 + 044 13°84 2274+ °042 | 5:37 
| 32 — 036 + 045 “80 133 + 039 3°45 | ‘097+°043 | 9-24 
a, — *291 + -039 7°37 ‘476+ °048 10°01 | °185+°046 | 3°98 
34 036 + °063 ‘57 | 558 + °051 10°87 | ‘59347054 | 10°93 
35 240 + 036 6°65 | "100+ ‘031 3°19 | ‘340+ °037 9°29 | 
37 “427 + ‘065 6°59 — *230+ °057 4:04 | °‘197+°065 | 3°03 | 
| 38 — +187 + 036 5°16 1170+ °044 26°65 | ‘983 + 043 22°66 
| 39 “307 + 049 6°22 “050 + -040 1°36 | *357 + 050 7°10 
} 40 “662 + °055 12°13 *453 + ‘044 10°26 | 1°115+°051 21°83 


ovules is higher in the developing series in 20 cases and higher in the eliminated 
series in 8 cases. This is a deviation of only 6 shrubs from the theoretical 14 of 
each class and hence cannot be regarded as very strongly significant. So far as 
the evidence goes, it indicates that ovaries bearing a smaller number of ovules are 
less likely to continue their development towards maturity than those having a 
larger number. 

Even if there were a strong selective elimination of ovaries with a smaller 
number of ovules per locule, we would not expect to detect the full intensity 
of such a selective elimination in a comparison of the two series A and B, for the 
simple reason that we are comparing an eliminated series of ovaries with a series 
in which further elimination is to be expected. 

If it is true that the elimination is not yet complete in the B series, and that 
the failure of ovaries to develop is not random, but selective, we would expect 
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to find differences between the B series and the matured fruits—series C—of the 
same prevailing sign as between the series A and B. 

The results obtained by taking the difference C—B bear out in a most striking 
manner the considerations drawn from the differences B—A, for only 3 out of the 
28 differences bear the negative sign ! 

Finally the most critical comparison is that between the sample of eliminated 
ovaries and the fruits which reach maturity, ie. C—A. The results bear out fully 
the conclusions drawn from the two foregoing comparisons, for only one of the 28 
differences has the negative sign. 

If any biologist wishes to test the number of cases in which twenty-seven 
positive and one negative sign would be expected according to purely chance causes, 
he might take 28 coins and toss them, noting the number of heads and tails he 
gets in each toss in order to determine the number of times that he will find 
27 heads and 1 tail in a few hundred tosses. Mathematical theory says it will be 
about 46 times in 100,000,000. 

The reader may take these chances for what he considers them worth. Person- 
ally I feel that the evidence is rather strongly in favour of the conclusion that there 
is a selective elimination of ovaries and that on the average it is those with the 
larger number of ovules which reach maturity. 

In Diagram 1 I have tried to make patent to the eye the difference between 
the eliminated and the matured ovaries for each of the 28 individuals and for all 
of the material examined. Here the circles show the mean of the eliminated 
ovaries and the solid dots the mean of those which reach maturity. The length of 
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Mean Number of Ovules per Locule for Individuals. 
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Diagram 1. Showing difference in mean of eliminated and matured ovaries for the 28 individuals and 
for the population. Circles=mean of eliminated locules for individuals; solid dots=mean of 
matured locules for individuals. Broken transverse line=mean of all eliminated locules; solid 
transverse line= mean of all matured locules, 
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the line connecting these shows the difference between the two samples for each 
individual. The transverse lines show the position of the means for the whole 
population, the broken one being for the eliminated and the continuous one for 
the matured ovaries. 


The diagram shows one further point—the great difference in the means of the 
individuals. The two transverse lines may be regarded as in a way smoothing 
these irregularities. The solid line is really the mean position of the solid dots on 
the seale while the broken line shows the mean position of the circles, each indi- 
vidual being weighted with the number of ovaries which it furnished. I think the 
differences shown on the diagram are conspicuous enough to carry considerable 
weight. 

Let us now consider the significance of the constants from the individual trees. 
We shall not expect this to be very high as judged by a comparison of the probable 
errors, because of the smallness of the samples of material taken. Unfortunately 
the excessive labour of the dissection and the counting of the ovaries, even in the 
mature state, precluded the examination of more extensive collections. 

The trustworthiness of the difference between any two constants is indicated 
by its probable error. The probable error of the difference between any two 
uncorrelated* constants, say # and y, may be calculated by the formula 

E(«- y)= V Bae 4 Ep. 

Conventionally, a constant is not considered significant unless it is 2.5 or more 
times its probable error. This is the degree of divergence from 0 demanded in the 
following discussion. For convenience of comparison the ratio 

Difference 
Probable Error of Difference 
has been tabled for the three comparisons for each tree (see the third column, 
Table IV.). 

For the comparison B—A it appears that there are 13 significantly positive 
and 5 significantly negative differences. For the comparison C—B there are 22 
significantly positive and 3 significantly negative differences. 

Finally for the most critical test of selective elimination, C—A, there are 25 
significantly positive and 1 significantly negative differences. 

The thing which strikes one particularly about these conclusions is the high 
number of cases in which the differences between the constants are significant with 

* I am afraid that here our constants are not perfectly independent. As Professor Pearson suggested 
to me, this is certainly not true where the size of samples is at all large as compared with the total 


number of ovaries produced, for a random excess in a first sample would be associated with a random 
deficit in the second. 


I do not believe that all of my samples combined equal one-tenth of the total ovaries produced, and 
I suspect that the influence of the correlation of the samples would not be large, but I have no statistical 
proof of this. Indeed in a form like Staphylea it would be difficult to get suitable data for the correction 
of the probable error of the differences in two constants, even if the formulae were available. I think 
these constants affect in no way the validity of the conclusions drawn in this paper, but they may 
render some of the probable errors somewhat questionable. 
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respect to their probable error. In 18, 25 and 26 cases out of the 28 the differ- 
ences exceed 2°5 times their probable errors, and often reach many times their 
probable error. I think relationships of significantly positive to significantly 
negative differences, such as 13 : 5, 22 : 3, 25 : 1 can leave no doubt in the mind of 
the most sceptical that there is a very potent selective elimination, by which the 
mean number of ovules per ovary is very materially increased. 

The combination of the material from all the individuals gives the results for 
the differences in the mean number of ovules per locule presented in Table V. 


TABLE V. 


Comparison of Means for total Material. 





Difference | Ratio of Difference tee a nee 
Comparison and to its = mage pian 
Probable Error Probable Error in Mean | 
| | 
tA. 1426 + 0129 11°054 1971 
C-—B -3693 + 0120 30°775 | 5°005 
C—A 5119+ 0121 42-306 7075 


Here the absolute differences are given in the first colum: , and the percentage 
increase of the higher over the lower mean in the last column. The results 
emphasize the conclusions drawn from the constants for the individual trees. 

The difference between the eliminated sample and the sample of ovaries 
remaining on the tree when they were taken is only about 2 per cent. of the 
former, but the difference is about 11 times its probable error and thus certainly 
significant. 

Doubtless the very slight divergence between the series A and B is due to the 
presence of a large number of ovaries in B that would a few days later fall from 
the tree and be classed with the A’s, for we note a much greater increase from B 
to C than from A to B, amounting to no less than 5 per cent. 

As is to be expected the difference between the eliminated and the matured 
fruits is greatest of all, the latter having about 7 per cent. more ovules per locule 
than the former. 

The reader will note that in the last two comparisons the differences are 31 and 
42 times their probable errors. 

As a further demonstration of the differences in number of ovules per locule or 
per fruit in the three series I give (Diagram 2), the data in the form of an integral 
polygon in which the frequency on each grade is the sum of its own frequency and 
that of all lower grades. This form is one of the most convenient for purposes of 
comparison. 


It is unnecessary to table the actual frequencies for they, can easily be com- 
puted by anyone who requires them from the totals of our tables. The percentage 
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frequencies for number of ovules per locule are given in Table VI, and for total 
number of ovules per fruit in Table VII. This method of presenting the data 
shows very clearly the selective elimination. For instance, in the comparison of 


TABLE VL. 


Percentile Frequencies of Ovules per Locule 
in three Series. 





| Number 








of A B Cc | 
| Ovules 
| 
2 = | ‘Ol =_ 
3 03 | 03 -— 
4 22 “35 04 
5 4°25 4:27 74 
6 27°86 25°19 15°78 
7 | 54°64 | 48°30 31°39 
8 91°50 88-00 59°79 
9 | 98:01 96°59 93°11 
10 | 99°94 | 99°50 99°78 
11 | 100:00 = |_—Ss« 99°98 99°99 
12 — 100°00 100°00 





TABLE VII. 


Percentile Frequencies of total Ovules per Ovary 
in three Series. 











Number 

of | A B C 
Ovules | | 

| 

| = 04 — 

12 | 05 04 04 
is | 05 12 04 
14 19 "24 04 
15 76 77 04 
16 2°10 2°60 18 

17 6°01 6°33 96 

18 16°13 15°58 6°92 | 
19 23°72 22°60 11°69 
20—| 33°17 30°30 18°19 
21 43°44 38°54 24°15 
29 | 57°42 47°99 33°28 | 
23 | 71°88 62°27 45°12 | 
2h OC 88°50 83°04 76°18 

| 2 | 98°75 89°86 84°95 
26 | 96°18 93°87 89°79 
DY 97°57 95°94 92°71 
28 98°90 97°73 94°75 
29 | 99°67 99°03 97°15 
30 | 99-90 99°80 99°63 

3 | 100-00 99°96 99:96 

| g2 | = 99-96 100-00 
33 - 100°00 — 
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number of ovules per locule we note that locules with 7 ovules or less form 54°6 per 
cent. of the collection, while in the collection of matured fruits they form only 
31°4 per cent. In all three series we are dealing with over 2000 fruits, or over 
6000 individual locules, so the results can be looked upon with considerable 


confidence. 


Diagrams 2 and 3 show these results graphically. 





90) 





Percentile Frequencies. 


/ 








en ee ee 
4 6 6 


3 7 4 Q 1 W 12 
Ovules per Locule. 

Diacram 2. Percentile curve showing frequency of different numbers of 
ovules per locule in the eliminated series and the series developing 
to maturity. Dots and broken lines=Eliminated series; dots and 
firm lines= Matured series. 


Problem 2. Is there a difference in the variability in number of ovules per 
locule and per fruit in the eliminated ovaries and those which develop to maturity? 

The results of the preceding section indicate clearly that the chances of an 
ovary developing to maturity are not independent of its number of ovules but that 
those with the smaller numbers are apt to be eliminated. 


It is quite thinkable that a strong selective elimination might exist without 
bringing about any sensible difference in the mean values of a character in the 
population before and after elimination. If the elimination were not of the larger 


or the smaller fruits as measured by the number of ovules per locule or per fruit, 
but of both largest and smallest fruits, i.e. of the extreme variates of both kinds, the 
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mean values of the number of ovules in the matured fruits might not differ sensibly 
from that for the original ovaries, but the variability in number would be reduced. 


On the other hand it is not impossible that a selective elimination might effect 
chiefly the modal region of the curve of distribution, materially flattening it off 
and so increasing the variability as measured by the standard deviation. To the 
biologist this seems theoretically less probable than the elimination of either the 
largest or the smallest ovaries, or the extreme variates of both types, but of course 
all possibilities must be tested out on the basis of actual data. 














Percentile Frequencies. 


Ovules per Ovary. 
Diacram 3. Percentile curve showing frequency of different numbers of 
ovules per ovary in the eliminated series and the series developing 
to maturity. Dots and broken line=Eliminated series; dots and 
firm line= Matured series. 


In comparing the means of the eliminated and uneliminated series, it was only 
necessary to determine the values for either ovules per locule or ovules per fruit, 
since one is necessarily three times the other. For the variabilities, however, it 
seems best to examine the standard deviation for both ovules per locule and total 
ovules per ovary, since one cannot readily be got from the other. 





We are dealing with integral variates, and hence Sheppard’s modification for 
the second moment was not applied, 


Biometrika vir 60 
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First, let us deal with absolute variabilities. 


Selective Elimination in Staphylea 


Standard Deviations. 








The standard deviations for the 


number of ovules per locule in the three collections are given in Tables I—III, 
and are compared in Table VIII. 


TABLE VIII. 


Comparison of Standard Deviations of Number of Ovules per Locule 





EN ee eee seo 


of 
Shrub 


Www WW WWN 
BIH AW HW 


Se Co So So So 
2 Sit Co H 


We 


Difference 
and 
Probable Error 


098 + °024 
069 + °037 
— 232 + 037 
"159 + 046 
— 156+ °032 
‘O73 + 041 
— 061+ °034 
*197 + ‘037 
057 + °041 
— "110+ °041 
— 185+ °036 
— 092 + '035 
— 097 + 042 
“059 + °045 
— ‘065 + *054 
039 + ‘047 
*065 + *042 
“092 + ‘032 
“101 + ‘030 
050 + 031 
015 + 032 
008 + -036 
"198 + 044 
— 184+ 025 
— ‘246 + 046 
020 + 026 
031 + 035 
— ‘010+ ‘039 


per locule, ie. B—A, C—B, 


frequent. 


To render the comparisons more critical, the probable errors of the standard 
deviations and the probable errors of their differences have been tabled. The ratio 
of the difference to its probable error is also given in the same manner as for the 


means, 


in the three Series. 


C—B 

Diff. Difference Diff. 

a and } ——, 
E diff. probable Error | ~ diff 
4°04 088+ °028 | 3°14 
1°87 — "196 + 035 5°61 
6°19 ‘291 + 030 9°54 
3°46 — "166+ °043 | 3°83 
1°75 049+ °032 | 1°56 
1°78 -—‘177+°038 | 4°63 
1°78 "204 + 037 5°50 
5°31 - 152 + 032 4°68 
1:50 | —-304+-036 8°45 
2°67 — °205 + 033 6°22 
5°18 ‘O51 + 025 2-02 
2°62 — *129+ -026 4°97 
2-28 | 056+ 039 1°42 
1°31 | —-234+ 038 6-20 
121 | -800+-064 | 12°51 
. | — -214+4-040 5°31 
1°53 ‘hs 272 + 035 7°74 
2°86 | — ‘072+ °033 2°20 
3°32 | —-094+-029 3°24 
1°57 — 083 + ‘031 2°69 
‘ — ‘074+ 027 2°72 
23 | —-081+-034 2°42 
4°67 — *296 + ‘036 8°13 
7°23 021 + 022 96 
5°37 — ‘025 + ‘040 62 
. “240 + ‘031 7°74 
90 044 + 028 1-56 
25 — 154+ 031 4°95 


Difference 
and 
Probable Error 


+ 026 
+ '033 
+ °042 
7 + °038 
+ 033 
+ ‘037 

"143 + 038 

046 + °035 
— °247 + 033 
— 315 + 037 
— 134+ °036 
— °221+ 031 
— ‘041 + 042 
— ‘175+ 037 

‘734+ ‘070 
—+175+ 039 
— ‘207 + 036 

020 + 031 

‘007 + °028 
— 033 + ‘030 
— 059 + ‘031 
— ‘073 + 033 
— ‘097 + °038 
— ‘163 + 026 
— ‘271+ °045 

‘260 + -031 

075 + (036 
— ‘163+ °036 


Diff. 
E diff. 





eae 


7°09 
3°79 
1 . 





The differences are taken in tle same way as for the average number of ovules 
and C—A. If selection tends to reduce variability 
the negative sign should predominate in the differences; if selection tends to 


increase variability the differences with the positive sign should be the more 
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Studying the tables first for the prevalence of positive and negative signs 
we find the results: 
B—A, 11 negative and 17 positive differences, 
C—B, 18 negative and 10 positive differences, 
C—A, 19 negative and 9 positive differences. 


Taking up next the significance of the differences for the constants for indi- 
viduals as measured by their ratio to their probable error, we find results as follows. 


For the comparison B—A only 12 of the 28 differences are significant as com- 
pared with 18 out of 28 in the case of the means. Of these 6 are positive and 
6 negative in sign. The difference C—B shows 20 significant constants, of which 5 
are positive and 15 negative. C—A shows 17 significant constants, of which 4 are 
positive and 13 negative. The first comparison indicates that selective elimination 
has no influence upon variability. The second and third both suggest that the 
variability is sensibly less after elimination has taken place. 

The differences in the variabilities for the lumped collections may now be 
examined in Table IX. 


TABLE IX. 


| | 
| Difference | 


‘ Diff. 
Comparison and Rafi 
Probable Error a. 
pes . ppl SF es 
3—A + ‘0756 + ‘0091 8°31 
Do all — 0938 + ‘0086 10°91 
Qa — 0182 + :0081 ) 


At first these results puzzled me very greatly; and the following interpretation 
is offered only tentatively and as a working hypothesis for further investigation. 

The collection B probably comes the nearest to representing the original popu- 
lation of ovaries of any that we have, since it still contains a high percentage of 
ovaries which will be eliminated later*. 

Now it is quite clear that if the ovaries with the lower numbers of ovules per 
locule are eliminated, their variability will be expected to be less than that of the 
original population. Likewise if the ovaries with the higher numbers of ovaries 
develop to maturity it is clear that their variability also should be less than that of 
the population as a whole. In short, both the samples of eliminated ovaries and 
the samples of matured fruits represent the more extreme portions of the curve of 
variation. If this be true, we should expect to find the variability of B greater 
than that of A or C, while an inevitable consequence of this condition of affairs 


* I am sorry not to be able to give definite quantitative data bearing upon this point, but this 
is quite out of the question in a ligneous plant like Staphylea, where the number of flowers produced is 
immensely large and where it is impossible to determine with certainty from the mature inflorescence 
the number of flowers originally produced. 


60—2 
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would be that the difference in the variability of C and A would be less than that 
of either of the other two differences. 

This is precisely the condition which we find; for the total material the differ- 
ences B—A and C—B are in opposite senses and clearly significant with regard to 
the probable errors of their differences, while the difference between C and A 
cannot safely be regarded as significant. 

Returning now to the results for the individual shrubs, we note that so far as 
the signs give evidence it is that A and C are less variable than B, just as our 
hypothesis demands. For the differences for individuals which are significant with 
regard to their probable errors, there are 6(A <B) and 6(A>B). For C—B 
there are 5(C > B) and 15(C<B). Finally C—A gives 4(C > A) and 13(C < A). 
These results contain only a partial substantiation for our hypothesis, but no 
evidence directly against it. 

Only further research will show whether these conditions are generally found. 

Take next the results for total ovules per ovary. The comparison is made in 
Table X. 

TABLE X. Comparison of Standard Deviations of total Ovules per Ovary 
in the three Series. 





B—A C—B C—A 
Number — —- 
of : 
Shrub Difference Diff Difference Diff. Difference Diff. 
and — and a and ~— 
Probable Error | “ “'l- | Probable Error | / 4iff- Probable Error | © diff. 
i “156 + ‘082 1°89 504+ °106 4°76 660 + °101 6°51 
| "367 + ‘099 3°71 — ‘770+ -080 9°65 — 040+ °126 "32 
13 — °733+4°135 5°43 "782+ °104 7°49 005 + *153 03 
14 ‘473 + °197 2°40 — 626+ °182 3°43 — -015+°156 10 
| 15 | — *250+°109 2°30 -493+°118 4°18 243 + 125 1°94 
16 ‘280 + °175 1°60 — *480+°165 2°90 — °020+°155 13 
oa — 035+ °133 ‘26 ‘876+ °153 5°72 053 + 163 “32 
18 "754+°141 5°36 | — -308+°136 2°26 045 + °130 34 
| 19 250 + °157 1°59 “617 +143 4°31 — 037+ °127 29 | 
| 20 — 043+ °152 28 — *511+°126 4°05 — 055+ °133 “41 | 
2 — *409+°139 2°94 178 + °102 1°74 — *231+°140 1°65 | 
22 — °305+°124 2°45 — ‘100+°192 1°09 — 041+°114 *B5 
23 1°650 + -271 6-09 1-900 + *260 7°32 — *250+:°163 1°53 
Qh 362+ °178 2-03 - *891+°149 5°97 530+ °134 3°94 
25 — ‘8114-231 3°52 - *587+°151 3°88 | —1°399+-°217 6°44 
| 2 | 085 + +180 ‘47 — "500+ 152 3°30 — °0414+°147 “28 
28 007 + 170 04 — °789+°130 6°05 - °078+°141 "DD 
29 | 2434127 191 | — -2444°127 1:92 — 0004-119 ‘00 
30 | 454+ °108 4°20 — 380+ 106 3°58 “074 + ‘090 83 
a 273 + °123 2-22 — *282+°123 2-29 008 + °114 07 
82 | "1724-124 1°39 257 +°108 2°38 | — -086+°116 ‘74 
33 | — -050+°138 36 | — ‘1214-129 94 | — -171+°129 1°33 
34 241+ °183 1-32 ‘781 +°138 5°67 — 540+°158 3°41 
35 — *467+°188 5°32 081 +°173 1°10 — 386+ 090 4°29 
Sv =| — 8124-194 4°19 — ‘331+ °154 2°15 — 1°142+°185 6°18 
388 | — 078+ "098 “80 "932 + °130 7°15 "854 + °132 6°46 
39 | — 0414-130 31 "112+°100 1-12 071 +°133 “53° 
40 033+ °151 22 — -460+°121 3°81 — *427+°138 3°82 


} 
| 
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The frequencies of positive and negative differences, both for all differences 
and significant differences, for the individuals are shown in Table XI. 


TABLE XI. 


| 
| Aut DIrrERENCES | Sieniricant DIFFERENCES | 
Comparison | —_— cence = —4 
| Positive Negative Positive Negative | 
| | 
B—A 16 12 5 5 | 
C—B 8 20 4 | 13 
C—A 9 19 2 | 


The comparison of the variabilities calculated from the ovaries for all 28 trees 
is made in Table XII. 
TABLE XII. 


Difference Diff. Percentage Increase 
Comparison | and ry | or 
Probable Error E diff. Decrease* | 
| | | | 
A + +2630 + ‘0410 6°41 +9-49 | 
C—B — 2235 + ‘0389 5°75 | — 7°37 
C—A + 0395 + °0387 1°02 +1°43 


From these tables the same conclusion is to be drawn as from that for the com- 
parison for the variability of the number of ovules per locule, namely that both the 
eliminated ovaries, A, and the matured fruits, C, show a less variability than the 
general population from which a portion of the ovaries have been eliminated. 
These results are very clear in the comparison for the totals of the three series ; 
they are not so evident on the small individual samples with their higher probable 
errors. 


Coefficients of Variation. 


In the preceding section it was shown that a selective elimination brings 
about a different mean number of ovules per locule or per fruit in the eliminated 
and the matured ovaries. The ovaries which develop to maturity have a 
higher average number of ovules than those which are eliminated. Possibly, 
then, our conclusions concerning variability may be open to question because of 
the fact that we are comparing the variabilities of samples with generally higher 
means with those with generally lower means. The absolute variability is gene- 
rally roughly proportional to the mean value of a character. 


* Increase or decrease from the eliminated to the uneliminated or from the less mature to the more 
mature. 
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To make our comparisons more comprehensive and to free them in so far as 
possible from the influence of the actual magnitudes, we may now make use of the 
relative variabilities as expressed by the coetficients of variation. 


The differences between the three series for the coefficients of variation for both 
ovules per locule and ovules per fruit appear in Tables XIII and XIV*. 


TABLE XIII. 


Comparison of Coefficients of Variation of Number of Ovules 
per Locule in the three Series. 








Bs | C—B C—A 
Number |——— — ns _ 
. — | Diff. | — Diff. —! Diff. 
ce) | " oO eg ce) sono 
Variation EK diff. Variation E diff, | Variation E diff, 
| 
| 
Jee Gee , ae & 4 
11 1:17 +°30 3°90 | 5 +34 2°50 2°02 + *32 6°31 
12 1604+ '54 | 2-96 + -49 7°57 | —2°114:46 4:58 
13 —3'16+°51 | 6:20 +°44 12°66 2°41 +°50 4°82 
14 2°19+°72 | 3°04 + °67 3°81 | — °36+°60 72 
| 45 — 944-40 2°35 | — +13+°38 "34 | ~1-08 + -40 2-70 
| 16 | 163+4°58 314 | —4:03+°52 7°75 | —2°84+4°49 9°80 
17 | —1:19+°43 2°60 1°62 + °45 3°60 | 90 + °47 1:06 
18 3°314°52 | 6°36 —2°85+°46 6°19 15+ 46 98 
19 | 67 £°57 1°17 ~4:'93+ 49 10-06 | —4°26+°45 9°47 
20 | 252+°61 4°13 —3°59+ °46 7°80 | —6°11+°55 11°11 
| 21 | —2°81+°47 | 5°98 34+ °32 1:06 | —2:46+°46 5°35 
| -1:27+'46 2°76 | ~1:93+°34 5°68 | —3:20+°40 8-00 
23 —1°76+°70 2°51 67+°64 | 1:05 | —1:09+°69 1°58 
| 2 102+ °65 1°57 | —4°02+°54 | 7-44 | —3°00+°52 5°78 
25 — 2°00 + 66 3°03 7°34+°70 | 10°48 5°33 + ‘80 6°66 
26 B31 + ‘69 “45 3°664°56 | 653 | —3°364°56 6-00 
ay 25 +59 42 | —4:424°47 9:40 | —4:17+°50 8°32 
29 1:96+°50 | 3°92 —2710+°50 | 4:20 | — 13+4°46 28 
30 1°43 + *52 2°67 | -2°14+°48 4:50 | — ‘714°45 1°58 
31 1°32 +45 2°93 | —2:09+-44 452 | — ‘77+°41 1°88 
32 24+ °42 57 | -1°124+°36 311 | — *88+°40 2°20 
33 72+°56 1:28 | —2°14+°52 4:11 | —1:42+°50 2°84 
34 2°74+°64 | 422 | -4:934+°51 9°67 | —2:20+-54 4:07 
35 — 2°69 + 32 8°41 18+ 28 64 | —2°52+°32 7°87 
37 —3°56+ °55 6°47 "03 + °47 06 | — 3°53 4°55 6°42 
38 69+ °44 1°57 1°51+°48 3°14 2°20 +°47 4°68 
39 O7 +44 16 "48 + °3B5 27 | 55 + °45 1°22 
40 —1'18+°53 | 2°23 | —2°494+°39 6°38 | — 3°67 +49 7°49 


The results may be tersely summarized for the individuals in Table XV. Here 
the frequencies for positive and negative differences are shown, both for all the 
differences and for probably significant differences, i.e. for those which are 2°5 or 
more times their probable error. 


* Naturally the significance of the sign as indicating the tendency of selection to increase or 
decrease variability is the same for V as for co. 
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For the total material we have, for locules, 








B—A = +°735 + '128, 
C—B =— 961 +116, 
C—A = — ‘225 + °112, 
and for total ovules per ovary 


B—A=+ ‘941 +190, 
C—B =—-1°615 + '175, 
C—A=— 674+ °176. 

These results are consistent throughout, for both ovules per locule and total 
ovules per fruit. They show that, both for the series of individuals and the three 
grand totals, the relative as well as the absolute variability is higher in Series B 
than in Series A or C. The variability of the eliminated series, A, seems to be 











TABLE XIV. 
Comparison of Coefficients of Variation of Number of Ovules 
per Ovary in the three Series. 
; atti a ee ox ens 2 
| | B—A C—B C—A 
foo | See See ae aera. BS Se eS ers ————— a. 
| of | wp Wil ee | 
Shrub | Coefficient Diff. Se Diff, | en Diff. 
0 — o —— oO —_ 
| Variation E diff. Variation | KE diff. Variation | E ditt. 
ee ———E ——E — = a 
ll | 61+ °33 1°85 1'88+°42 | 4:47 2°50 + “40 | 6°25 
12 2°21+ °73 3°03 — 432+ °67 6°44 —2°11+°58 | 3°63 
13 | —3°26+ ‘59 5°52 4°61+°50 9°22 1°35+°70 | 1°93 
1h | 2221+ -02 217 | — 3:19+°93 | 3°43 | — 984-81 1-21 | 
15 | —1'20+ °45 2°67 1°44+ °46 3°13 24+ °50 “48 
16 144+ 83 1°73 | — 3°514°75 4°68 | —2°08+-69 3°02 
7 ~1°72+ -56 307 284460 | 4°73 112466 | 1°70 | 
18 3904 °65 6:00 | — 203463 | 3°22 1:864°58 | 3-21 | 
|} 19 | 107+ °73 1:46 | — 3424-65 | 5°26 | —2°35+°58 4°05 
| 20 | — 864 -74 162 | - 2904°58 | 5:00 | -3°764-64 | 5°87 | 
| 21 | -2-08% -62 3°35 ‘52+ 43 121 | -1:56+°62 | 2°52 
22 | —1°38+ 54 2°55 | — °62+°39 1:59 | —2:00+°50 3°33 
| 23 | 88241-52 16°96 | —10°494°44 7:28 | —1°66+°89 1°86 
| 2h | 1°85+ -84 2°21 | — 4:67+°70 6°67 | —2°83+°62 4°56 
2 | 416+ 94 4°42 | — 2°744°55 4°98 | —6:90+°88 7°84 
6 ‘3+ -86 ‘97 | — 281+°71 3:96 | —2°594+-69 | 3°80 
27 — 46+ -79 58 | — 4:04+°58 6:96 | —4:50+°65 6°92 | 
| 29 166+ °65 2°55 — 1°98+°65 3°04 — °383+°58 ‘57 | 
$0 236+ ‘61 3°87 — 2°47+°58 4°26 — *12+°49 "24 | 
31 | 1°75+ °59 2°96 — 2°03+°58 3°48 — °27+°52 *52 | 
32 | 78+ “54 144 | — 1234-47 | 2°62 | — -454°50 ‘90 | 
33 i 18+ °7i 2°53 | — 1°28+°67 191 | —1:10+-64 | 1°72 
34 | 1°08+ 87 1°25 | — 4:°20+°64 656 | —3°12+°74 4:22 | 
35 | —2:22+ -39 5°69 29 + 32 ‘91 | —1:94+°39 4:97 | 
37 | -—3°72+ °77 4°83 — 1°04+°60 1°73 —4°76+°75 6°35 | 
38 «=| — (18+ °56 32 3°09 + 66 4°53 291+ 66 | 4°41 | 
39 | — ‘39+ °55 71 *41+ °41 1°00 02+ °56 “04 
40 | - 64+ -69 93 | — 2344-51 | 4°39 | -2-984-63 | 4-73 | 
' 











476 Selective Elimination in Staphylea 


somewhat higher than that of the matured series, C, but the difference in the 
relative variability may be in part due to the fact that the mean of A is low as 
compared to that of C. 


TABLE XV. 
| itr Ese ee 2 
DIFFERENCES FOR DIFFERENCES FOR TorTaL 
OvuLes PER LocuLE OvuLEes PER Ovary 
Series 
Compared ype esa ey aa ae 
Positive Negative Positive Negative 


| Differences | Differences | Differences | Differences 
| 





| 
| All Differences | | 
| B—A 18 10 16 12 
C—B 1] oo s 20 | 
C—A 7 21 | 7 21 | 
Significant | 
| Differences | 
| BA 10 8 | 7 
o—38 | 6 16 5 17 
C—A 5 14 2 | 15 
| 


The results for variability are not so conclusive as those for type, but I think 
we may safely conclude : 

(a) That both the eliminated ovaries and the matured fruits have a lower 
absolute and relative variability than the original population of ovaries. 

(b) This is due to the ovaries with tlie lower number of ovules failing to 
develop, and to those with the higher numbers having the better opportunity to 
reach maturity. 

Problem 3. Is there a difference in the radial asymmetry of the eliminated 
ovaries and those which develop to maturity ? 

From the nature of the thin-walled, inflated “bladder” fruit we could not 
expect it to be other than somewhat irregular in form. The externally visible 
irregularities in the evaginations of the three locules are of such a nature that 
no biologist would describe the fruit as asymmetrical. Externally there is no 
conveniently measurable character from which the amount of irregularity might be 
measured. The only character which can be determined with certainty is the 
number of seeds developing or the number of ovules formed per locule. 

Opening the three locules of a fruit quite at random—there being no external 
character to indicate any differentiation whatever between them—we may find 
such numbers of ovules per locule as: 








11—11—11, 
10—11—10, 
8— 9—11, 
9— 7—11, 


1l— 7—6. 
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The first of these is clearly enough a radially symmetrical fruit with respect to 
number of ovules per locule. Any fruit in which the number of ovules is the 
same in each locule, as 6—6—6, or 7—7—7, is radially symmetrical with respect 
to this character. Those in which the number of ovules differs from locule to 
locule are quite as obviously irregular or radially asymmetrical with respect to 
number of ovules. 

Admitting that fruits in which all the locules produce the same number of 
ovules are radially symmetrical, while those with at least one of the locules 
differing from the others in its number of ovules are radially asymmetrical, we 
may be either (a) content to divide our ovaries into two classes, radially sym- 
metrical and radially asymmetrical, or (b) get some measure of the amount of 
asymmetry in individual fruits so that the asymmetrical fruits may be subdivided 
for further analysis. 

The measure of asymmetry must be one for the individual fruit, not for 
a population. The measure must also be independent of the order in which the 
three locules of the fruit are taken, for there is no constant differentiation between 
them and they may be opened in any order. Two measures have occurred to me. 

First, the sum of the positive difference between the number of ovules in the 
locules of an ovary might be used. In a fruit of the formula 


(a) (6) (c) 
Fuk, 


we have, taking all possible differences, 


a—b=-l1, b—c=+2, 
a-—-c=+1, c-—a=-l, 
b—a=+1]1, c—b=-2, 


Sum of positive differences = 4. 

Comparing a fruit of the formula 7—8—7 we find the sum of the positive 
differences = 2. The first is more asymmetrical than the second. This is obvious 
in the present case from mere inspection. 

The second measure is the square root of the mean square deviation of the 
number of ovules per locule from the mean number in the whole fruit. For the 
first illustration the mean number per locule is 7 and the deviations are 

a=0, a? =0, Coefficient of asymmetry = 73 =°8165. 

b=+1, B=l, 

c=-l, c=1. 
For the second illustration: A =7'3333, the deviations are 

a = — ‘3333, b=+ 6666, c= — ‘3333, 
and the coefficient of asymmetry is 
fs + 6666? + 3333" _ 
3 
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This measure is, in short, the standard deviation of the locules of the fruit 
around their own mean. According to this measure the asymmetry of a fruit 
with the formula 7—8—7 bears the ratio to one of the formula 7—8—6 of 
‘471 : ‘817 instead of 2 : 4 as indicated by the interlocular difference method. 

Since the locules of a fruit are undifferentiated it does not matter how we 
arrange them for convenience of treatment. I have always copied the results of 


original countings off in a descending series. The two constants for a few 
illustrative ovaries are: 


Formula Interlocular Coefficient of 

Difference Asymmetry 
l1l1—11—1l = Oo = 0000 
14 9 @ pa an 
n—¢t = etx 9498 
awe 6 = 1-2472 
a= 2 6 = 1°4142 
ag eget Al 8 = 16330 
+ gh ag a 8 = 16997 
a s+ 10 = 1°8856 
— 2 = 10 = 2-0548 
1— 7—6 = o = $-1608 


These illustrations show that while the interlocular difference does not 
distinguish between the amount of irregularity of certain fruits, the coefficient of 
asymmetry, as I have called the standard deviation of the locules of a fruit 
around their own mean, does. By mere inspection I am quite unable to decide 
whether a fruit of the formula 10—6—6 is more irregular, or radially asymmetrical, 
than one of the formula 10—7—5. Both have the same mean number of ovules 
per locule and the interlocular difference is the same for both, but the coefficient of 
asymmetry is slightly higher for the formula 10—7—5. 


I have adopted the coefficient of asymmetry in this paper for the following 
reasons. 


(a) It is merely the standard deviation—so universally employed in modern 
statistical work-—of the number of ovules per locule in an individual fruit. 


(6) It differentiates—whether with quantitative accuracy or not—between 
degrees of asymmetry not distinguished by the interlocular difference method. 


(c) From data tabled in coefficient of asymmetry classes the interlocular 
differences may be obtained by the use of a table. The converse is not true. 


There is one patent objection to the use of the coefficient of asymmetry as 
defined here: it is not independent of number of seeds per locule. There are 
reasons for regarding fruits of the formula 10—6—8 and 5—3—4 as equally 
asymmetrical, but according to our coefficient their asymmetries bear to each 
other the ratio of the square roots of 8/3 and 2/3. 


The simplest method of freeing our constant from the influence of the absolute 
number of ovules would seem to be to take the ratio of the coefficient of 
asymmetry to the mean number of ovules per locule. The reason that I have 
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not done this is that I have used the coefficient of asymmetry in a large series of 
investigations on the relationship between asymmetry and fecundity in this and 
other species*. In them I coukd not well use the relative coefficient because of 
the danger of introducing spurious correlation in determining the interdependence 
of asymmetry and fecundity. It seems desirable, for reasons which will be 
apparent when my data for fecundity are ready to publish, to use the same method 
of describing the irregularity of the ovaries in both studies, 

The chief objection to the use of the coefficient of asymmetry will be removed 
if we show that the magnitude is not greatly dependent upon the number of 
ovules per locule. This must be done for every species investigated. For the 
present series I have discussed this point in Section VI (see page 493) where it 
has been shown that the correlation between the coefficient of asymmetry and the 
number of ovules per locule is very slight. 

Table XVI gives the interlocular difference for all the coefficients of asymmetry 
found in the present study. Anyone who desires may calculate all the constants 
dealt with in terms of interlocular differences by the use of this table and the 
tables of data. 

TABLE XVI. 
Interlocular Differences for Coefficients 
of Asymmetry. 


Coefficient 
of 
Asymmetry 


Interlocular 
Difference 


0000 0 


‘4714 2 

*8165 4 

“9428 4 
1°2472 6 
1°4142 6 
1°6330 8 
1°6997 8 
1°8856 10 
2°0548 10 
2°1602 10 


To determine whether there is a selective elimination depending upon the 
asymmetry of the fruit, the asymmetry of each individual fruit is determined, and 
the distribution of the asymmetries of the eliminated and uneliminated ovaries 
compared, The determination of the asymmetries of individual fruits, a rather 
laborious process at best, is facilitated by tabling the coefficients for all the 
different ovule formulae at the outset. 

A standard deviation calculated on only three “observations”—the three 
locules of an ovary—is of little value as a means of predicting the standard 


* I hope that some portions of these investigations on fecundity will be ready for publication in 
a few months, 
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deviation of another fruit of the same species because of the large probable error, 
but as a description of the irregularity of an individual fruit it is perfectly trust- 
worthy to as many places as we care to carry our arithmetic. I believe the four 
places to which I have tabled are arithmetically correct. For an individual fruit 
or very small sample this refinement is quite without significance, but where 
several hundreds of fruits fall in the same class, as is the case in this work, and 
especially in large series of data on other problems concerning asymmetry to be 
published later, the final decimal places may be of significance. 


TABLE XVII. TABLE XVIII. 
Mean Asymmetry for Individuals. Comparison of Asymmetry. 
Shrub | A B Cc Shrub B—A C—B C—A | 
wae | | m ii 
| 

11 | 305 373 325 11 ‘068 — ‘048 “020 

y | +488 "456 429 12 — 032 - 027 — 059 

3 | 479 346 536 13 = 433 “190 ‘057 
14 325 “400 389 14 ‘O75 — Qll ‘064 
| 15 478 470 424 15 — ‘008 — 045 — 053 
| 16 | *457 ‘470 ‘328 16 012 —"142 | “129 
17 “462 "492 “441 | | 17 ‘030 —051 | —-'021 

18 *397 "435 270 | 18 ‘038 —'165 | -— ‘127 

19 539 538 296 19 | —-001 --242 | —°243 
20 | "644 “494 374 20 —*150 ~°120 | —'270 

| 401 267 | °339 | ?1 — 134 — ‘028 — 162 
a’ <i 419 387 185 22 — 032 — "202 — "235 
23 426 "B59 380 23 - 067 021 — 046 

24 49) *454 107 a4 — ‘036 — ‘047 — 083 

25 487 676 340 25 ‘188 — 335 — ‘147 

26 500 470 359 26 “030 - "111 141 

27 459 514 332 S7 056 — 182 — 126 

29 | +422 “430 466 29 ‘008 “036 043 

30— | $897 366 402 30 — ‘031 037 “005 

31 475 436 400 31 039 — +036 — ‘075 

32 372 *360 287 32 — 012 — ‘073 — ‘085 

33 435 ‘491 417 33 056 — 074 018 

34 599 549 362 34 — ‘051 - *])87 — 238 

35 | 421 326 275 a5 — ‘095 — O51 — ‘146 

37 589 D34 626 37 054 091 — °037 

58 |} *285 “361 *339 &8 ‘O77 — ‘023 “054 

39 «| 36 *403 "425 39 042 021 063 
40 =| *496 463 “380 | 40 — 033 — 083 — ‘116 


The frequencies of the different coefficients of radial asymmetry for each of 
the individual trees are given for the three collections in Tables XXXIV— 
XXXVI. 


The comparison may first be made on the basis of the mean asymmetry of the 
eliminated and the developing ovaries. 


Table XVII gives the mean asymmetry for the three collections from the 
The excessive labour involved in calculating the probable errors for 
these 84 means, where the standard deviations must be calculated for classes so 


28 trees. 
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irregular as our asymmetry grades, has led me to omit the probable errors for the 


individuals, The results for Series A and C are shown graphically in Diagram 4. 








Asymmetry for Individuals. 











Mean Radial 
< 1 ae 
=) 
e 
e 
e 
e 
° 
e 
e 
© 








’ 22 23 24 25 26 27 29 30 2 33 34 35 37 38 39 40 
Serial Number of Individual. 

Diagram 4. Showing mean difference in asymmetry of eliminated and matured ovaries for the 

28 individuals and for the population. Circles=mean asymmetry of eliminated 


ovaries for 
individuals ; solid dots=mean asymmetry for matured ovaries for individuals. 


Broken transverse 
line=mean asymmetry of all eliminated ovaries; solid transverse line=mean asymmetry of all 
matured ovaries. 


The usual comparisons (B—A, C—B, C—A) are given in Table XVIII. 
The differences are taken in such a way that if selection tends to eliminate the 
less symmetrical pods the differences will have the negative sign. For the 
differences we find: 


3—A, 12 positive and 16 negative, 
C—B, 6 positive and 22 negative, 
C—A, 7 positive and 21 negative. 
The mean difference per individual is: 
B—A = — ‘0064, 
C—B = — ‘0675, 


C—A =— 0804. 
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For the grand totals the results are: 


Mean asymmetry, Lot A =°4515 + ‘0051, 
Mean asymmetry, Lot B= 4415 + 0048, 


Mean asymmetry, Lot C = °3724 + 0045, 


B—A =— 0100 + 0070, 
C—B = — 0691 + 0066, 
C—A =— 0791 + 0068. 


The first difference is probably not significant ; the second and third are over 


ten times their probable errors. 


The difference in the distributions of asymmetry for the eliminated and the 
matured series is rather strikingly shown when the frequencies for the two series 


are reduced to percentages and plotted in a polygon. 


In Diagram 5 the 


relatively few asymmetries of 1:2472 and over have been graphed together. 


| | 
“0000 “4714 "6165 “9428 1°2472+- “0000 


Series A. Eliminated, 


| 
| 


"4714 "8165 ‘9428 1°2472+- 


Series C. Matured. 


Diacram 5. Percentage frequencies of chief asymmetry classes in the eliminated 


and matured series. 


Whatever way we arrange our data the conclusions flowing from them are 
the same—that the radial asymmetry of the ovaries which develop to maturity is 
less than that of those which are eliminated. 














| 
| 
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Problem 4, Is there a difference in the proportional number of locules with 
“odd” numbers of ovules in the eliminated ovaries and those which develop to 
maturity ? 


Locules with 4, 6, 8,10 or 12 ovules may be conveniently designated as “even” 
while those with 5, 7, 9 or 11 ovules may be described as “ odd.” 


With respect to this character the fruits fall into four classes : 
3 “even,” 
2 “even” +1 “odd,” 
1 “even” +2 “odd,” 
3 “odd.” 


To the reader it may at first seem to be mere juggling with numbers to seek 
to determine whether there is any selective elimination of fruits depending upon 
the number of “odd” or “even” locules per fruit. Nevertheless purely embryo- 
logical considerations were the incentive to the work. 


I have not been able to find any account of the development of the ovary of 
Staphylea, but from the appearance of the matured organ and from our knowledge 
of the embryology of other forms it would seem quite likely that the ovules of 
a locule belong to the two margins of the same carpel. If this is the case, locules 
with “odd” numbers of ovules must have a different number produced on the two 
margins, while those with an even number almost as surely have the same 
number on the two carpellary margins. If this be the case no one will deny that 
structurally “even” locules are more perfectly formed than the “odd,” being 
bilaterally symmetrical with respect to the number of ovules found on the two 
margins, 


The morphology of the fruit can be thoroughly understood only after an 
embryological investigation, but in the meantime the sensible thing to do seems 
to be to investigate the question of selective elimination with respect to this point 
on the basis of the material in hand. 


A reason for investigating this point lying within our material itself is the 
different frequency of fruits with “even” and “odd” locules. For the totals of 
the three collections the number of fruits belonging to each class is seen in our 
Table XIX. The predominance of locules with an even number of ovules 
suggests that this type may represent a point of morphological stability*. If such 
be the case the determination of whether there is a tendency towards the 
elimination of those with the larger number of “odd” locules is obviously 
important. 


The test may be made in two ways: (a) by the comparison of the percentage 
of “odd” locules in the three series, and (b) by calculating the percentage of each 
of the four different ovule formulae in the three collections. 


* Compare Sir Francis Galton’s suggestion concerning points of organic stability in natural 
inheritance. 
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The percentage of total locules with “odd” numbers of ovules in the three 


series is ° 


Percentage “odd” in Lot A = 37°407 + °412, 
Percentage “odd” in Lot B= 36:065 + ‘377, 
Percentage “odd” in Lot C = 25°185 + °325. 


The probable errors for these frequencies were calculated from the formula 

Ey, = 67449 if Ys X (1 _ Ye) 

m / 

where y,; is the frequency in any one of the four classes, and m is the total 

population. From the probable error for the absolute frequencies so obtained the 

probable errors for the percentage frequencies were got by taking the ratio 
Ey,/m. 


Comparing the series among themselves as usual we get : 


B—A=— 1:342 + °558°/,, 
C—B = — 10'880 + -498°/,, 
C—A = — 12:222 + 524°/,. 


In all these cases there is an actual decrease in the relative number of “ odd” 
locules as elimination progresses. The decrease of 1°3 per cent. noted from A to 
B is perhaps not significant, being only about twice its probable error, but the fall 
of 11 per cent. from B to C is about 22 times its probable error and the 12 per 
cent. decrease from A to C is over 23 times its probable error. 

There can be no reasonable doubt, therefore, that “odd” are more likely to be 
eliminated than “even” locules. Naturally one locule cannot be eliminated 
alone, but it has seemed interesting to work through the results for individual 


locules. We turn now to the more difficult, but more significant, comparison of 


individual formulae. 
The frequencies, absolute and relative, appear in Table XIX. Diagram 6 shows 


TABLE XIX. Frequency of Ovaries of Different Locular Composition 
in the three Series. 





A B C 
Formula | £ 
| 
f f f | ] 
oe | 
3 “even” me 595 28°401 741 30°061 1265 | 46°783 
2 “ even,” l “odd” 7388 37°613 936 37°970 928 | 34°320 
L “even,” 2 “odd” 573 27°351 633 25-680 118 | 15-458 
3 “odd” 139 6°635 155 6288 93 | 3°439 
| ah ae : int ee 
| 
Totals are — 2095 100-000 2465 99°999 2704 | 100°000 
| 
| } | aot 7 | 
| 
Even Locules — 3934 62°593 4728 63°935 6069 74:°815 
Odd Locules vt 2351 37°407 2667 36:065 2043 25°185 | 
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the results graphically. The three comparisons are made on the basis of the 
percentage frequencies in Table XX. 


TABLE XX. 


Formula B—A C—B C—A 


3 “even” ‘ 


+1°660 + 16°722 +18°382 | 


| | 
re aS Mater eee 
| 


2 “even,” 1 “odd” | + °357 | — 3650 | — 3-293 
1 “even,” 2 “odd” —1671 | -—10°222 —11°893 
| 3“odd” | — -347 | - 2-849 | — 3-196 


ae a 


| 
le 
% a 
= 








25 | | i= 
| 
| 
| 
| 
| 
| 
| 














sE 2E+O E+20 30 3E 2E+O E+20 30 sE 2E+O E+20 30 


Series A. Eliminated. Series B. Developing. Series C. Matured. 
Dracram 6. Percentage frequencies of four types of locular composition of ovaries in 


the three series. 


In both of the critical comparisons it appears that elimination brings about an 
increase in the per cent. of ovaries with all three locules “even” and a decrease in 
the per cent. of those which contain one or more locules with an odd number of 


ovules. The differences are also conspicuously large, amounting to 18°4 per cent. 
in the case of the “all even” class. 


In view of the results secured for the individual locules, it has not seemed 
necessary to calculate the probable errors of their differences. A better test would 
seem to be that for the probability against a given system of deviations from 


theory, in the case of a correlated system of variables having arisen through the 
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errors of random sampling*. In making this test we calculate the frequencies to 
be expected in the classes of a series n if they occurred in the same proportionate 
frequency as in the series m, and determine the chances against this being due 
merely to the errors of random sampling from Elderton’s Tables. 

Only the critical comparison C—A need be tested. Table XXI gives the 
data. y* equals 510°8, a value over seven times as large as Elderton has thought it 


TABLE XXI. 











| | | 
F 1 Observed | Calculated , | (m.-—m,’)? | 
| ormula Frequency | Frequency m,—m,’ |: a 
3 “even” ous 1265s 768 -—497 | 321°6 
2 “even,” 1 “odd” 928 | 1017 + 89 | 7°8 
1 “even,” 2 “odd” 418 740 | +322 140°1 
3 “odd” 93 | 179 + 8 | 41°3 
| 
ree ——— 
Totals i 2704 2704 0 510°8 





worth while to table. The chances are thus many millions to one against such a 
divergence as this occurring by the accidents of random sampling. 

Why locules with an odd number of ovules should be inferior to those with an 
even number, or fruits composed in part or altogether of “odd” locules should be 
inferior to those composed chiefly or altogether of “even” locules, our data do not 


permit us to determine. There can be no reasonable doubt of the fact of the 
relationship for our material. 


Problem 5. Is there a difference in the percentage of 2- and 4-merous 
ovaries among the eliminated ovaries and those which develop to maturity ? 

The answer to this question cannot be expected to be very conclusive 
because of the rarity of the anomalies in the series under consideration. The 
available data are given in Table XXXVII. The reason for the frequencies 
for 3-merous fruits being higher in Series C in this Table than in those given 
before is that, to determine more accurately the percentage of anomalous fruits, 
larger samples were examined for this point without counting the numbers 
of ovules. 

For present purposes using only the grand totals, we find the frequencies for 
2-loculed fruits : 

Series A= 2-477 per cent. 


Series B= 2°321 per cent. 
Series C= 1488 per cent. 
B—A=— ‘156 per cent. 
C—B=-— ‘833 per cent. 
C—A=-—- ‘989 per cent. 


* Pearson, K.: Phil. Mag. Vol. u. pp. 157—175, 1900. 
+ Elderton, W, P.: Biometrika, Vol. 1. pp. 155—163, 1901. 
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For 4-merous fruits : 
Series A= 1:376 per cent. 
Series B= 2-243 per cent. 
Series C= 2°461 per cent. 
B—A=+ ‘867 per cent. 
C—B=+ ‘218 per cent. 
C—A=+1°085 per cent. 
No probable errors have been calculated for the reason that our percentages 
are so low that the familiar theory of probable errors cannot be applied. 


It is quite impossible to say, therefore, whether or no these differences are 
really significant, but so far as our data go they indicate that the 2-merous ovaries 
are more likely to be eliminated than the 3-merous, while the 4-merous ones are 
less likely to be eliminated. 


V. Discussion oF PROBLEMS ON THE Basis oF ACTUAL Data For 1906. 


Problem 1. The seriations of the number of ovules per locule are found in 
Table XXXVIII; those for total ovules per fruit in Table XXXIX. 


For constants we get: 

Series A ( 6—10 mm.), Mean = 7:232 + 029. 
Series B (11—15 mm.), Mean=7-660 + °014. 
Series C (16—20 mm.), Mean =7°821 +016. 

B—A = + °428 + 032. 

C—B = +161 +:021. 

C—A = + '589 + 033. 

The increase from A to B is 13°4 times its probable error, that from B to C 
is 8°6 times its probable error and the difference between the mean of the smallest 
and the largest fruits is 17°83 times its probable error. The increase from the 
smallest to the largest fruits is 8°14 per cent. of the value for the former. 

These results certainly substantiate in the most unequivocal way, those secured 
for means from the 1908 series, and are the more interesting from the fact that 
the material was collected not merely for a different purpose and with no idea of 
testing selective elimination, but by a quite different method. 

Problem 2. The absolute and relative variabilities for both ovules per locule 
and total ovules per ovary may again be considered. The complete answer to this 
problem is given in Table XXII. 

As the fruits become larger the variability becomes less, For the absolute 
variabilities the probable errors show that there is a very high probability that 
every difference is significant. Probable errors have not been calculated for the 
coefficients of variation, but the differences are all relatively large and consistent. 
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TABLE XXII. 


Comparison of Absolute and Relative Variabilities 
for 1906. 


Standard Deviation Coefficient 
Material and | of 
Probable Error Variation 
Ovules per Locule | 
Series A one “968 + -020 | 13°382 
Series B bis ‘774+ 010 | 10°104 
| Series C cle 685 + ‘011 8°762 
B—A a5 — 1944-022 | — 3278 
; ces ea — 0894°015 | — 1:342 
C—A mat — ‘283+ °023 — 4°620 
Ovules per Ovary | 
| Series A 2497+ °090 | 11°528 
Series B 1°895 + 042 | 8°247 
Series C se 1594+:046 6°795 
3— A me — ‘602+ °099 | — 3281 
C—B aoe — ‘301+ ‘062 — 1°272 
C—A he — -903+4°101 — 4°553 


Problem 3. The frequency of the different asymmetry classes in the three 
series is given in Table XL. The three differences of the relative frequencies are 
given in Table XXIII. 

TABLE XXIII. 





B—A | C—B | C—A 
| 
0000 +3°751 +8°495 12-246 
‘4714 +3064 —9°790 — 6°726 
8165 — 4°87 +2°067 — 2°810 
9428 - 1-790 — 1-077 — 2867 
1°2472 — — - 





e ’ = \ 


The results show that as the size of the fruit increases the symmetrical fruits 
form a higher percentage of the entire population, while the asymmetrical fruits 
of the various grades become relatively fewer. 
irregularity is to be expected. 


Our numbers are small and some 


The average asymmetries are : 
A, 6—10 mm. long, 
B, 11—15 mm. long, 
C, 16—20 mm. long, 


Mean = ‘3822 + ‘01 
Mean = ‘3381 + ‘0€ 191. 
Mean = ‘3025 + ‘012 


B—A = — 044 + 018. 
C—B = — ‘036 +:016. 
C—A = — ‘080 + 020. 
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There is clearly a decrease in asymmetry from A to B and again from B 
to C. Owing to the smallness of our series the probable errors are high, but 
all the differences are over twice their probable errors, and the critical difference 
C—A is four times its probable error. 

Diagram 7 shows graphically the frequency of the four chief asymmetry classes. 
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‘0000 *4714.—Ss “8165 = *0428-++- > ou wes comet “0004714 “8165 “0428-4 
Series A. Youngest. Series B. Intermediate. Series C. Oldest. 
Dracram 7. Percentage frequencies of chief asymmetry classes in the three series of 
developing fruits for 1906. 

Problem 4. The frequency of the four types of ovaries with respect to the 
number of “odd” locules produced is shown for the three lengths of fruit in 
Table XXIV. and graphically (with percentage frequencies) in Diagram 8. 


TABLE XXIV. Frequency of Ovaries of Different Locular Composition 
in the three Series. 


6—10 mm, =A 11—15 mm, =B 16—20 mm. =C | 
Formula a | 
| | f | %, f 1 f ‘ | 
| 3 “even” .. | 58 | 33-333 169 35°957 125 45°787 | 
2 “even,” 1 “odd” 68 | 39:080 162 34468 85 31136 | 

1 “even,” 2 “odd” 35 | 20-115 107 22-766 51 18°681 


3 “odd” 13 | 7°471 32 | 6-809 12 | 3:396 | 
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The reader will note, by comparing the percentage frequencies for the smallest 
and the largest fruits, that the 3 “even” class is much more abundantly repre- 
sented in the largest than in the smallest fruits. Those with one or more “odd” 
locules are necessarily less abundant. 


The percentage frequency of “odd” locules in the three series is: 
Lot A, 6—10 mm. = 33°91 +1°40°/,. 
Lot B,11—15 mm. = 33°48 + ‘85° 


°° 


| 


Lot C, 16—20 mm. = 27°23 4+1°05°/.. 
The difference, C—A = — 6°68 + 1°75°/,. 
In predicting the frequencies of the four classes in Series C from the distri- 


bution found in Series A, to ascertain whether the former could possibly be 


regarded as arising from the latter through the errors of random sampling, 
we find 


x= 207, P =00012. 








ee ae 


3E 2E,10 1E, 20 30 3E 2E,10 1E, 20 30 sE 


L 








2E,10 1£,20 380 
Series A. Youngest. Series B. Intermediate. Series C. Oldest. 


Dracram 8. Percentage frequencies of four types of locular composition of ovaries in the 
three series of developing fruits for 1906. 


Thus deviations of observation from theory so great as those found here would 


be expected to arise through the errors of random sampling only about 12 times 
out of 100,000. 
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Allowing full weight to the fact that the frequencies in series A have their 
probable errors, I think we can maintain that series C and series A are differ- 


entiated from each other with réspect to the number of “ odd” and “even” locules 
in their ovaries. 


Problem 5.° This problem cannot be discussed, since no data were collected 
for other than 3-merous fruits. 


VI. INTERRELATIONSHIP OF CHARACTERS CONSIDERED. 


For convenience and clearness our characters have been treated in the fore- 
going discussion of problems as though they were quite independent. 


Analysis can be carried somewhat further. Ifa character is quite uncorrelated 
with others, we shall be confident in regarding any selective elimination which we 
find associated with this character as arising through some unfitness for continued 
development associated with it. If, on the contrary, it appears that a character is 
correlated with some other, we cannot know without further evidence which of a 
pair is potent in producing a selective elimination. 


These points will best be made clear by an examination of the three relation- 
ships which seem worth consideration for our material. 


1, The Correlation between the Coefficient of Asymmetry and the Number of 
“odd” Locules per Ovary. 


The fact that both radial asymmetry as measured by the deviation of the 
number of ovules per locule from their mean and the composition of the ovary 
with regard to the “even” and “odd” nature of the individual locules seem to be 
of significance in determining whether it shall develop into a mature fruit, does 
not necessarily prove that with both of these characteristics there is associated 
some functional unfitness for development. 


There is necessarily some correlation between the two characteristics in any 
fruit. If the three locules have all the same number of ovules and they be even, 
say 8—8—8, the asymmetry will be 0000. If it be “2 even, 1 odd,” say 8—8—7 or 
9—8—8, the asymmetry will be -4714. Thus a fruit composed of “odd” and “even” 
must necessarily be somewhat asymmetrical while one of “all even” or “all odd” 
may be perfectly asymmetrical. If an “all even” or “all odd” fruit is radially 
asymmetrical with respect to number of ovules per locule, it may be more asym- 
metrical than a fruit of locules of both types, since there must be a difference of 
two ovules between two different “even” locules, while an “even” and an “odd” 
may differ by only a single ovule, For instance, 8—8—6 is more asymmetrical 
than 8—8—7. 
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As the reader will see, the problems are very complex. Even the relatively 
large series of data available for this research are, I fear, insufficient for a detailed 
study of the questions of the interdependence of these two characteristics of the 
fruit and their fitness. Later I hope to publish extensive data already on hand 
that may throw some light on these questions. For the present I shall only show 
the correlation between the coefficient of asymmetry and the number of “odd” 
locules per fruit. The data for the three series appear in Tables XLI—XLIII, 
and the constants by the product moment method* are given in Table XXV. 


TABLE XXV. 





Correlation and 
| Probable Error 


— —__—__—_—_ —] 


Series A... *2247 + 0140 
Series B ee +3296 + 0121 
Series C ood 3282 +: °0116 


Clearly there is a definitely significant, though not very large, interdependence 
between the two characteristics. 


It may be possible, therefore, that one of these characteristics has compara- 
tively little independent significance in determining the fitness of an ovary for 
continuing its development. If only one of the two characteristics is thus potent 
while the influence of the second is only apparent because of its dependence upon 


the first, we cannot determine at present which is independent and which is 
dependent. 


Personally, I believe that the coefficients of correlation between them are too 
low to account entirely for the results for elimination that have been secured, and 
that therefore they are both concerned in determining whether or no any ovary 
shall develop. My reasons are as follows: 


If the elimination of the ovaries of one characteristic, a, of the two here 
considered be not primarily due to their possessing their first characteristic but 
merely apparently connected with it because it is itself correlated with the second 
characteristic, b+, we should expect no selective elimination with respect to the 
dependent character within the subgroups of the independent character. Con- 
cretely, if the elimination of ovaries with a larger number of “ odd” locules is due 
solely to the fact that these ovaries are also more radially asymmetrical, we would 
expect to find no elimination with regard to number of “odd” locules when we 
work within the same asymmetry class. Similarly for elimination with respect to 
radial asymmetry. Tersely: for constant a is there an elimination depending on 6; 
for constant 6 is there an elimination depending on a? 


* Calculated without any combination or modification of classes for the coefficients of asymmetry. 
+ Pearson, K.: Phil. Trans. A., Vol. cc. pp. 18, 19, 1902, has termed this indirect selection. 
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Data are available in Tables XLI—XLIII. Only the critical comparison 
C—A need be made. Table XXVI. gives the difference for the mean number 





TABLE XXVI. 

Radial Mean “Odd” 

| Asymmetry Locules, C—A 

“0000 — °335 
"4714 — 106 
“8165 — 115 
“9428 — 454 

1°2472 — °157 


of “odd” locules for the first five radial asymmetry classes. Beyond this the 
frequencies are too few to be trustworthy. Table XXVII. gives the difference 


TABLE XXVII. 
' 


Character of Ovary 


Mehn Radial 
Asymmetry, C—A 


3 “even” eee —°017 

| 2 “even,” 1 “odd” | — ‘019 
1 “even,” 2 “odd” — 034 

3 “oda” — ‘077 


between the mean radial asymmetry of C and A for the four classes of ovaries 
with respect to number of “odd” locules per ovary. 


I attach no importance to the numerical value of these differences, for I have 
not calculated their probable errors, but the negative sign throughout seems to 
me fairly satisfactory evidence that both of these characters are to some extent 
of independent significance in determining whether an ovary shall develop to 
maturity. 


2. The Correlation between the Coefficient of Asymmetry and the Number 
of Ovules per Locule. 


Both the actual number of ovules per locule—or per ovary—and the radial 
asymmetry of the fruit with respect to the number of ovules per locule seem to 
have significance in determining the fate of an ovary. Are they really both 
significant, or is one dependent upon the other ? 








This is a fundamental physiological question. I believe a trustworthy answer 
is given by the coefficient of correlation between the asymmetry of the fruit and 
the number of ovules per locule. Series B is the most like the original popula- 


Biometrika vi 63 
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tion of ovaries and serves for the determination. The coefficient calculated from 
the surface given as Table XLIV is 


r ='0221 +0078. 


Possibly the coefficient is statistically significant, being about thrice its probable 
error, but it is clearly of so low an order that it can have no appreciable biological 
significance in the problem now in hand. 


Apparently, therefore, both number of ovules and radial asymmetry in the 
distribution of the ovules in the three locules are independent and physio- 
logically significant factors in determining whether an ovary shall or shall not 
reach maturity. 


€ 


3. The Correlation between the Number of “odd” Locules per Ovary 
and the Number of Ovules per Locule. 


Again, series B may serve as material. The correlation surface is shown as 
, 
Table XLV. The constant is 


r = —'0368 + ‘0078. 


With a correlation so low as this it seems reasonable to think that, in their rela- 
tion to selective elimination, these two characters are practically independent. 


Vil. RECAPITULATION. 


1. The foregoing pages contain an account of a quantitative study of the selec- 
tive elimination occurring in the ovaries of Staphylea trifolia between the time of 
flowering and the maturing of the fruit. The scope of the paper is limited strictly 
to the presentation of observed data and their description by means of biometric 
constants. No biological theories to explain why one ovary should be less fit than 
another are suggested, and no hypotheses concerning the influence of the selective 
elimination demonstrated upon evolutionary progress are brought forward; there 
will be ample time for this after the much more extensive data already collected are 
published. A comparison of these results for the selective elimination of organs 
with those secured by others for a selective elimination of individuals may also be 
profitably postponed. 

2. The ovaries with relatively low numbers of ovules are more extensively 
eliminated than those with high numbers. The mean of the population re- 
maining after selection is about 7 or 8 per cent. higher than that of the 
eliminated individuals. 


3. The variability of both the eliminated ovaries and the series remaining 
after elimination is less than that of the original population. 


This condition 
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would be expected as the result of the elimination of the ovaries with the 
generally lower number of ovules, and the continued development of those with 
the generally higher numbers, a8 described immediately above. 


4, The ovaries which remain after elimination are more radially symmetrical 
(as measured by the coefficient of asymmetry, or the standard deviation of the 


number of ovules per locule around their own mean) than those which are 
eliminated. 


5. Ovaries with one or more locules with an “odd” number of ovules are 


more likely to be eliminated than those with all the locules bearing an “even” 
number. 


6. Dimerous ovaries seem less likely, and tetramerous ovaries more likely to 


develop to maturity than the normal trimerous ones, Further data bearing on 
this point are needed. 


7. Radial asymmetry and the composition of the fruit with respect to number 
of “odd” and “even” locules are not independent, but correlated characters. Evi- 
dence is adduced to show that both are to some extent of independent significance 
in determining whether or no an ovary shall continue its development. 


8. Neither radial asymmetry nor the composition of the ovary with respect 
to number of “odd” and “even” numbers of ovules per locule is closely enough 
correlated with number of ovules per locule to modify in any way our conclusions 
concerning them as independent characteristics in their relation to selective 
elimination. 


BIOMETRIC LABORATORY, 
University CoLiecr, LoNDON. 
January 31, 1910. 


y 
VIII. AppeNDIX oF TABLES OF DarTa. 


In publishing these data and the constants calculated from them I have 
retained for the individuals the original numbers of my notes, 11—40, since I shall 
later treat other features of the same individuals and it may be a convenience to 
the reader who cares to make comparisons, as well as to myself. 

In titling the Tables the following symbols were used : 

L=Number of ovules per locale. O =Total number of ovules per ovary. 
N=Number of locules per ovary. LC=Locular composition, ie. the com- 
position of the fruit with respect to number of locules with “odd” and “even” 
numbers of ovules. A=Asymmetry, as measured by the coefficient of asymmetry. 


A, B and GC, following the year of collection, refer to the categories explained 
on p. 456. 


63—2 
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Selective Elimination in Staphylea 


TABLE XXXVII—N 1908, A—C. 











Locules per Ovary in three Series. 
| 
A B Cc 
Tree | | | | bs | | | S 
2/3 |4)5|2] 814) 5) 2 | ae ae. 
| 
| ww baal 
11 | 32| 100} 1] —| 14] 100 -~| 5] 104| — | — 
12 | 100 | — | — | — | 100) —| —| —| 109) 1) — 
13 |—| 50} 1|—|— | 100 3 | 2) 118 - 
14 9 50}; —|—| 10 | 50 14 75 — | 
| 15 |—|100| —|—J] 1} 100 — | — 1.306 | — | 
. oe fo 90 | oe | 5| 90 | — at fhe 
ee 3 — | 100| — Bow oe | 100 -|-— - | 109 1 | 
18 = i | 2 | 100 6] — 1} iki 3 | - 
19 5 90 | — 6 80 =f 11 153 
20 _— 80 — gO | — —|lll | — 
21 . Se | i) oe) 2) =P — ee St = 
2 |—| 55| 3| —| 65} 1) —| 102 | 1 
D3 es 50 2 | 50 1 | |— 61 l 
| 2h — 60 5 | 2 65 - i}—j}110); —| — 
| 25 - 65 - | - 100 1 we 163 4 
| 26 — 50 50 | — -|— 76) — ~ 
rY 65 | — | 80 |} —| — | 122} —| —| 
29 | — | 100 100 A - 123 | 2\|—| 
30* eo | “- 85 1} — 139 | 2 | 
| $1 —|100;} 2/|- 100 l 105 3 | | 
| 32 -| 6 | 3| 1 100 | 9 184 | 26 | — | 
; 5| 90|—| o) 6 | = 13 | 134 | — 
3, - 50 3 | -}| 90) 16) 2 |- 118 | 16 
35 - 100 | 10 100 | 13 | 126 18 
if 100 | — | - 100 | — - | 124 —|}— 
38 2 | 100 - 18 | 100 3 | 153 - 
29 10 100} 3)—| 1]135| 3] —|] 
40 65 100 l 158 | | 
| 
54 2095; 30 l 60 | 2465) 58 ) 52 | 3355) 86 


* Shrub 30, Series C has one fruit intermediate between 3- and 4-loculed condition. 














TABLE XXXVIII—L 1906. 
Ovules per Locule. 
5 6 ri 8 9 | 10 11 | Totals 
6—10 12 | 127 142 214 93 | 4 a 522 | 
11—15 2 104 393 810 76 24 1 1410 
16—20 2 28 168 553 53 | 15 819 | 
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TABLE XXXIX—O 1906. 


Total Ovules per Ovary. 








| 17 | 18 | 19 | 20 | 21 | 22| 23 | a | 25| 26| 27 | 28 | | 














h y | 26 y | 2 | 29 | 30 | Totals 
A 7 /15]19|18| 21 | 21 | 21 | 35 4) 418 \—|-) = oe 
B —| 6|11| 32] 40] 77| 107/146 | 92/10] 7| 4} 7]1]4 470 
C —| 1] 3) 3}18| 41] 54/19/15] 6) 6] 5) 1]1] avs | 
| | ! ee, 
TABLE XL—A 1906. 
- | 
A B c 


| a oh ae 4 2 , ; 


-0000 59 33°908 | 177 37°659 126 46°154 


‘4T14 85 48°850 244 51°914 | 115 42°124 
“8165 17 9°770 23 4°893 19 6°960 
“9428 12 6896 | 24 5°106 11 1-029 
1:2472 l “574 2 "425 2 733 
es! = | ee Bee, eee. = 
Totals | 174 — | 470 lb 273 5 | 


TABLE XLI—ALC 1908, A. 


Coefficient of Asymmetry. 








: | "0000 | 4714 | *8165 | 9428 | 1°2472 | 1°-4142 | 1°6330 | 1°6997 | 1°8856 | Totals 
Ss | 3 even ie 4 462 os a ne = ons 2 a 1 595 
= & | 2 even, 1 odd - 614 138 - 21 14 1 ~ 788 
5 &.| 1 even, 2 odd — 443 95 | 3 8 — 5 — 573 
me £1] 30dd ies 103 — — 35 a _ ] -- _ 139 

eal 
2 | Totals wt 565 | 1057 | 233 165 43 | 22 3 6 1 2095 
{ ? Cee, "sees" ke ean ' 








TABLE XLII—ALC 1908, B. 


Coefficient of Asymmetry. 

















-0000 ‘4114 | 8165 9428 | 1-2472 | 1-4142 | 1-6330 | 1-6997 | 1-8856 | 2-0548 | 2-1602 | Totals| 
- ! | | \ \ | 
iS | 3 even -- | 610 | — | — {|13| — — a 17 1 | — 741 | 
“= £ | 2 even, 1 odd -- 739 | 146 — | 85 9 — 6 — ee 936 | 
2 6} leven,2odd | — | 451 | 1422 | — 28 7 — 2 _ Tet 633 | 
WH) 3odd ..f io; — | — | 42} — | — 3 = eee ee Ce iat 

° | 
© se | | | © 9 5 
Totals ... | 720 | 1190 | 288 | 172 | 63 16 | 3 tei 1 | 2465 

\ | 1 = ino ae 
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Selective Elimination in Staphylea 


TABLE XLIII—ALC 1908, C. 


Coefficient of Asymmetry. 









































. | “0000 | “4714 | *8165 | 9428 | 1°2472 | 1°4142 1°6330 | Totals 
"S| B3even .. 4) 1001; — - | 263 - 1 1265 
ra ‘a | 2 even, 1 odd —- 744 160 20 4 — 928 
2 a.| 1 even, 2 odd - 340 66 . ll l a 418 
4 & | 3 odd 76) — - 17 - — -- 93 
7 | Totals 1077 | 1084 | 226 | 280 | 31 a 2704 
TABLE XLIV—AL 1908, B. 
Ovules per Locule. 
2 3 h 5 6 8 9 10 | 11 | 12 | Totals 
Ss | 
=| 000 }— | — 18 | 555 | 255 | 1242 | 57] 33 -|—f 2160 | 
s | ‘4714 —| 5 | 169 | 652 | 1089 | 1191 | 381 | 81| 2] —J] 3570 
= | °8165 7| 48/161] 248| 232] 1296| 33] 8| 1 864 
>| °9428 at ah St) ee 66 | 206| 30] 40 | — 516 
<q | 12472 -) 2) ). Bi 36 44; 28) 18] 8] 1 189 
a, | 1°4142 l 6 3 6 16 6 5| 5 | 48 
2 | 1°6330 - l - 3 _ | = 2 | 9 
S | 1:6997 | 1 l - | 1 3 1 4 ae ee 24 
2 | 1°8856 Z 2 - | —] 3 
os | 270548 | — | — 2 l 2 l 2| 3 9 
$ 2°1602 ~ - l l - — — lj - 3 
oO 
Totals | 1 1 | 24 | 290 | 1547] 1709 | 2936 | 635 | 215 | 32 | 5 7395 
TABLE XLV—LCL 1908, B. 
Ovules per Locule. 
2 B Y 5 6 7 8 9 10 11 12 {Totals 
‘3 | 3 “even” — 7 694 1448 | 74 | — | — | 2223 
‘a | 2 “even,” 1 “odd” | 1 l 15 | 104 | 610 | 607 | 1144] 214 97 10) 5 | 2808 
S,| 1 “even,” 2 “odd” 2 | 141 | 243 | 778 | 344] 331] 44 / 16 | — | 1899 
gS | 3 “odd” 45 - 24 90 | — 6|—] 465 
° 
i) : aie ss | i 
Totals 1 | 1 | 24 | 290 | 1547 1709| 2936) 635 | 215 | 32 | 5 | 7305 
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A BIOMETRIC STUDY OF PHAGOCYTOSIS WITH 
SPECIAL REFERENCE TO THE “OPSONIC INDEX.” 


SECOND MEMOIR. ON THE DISTRIBUTION OF THE 
MEANS OF SAMPLES. 


By M. GREENWOOD, Junr., anp J. D. C. WHITE. 


(From the Statistical Laboratories of the Lister Institute of 
Preventive Medicine and the London Hospital.) 


In a recent memoir* we communicated the results of an analysis undertaken 
with the object of throwing light on difficulties associated with the ordinary 


method of estimating opsonic indices. The main results of that investigation can 
be summarised in the following way: 


(1) Phagocytic counts have pronouncedly skew frequency distributions and 
are good examples of Pearson’s skew curves. 


(2) The means of small samples have likewise a skew distribution. 


(3) The use of a thicker bacillary emulsion while diminishing the skewness, 
does not, under possible experimental conditions, do so sufficiently to allow of 
testing on the basis of a Gaussian curve. 


With respect to (2), we pointed out that results based, as were ours, upon 
samples extending in no case to more than 2000 cells, were provisional and that 


we intended to re-investigate this part of the subject with the help of more 
adequate material. 


The present memoir comprises a study of the distribution of a large sample 


of phagocytic cells, the distributions of the means of sub-samples and some 
preliminary contributions to the problem of mean distributions in general. 


The basis of this work is a count of 20,000 cells made by one of us (J. D.C. W.) 


who had had 18 months’ previous experience in the Inoculation Department of 
the London Hospital. 


* Biometrika, Vol. v1. p. 377. 


i= 


my 


t 


SSS 








506 On the Distribution of the Means of Samples 


In undertaking this large count, those methods were adopted which are usually 
advised in the preparation of films. The mixture, to be subsequently placed on 
the slides, was made up of three equal parts, (1) serum, (2) corpuscles, (3) bacterial 
emulsion. The first two were derived from the operator and, in the case of the 
corpuscles, care was taken that they should be equally distributed, the upper 
portion of the centrifugalised blood being taken off and thoroughly mixed in a 
small tube. The bacterial emulsion was prepared as usual and carefully mixed so 
as to secure the greatest possible homogeneity. Since we desired a large count of 
a homogeneous character, not a comparison of different mixtures, it was decided 
to use one pipette only, but of larger calibre than usual, in order to secure the 
greater amount of mixture requisite for the preparation of a large number of 
slides. This pipette was throttled for ease of manipulation. The proportions of 
the three materials were carefully mixed and the pipette placed in an incubator 
for fifteen minutes, the emulsion having been designed to give an average of three 
to four bacilli per cell. The pipette was slightly rotated from time to time in 
order to keep the corpuscles from settling. Sufficient slides being in readiness, 
the pipette was withdrawn from the incubator, the contents again thoroughly 
mixed, and films were prepared as quickly as possible. These were then stained 
with Aniline Fuchsin and Methylene Blue; thereafter cover-slips were applied 
with Canada Balsam to preserve the slides during the process of counting. The 
general quality both of films and staining was exceedingly good, very few slides 
being of an inferior character. 


The counting of so large a number of cells as 20,000 was, needless to say, 
sufficiently laborious, especially in regard to the leucocytes containing more than 
10 bacilli apiece. Reasonable accuracy was found possible up to 15 and the few 
cells containing more than this were marked 16. A mechanical stage was used— 
the only method, in the operator's opinion, which can secure that the same cell 
shall not be counted twice over. The only cells omitted were those the outline of 
which was indistinct, or where for any reason, e.g. the clumping of cells or bacilli, 
it could not be determined how many bacilli lay in an individual cell; the possible 
error introduced by this latter criterion will be discussed below. 


Definite fragments of bacteria were counted as bacteria; so also were bacilli in 
definite contact with the periphery of cells. In order to avoid mental prejudice, 
addition of the rows of figures was postponed generally for days and always at 
least until the end of the day’s work. In this way 20,000 cells were counted in 
the course of about a month, with occasional intervals. The sheets of figures were 
then taken and the totals of successive twenty-fives, fifties and hundreds were 
determined ; then the frequency of each number per cell on each sheet of 500 cells 
was tabulated, this latter process incidentally eliminating any arithmetical mistakes. 


The frequencies on each of the forty sheets were then added together and the total 
frequency obtained. 


The actual distribution of the twenty thousand cells is set forth in Table I. 
Proceeding to fit a curve, without using Sheppard’s corrections, the constants of 
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Table II were obtained and Graph 1 shows the theoretical curve and observations, 
It is evident that, although the curve represents the general trend of the observa- 
tions with fair accuracy, the goodness of fit as evaluated in the ordinary way is 
poor. This matter appears to require comment and explanation. 


TABLE I. 


Actual Frequencies of 20,000 Cells, together with Corrected 
Values determined on the Basis of a Random Count of 
2,000 Cells from the same Material. (See p. 509 infra.) 











| | 
Becillt per Saaient Corrected 
| Gell | alue Value 
0 1428 1402°75 
1 2633 2636°57 
2 3351 3355°55 
3 3556 3560°83 
4 2807 2810°81 | 
5 2038 2040°77 | 
6 1430 | 1431-94 | 
Y 1031 | 1032-40 
8s | 658 658°89 | 
9 407 —s | 407°55 
10 271 271°37 | 
11 163 | 16322 | 
12 | 103 | 103-14 | 
e 61 | 61:08 
1h 37s 37°05 
15 17 17°02 | 
16 9 | 901 
Seen arene 
| Totals | 20000 | 19999°95 


It will have been noticed that for the purposes of this count we did not enumerate cells 
which contained clumps of bacteria. Our reason for this omission was that, to the best of our 
knowledge, workers ordinarily disregard cells the contents of which cannot be resolved into 
definite bacilli or fragments of bacilli, and we were most anxious to use data so far as possible 
identical with those employed in actual practice. It was, however, subsequently pointed out to 
us, that this process was open to serious objection for the following reasons. In the first place, 
some workers of great reputation never employ for diagnostic purposes mixtures containing an 
appreciable number of clumps, say more than two or three per cent. In the second place, the 
exclusion of clumps would tend arbitrarily to heighten the proportion of cells containing 0 bacilli, 
since cells containing clumps must be regarded as possessing a phagocytic power at least equal 
to those which contain one bacillus apiece; indeed, Harvey and McKendrick have adduced 


evidencet which supports a belief that it is rather the number of bacteria than the number of 
acts of ingestion which should be taken as a measure of phagocytic power. 


* 1428 cells with 0 bacilli per cell will occur not in 20,000, but in 20,360; hence in 20,000 we shall 


have not 1428 but 1402-75, the excess is distributed proportionally among the other groups: see p. 509. 
+ Biometrika, Vol. vit. p. 64. 
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TABLE IL. 
3500! oa ce 
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Number of Bacilli. 
GraprH 1. Curve of 20,000 cells. 


Constants of Fitted Curve. 


o=2°6031, 8,=10365, B,=4:2622, xo=—1°6752, r—22°8160, b=41°5558, m,=1°3806, 
m= 19°4354, a=2°7561, a,=38°7997, %4%=3478°81, Skewness=6069, Mode at 21000, 
Mean at 3°6797, Start of curve at —°6562, End at 40°8996, y?=862, P very small, 


3478°81 { 1 . eo L 19-4354 
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Recognising the force of this objection, it was decided to examine the effect of including 
clumps in the count. To this end, another 2000 cells were counted on the slides already used 
and clumped cells were included. Of the 2000 cells counted in this way 36, or 1:8 per cent., 
contained clumped bacteria. Now if we allow, as seems reasonable, that this latter enumeration 
was a fair sample of the population, then the frequencies in our 20,000 cells would be changed 
in the manner shown in the second column of Table I. We have distributed the clumped cells 
among the other frequencies above zero in accordance with the proportions of those frequencies 
in the original count, since the small number of clumped cells which were present in the 2000 


did not allow us to determine the true proportion of clumps, whether mainly of two, three, 
or more bacilli*., 


The alteration which is effected by this correction is so trifling, that our erroneous method of 
counting, if it be considered an erroneous method, can hardly have been the cause of the poor 
fit which resulted, and we have not thought ourselves justified in discarding the data originally 
collected. Some other source of the poor fit must be looked for. 


Other possible explanations are: 


(1) The existence of heterogeneity in the material dependent on the fact that cells from 
different parts of the films are not, or may not be, strictly comparable with each other. 


(2) An artificial heterogeneity dependent on the process of counting, other than that already 
discussed. 


It has been noticed before that a large sample of material, even when adequately described 
by a frequency curve from the diagrammatic point of view, fails to satisfy the approved test. 
Indeed, many statisticians have been, as Pearson has remarked, far too easily satisfied with the 
test of mere inspection. 


Elderton writes: “I have found in applying the test, that when numbers dealt with are very 
large, the probability is often small, even though the curve appears to fit the statistics very 
closely. The explanation is that the statistics with which we deal in practice nearly always 
contain a certain amount of extraneous matter, and heterogeneity is concealed in a small 
experience by the roughness of the data. The increase in the number of cases observed removes 
the roughness, but the heterogeneity remains. The meaning, from the curve-fitting point of 
view, is that the experience is really made up of more than one frequency curve, but a certain 
curve, approximating to the one calculated, predominates t.” 


It will have been noticed that the poverty of fit is mainly due to the cells containing 
one bacillus being in defect and those containing three in excess; these two groups have added 
nearly 58 to the value of x2. Now the work of counting is excessively monotonous, and after 
going through some hundreds, it seems impossible to escape an impression that a certain 
measurement, say 3 bacilli per cell, is modal. Hence a tendency will arise to place any doubtfuls 
within that particular class. If the count be limited to one or two thousands, the error so 
introduced may not appreciably affect the result, but it will do so if the data run to many 
thousands, since the same percentage deviation in a large as in a small experience has naturally 
a much greater influence on the fit. 


* The 36 cells which contained ‘“‘ clumps ” were actually distributed as follows : 


Number of Number of Number of Number of Number of Number of 
Bacilli Cells Bacilli Cells Bacilli Cells 
1 2 4 4 7 5 
2 6 5 8 8 1 
3 8 6 0 9 2 


+ Frequency Curves and Correlation, by W. Palin Elderton, p. 142. 


Biometrika vir 65 
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The only prejudice to which the operator seems to have been—and thinks he was—subject 
was the idea which impressed itself on his mind during the count that the number of 3 bacilli 
per cell occurred with the greatest frequency; this idea may have inclined him to welcome 
threes with undue alacrity. It is possible that this may have helped in some degree towards 
producing an excess of threes in the final result. The operator believes that, at least in his 
case, the auditory sense is more impressed with the sound of numbers containing e’s and 7’s, 
3, 5, 9, than with others; he thus tended to believe that these numbers occurred more frequently 
than was actually the case. Evidently the truth of the matter cannot be exactly known until 
several independent observers have undertaken the same large count. 


We have already mentioned the difficulty regarding clumps; another point is as to whether 
some gross heterogeneity were produced either by notable differences in various slides or by 
variations in the method of counting on different occasions, differences of which the operator 
was unconscious. In making the frequency tables just considered, there emerged necessarily a 
rough representation of a frequency curve for each five hundred cells. 
a rough similarity and each one a rough regularity. 
appreciably less than in others. 


All these curves showed 
In some cases, however, the regularity was 
On going through the data a second time, the groups of 500 
which showed the greatest irregularities were excluded, on the supposition that they might 
represent a somewhat lower grade of accuracy in counting and might have introduced a gross 
heterogeneity into the total count. The application of this criterion led to the exclusion of 
5000 cells. The remaining 15,000 cells were then analysed as before. 


Tables III, III a and Graph 2 indicate the results of this process. There is a slight improve- 
ment in the fit, but the change is not pronounced. 


It can hardly be said that the results justify 
the exclusion. 


TABLE III. 
Actual Distribution of 15,000 Cells. 


Bacilli Number of Bacilli Number of 

per Cell Cells per Cell Cells 
0 1092 9 306 
1 2008 10 203 
2 2585 11 126 
3 2632 12 69 
4 2096 13 44 
5 1518 14 24 
6 1035 15 12 
7 760 16 6 
8 484 





Although we should naturally have preferred to obtain closer fits, the result 
does not seriously affect the value of our material for the purposes of the present 
investigation. It is clear that the heterogeneity is not caused by a compounding 
of Gaussian distributions, and that we have to deal with markedly skew variation. 
We sincerely trust that ome other workers may be tempted to repeat the experi- 
ment. We now turn to our main inquiry, viz. the problem of sampling. 


The 20,000 cells were arranged in consecutive 25’s, 50’s and 100’s, and grouped 
In each case the group is an inclusive one, thus 76—83 


as shown in Table IV. 
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By Bo Ke r b my, Mg ay Ay 
10381 4:2346 —1°5233 204319 38°1013 1:2784 1771535 26426 35°4587 
Yo Skewness Mode Mean Start End P 
2626°615 “61999 2°0360 3°6421 — 5533 37°4947 “00001 
a \12784 2 17-1535 
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TABLE IV. 
Actually Observed Frequencies of the Samples from 20,000 Cells. 


SampLes oF 25 Sampues or 50 Samputes or 100 
| 


Number of Number of Number of Number of Number of Number of 


Bacilli per Sample | Samples — Bacilli per Sample; Samples — Bacilli perSample | Samples 


4h— 51 l 11 23 2 267 —282 2 
5Y 1 


2—123 

g 59 14 24—135 9 283—298 9 
60— 67 38 136—147 21 299—314 4 
68— 75 74 148—-159 30 315—330 13 
716 89 128 160—171 58 331 —3B46 29 
84 91 145 172—185 83 347 — 362 36 
92— 99 147 184—-195 71 363—378 34 } 
100—107 113 196—207 64 3879—894 34 
108—115 74 208 219 31 395-—410 5 
116—123 37 220 —231 15 411—426 11 
124—131 19 232243 10 427 —J42 | 5 
182 —139 7 244 —255 3 443—458 3 
140—1 47 2 256—267 l 459-—h7 4 1 
148—155 0 268—279 0 475—4I90 l 
156—163 0 280— 291 l 491—506 2 
164---171 0 292303 l 507—-522 } 1 
172— 179 l 

Totals 800 400 — | 200 





includes samples with 76 and samples with 83 bacilli per 25 cells. For convenience 
of calculation, we took an arbitrary origin and unit of grouping in each case. In 
obtaining the moment coefficients Sheppard’s corrections were used in all cases; 
Table V gives the constants and equations, Graphs 3, 4, 5 illustrate the curves. 


TABLE V. 
Curves Resulting from an Analysis of the 20,000 Cells. 


. .| Meg i arms | } 2 (in terms 
Muuhae of [ean (in terms | Mode (in terms 





Material 2 2 of Bacilli per of Bacilli per | % (in nas By Be 
Samples Sample) Sample) of Bacilli) | 
; ae 
Samples of 25 800 92°04 89°7488 16°7184 "1033 -3°4321 
Samples of 50 109 184°07 181°3964 25°7496 "1467 | 4°0907 
Samples of LOO 200 367 °62 360°2632 39°9824 3004 8 4°3500 
Material Kg r m v a B.-3 | Skewness | P (Fit) | 
‘ P zs hy med es | 
Samples of 25 1436 25-2068 13°6033 |-10°3220 9°5149 -4321| °1371 ‘9765 | 
Samples of 50 0665 10°1436 6°0718 -— 2°7070 | 6-2691 1°0907, +1066 8673 


Samples of 100 "1367 10°1723. =6°0862 — 4°0483 | 7:0316 1°3500 "1840 1727 
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9 -13" 299) -_ a“ 
ae eae on 2 a im de 13°6034 eae ‘saa. 
4 y L” * 951492 


Origin at 2°6712, where 39°5 bacillj are taken as 0 and 8 bacilli per 25 cells form the statistical 


unit of grouping. 
x? —6°0718  2°7070tan-1 _” 
y=60°047 [ | xe SINE 





* 626912 
Origin at 4°87446, 105°5 bacilli taken as 0 and 12 bacilli per 50 cells as the statistical unit of 
grouping. 


L 
=18'117 iene 
Y si [2 + 70316 


Origin at 4:021598, 258°5 bacilli taken as 0 and 16 bacilli per 100 cells as the statistical 
unit of grouping. 


2 —7'0316 4°0483tan—1__* 
7-0316 | 
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Opsonic Index. 


Grapno 3. Curve of Samples of 25. 
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Opsonic Index. 


Grapu 4. Curve of 400 Samples of 50. 


The figures below the abscissa represent opsonic indices in terms of the mean of 
the whole series. That is to say, each group unit was translated back into physical 
terms and then divided by the mean number of bacilli per sample for the whole 
count. For instance, in Graph 3, the abscissa ‘9 was really 5 in statistical units 
corresponding to 39°5 (the origin) + 5 x 8 bacilli per 25 cells. The mean is 92°04, 
therefore the opsonic index in terms of the mean corresponding to this abscissa 
is kd = 86 *. 

92°04 


* The indices as plotted are only given to the nearest first decimal. 
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Objections may fairly be urged against this process. As we have previously remarked, the 
mean is theoretically not the best constant to use for descriptive purposes in the case of 
markedly skew frequencies and we suggested the mode as being more satisfactory. Harvey and 
McKendrick, in an important contribution to the subject*, have alluded to this point and remark 
with justice that the determination of the true mode requires calculation beyond the range of a 
laboratory worker. We have considered whether an approximation of sufficient accuracy might 
not be found. Pearson has shownt that in many cases a good enough value for the mode is 
given by: 

Mode — Median= Twice the Distance from Mean to Median. 


Unfortunately, however, the median cannot be determined with sufficient accuracy in the 
case of discontinuous variates such as phagocytic counts. Thus, we have experimentally no 
finer measurement than the integral number of bacilli per leucocyte ; if we mass in a group all 


the m individuals containing the same number of bacilli and know that the median is the nth 

a : : n ; : 

individual in that group, all we can do is to add — of the unit of grouping to the next lower 
m 


integral group measure. This is not nearly close enough in the case of such skew distributions ; 
and we found that the value of the mode thus obtained differed materially from the real mode. 
It was not therefore possible to use this process with success in the case of our sample curves. 
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Opsonic Index. 
GrapH 5. Curve of Samples of 100. 


* Biometrika, Vol. vit. p. 64. + Pearson, ibid. Vol. 1. p. 260. 
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We next treated, for purposes of comparison, the selected count of 15,000 cells 
in a similar fashion. Table VI contains the actual frequencies, Table VIA the 
constants and equations. Let us now consider in detail points suggested by these 
analytical results. The values of P show that in all cases the fit is close except in 
that of the sample hundreds. It appeared to us probable that in this case we had 
adopted too fine a unit of grouping, so that we re-worked, using a coarser unit 
which gave 11 groups, with the following result, 


we a \ 3-7098 tan-1 (x/4-4732) 
y = 29°97% = 2 
7 134\1 + paraay) *° , 


which gives P = ‘232. 

This finding confirms the suggestion to some extent, but it must be admitted 
that the 100’s are still definitely less satisfactory than the smaller samples; the 
explanation of this is not clear. Table VIB, which gives all the actual and 
calculated frequencies reduced to 2 common scale, will show the reader how little 
change was effected by the rejection of 5000 cells. 


TABLE VI. 


Actually Observed Frequencies of the Samples from 15,000 Cells. 


SamMpLeEs OF 25 SamMpies or 50 SamMpies or 100 


Number of Bacilli | Number of Number of Bacilli| Number of | Number of Bacilli | Number of 


per Sample Samples per Sample Samples per Sample Samples 
44— 61 l 112—123 1 267 --—282 2 
52— 59 11 124--135 8 283—298 8 
60— 67 29 136—-147 15 299 —314 3 
68— 75 62 48—159 25 315—330 10 
76— 83 98 160 —171 18 331—346 21 
84— 91 115 172—183 66 347-—362 27 
92— 99 107 184 195 50 363—378 30 
100—107 78 196—207 47 379—394 27 
108—115 56 208—219 22 395—410 10 
116—123 23 220—23. 9 411 —426 6 
124—131 15 232—243 6 427 —442 3 
132—139 4 24—255 l 443-—458 1 
140—147 1 256—267 l 45I—474 0 
268—279 0 475—490 0 
280—291 l 491—506 2 


Totals 600 — 300 — 150 
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TABLE VILa. 
Curves Resulting from an Analysis of 15,000 Selected Cells. 


} ’ o | 
| Mean (in terms | Mode (in terms 


Number of } 














Material “3 of Bacilli per of Bacilli per |% (in terms | | % 
| Samples Sample) Sample) | of Bacilli) | Pr Ba 
| ats or Tae 2 2 oe r aries : 
Samples of 25 | 600 | 91°1133 89°0117 | 16°3123 0457 | 2°8211 | 
Samples of 50 | 300 182-2200 179-4383 | 24°5984 | -0980 | 3-7994 . 
| Samples of 100 150 363 °5667 360°1649 | 37°0043 "1019 | 4°5418 
| : 
\ : | | 
: Material Ke r m v a B2-3 Skewness | P(Fit) 
( | 
H 2. are of e oes ened eae | 
} Samples of 50 +°0582 12°4217 | 7°2109 | —3°0885 | 6°7230 + ‘7994 1131 *7365 








| 

Samples of 25  ~-0702) — i — |— +1789) -1288 | -9966 | | 
) 

Samples of 100 | +-0292 | 7:4301) 2°7151 | —1-2890 | 5-7783 41°5418) -0920 | -2164 | 
| 


ee rey A x  \T0048 ' aw \12°5157 
uations : y=114°72 ss zu = 
‘a . y ( rq 1843) ( reas) 
Unit of grouping 8 bacilli per 25 cells. 39°5 bacilli per 25 cells taken as zero. 
a2 ]-72109 30885 tan-? 50 
y= 45°045)} 1 ——; xe 5°72 
y | * sia | ‘ 


Origin at 4°7217, where 12 bacilli per 50 cells form the unit of grouping and the zero is at 


105°5. 
a2 —47151 1-2890 tan-! __™ 
y=26°445] 1 + - xe 5°7783 


or" 783? 


Origin at 55642, where 16 bacilli per 100 cells are the unit and the zero is at 258°5. 


Next as to the agreement between the theoretical constants of the curves of 
means and their empirical values as above determined. Table VII gives the 
calculated values of B, and B,—3 [assuming that where n is the number of 
individuals in the sample and £,, 8, the constants of the original distribution 
B,=% and B,-3=%—*) together with the actual values. 

In only one case, that of the 25’s from the selected 15,000, is the agreement at 
all possible. “Student,” in a very interesting communication*, has suggested that 
the divergence noted by him in our previous curve of 80 means, depended on the 
existence of homotyposis in the samples. This suggestion no doubt contains some 
truth and could be satisfactorily tested by mixing up our individual records of 
leucocytes and drawing really random samples. We have not felt justified in 
delaying publication of our results until we could find leisure to undertake this 
very laborious task, particularly as other factors are involved. This is clear from 
the fact that the divergence is more marked in the case of the larger samples 
which must be less influenced by homotyposis. 


* Biometrika, Vol. vu. p. 210. 
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TABLE VII. 


) — es 
Number in Sample B, (empirical) | B, (caleulated) |B —3 (empirical) | By —3 (calculated) 
| | 


| (20,000 Count) 


| 25 hi e "1033 0415 "4321 "0505 
50 =a "1467 “0208 1-0907 “0252 
100 (ist curve) ... | 3005 0104 1°3501 0126 
100 (2nd curve)... 3142 0104 | 14601 0126 

| (15,000 Count) 
25 ‘ we “0457 “0415 — °1789 0494 
50 = exe | “0980 “0208 “7994 “0247 
100 na a | 1019 0104 | 15418 0124 


We can now examine the accuracy of the opsonic method in the way attempted 
in our former paper, viz. we can inquire what the chances are of obtaining from a 
“population” of means, samples giving different indices in terms of the real mean 
of all such samples. .In similar cases, it is usual to divide the frequency curve 
into a series of equal areas, setting up ordinates on either side of the mean or 
mode. Since this elaborate graphical method requires the services of a specially 
skilled draughtsman, and considerably increases the cost of production, we decided 
to adopt the plan followed in our last paper. A considerable number of additional 
ordinates were calculated for each sample curve and we determined the area from 
the beginning or end of the curve up to specified ordinates which corresponded to 
indices of ‘6, ‘7, ‘8, etc. in terms of the mean. From these areas the Tables of 


Chances (VIII, IX, X) were deduced. 


TABLE VIII. 
Chances of Obtaining Certain Deviations in Index Values. 


j 
SaMPLEs or 25 


Opsonic Index in 
terms of the 


athe Fraction of the Total Odds against 
es Area bounded by the | the occurrence of 
corresponding ordi- such a Deviation 
nate (Total Area=1) | or a Greater 
| 
6 “0089 111 tol 
— 0396 24°3 to | 
8 "1276 6°8 to 1 
9 *2973 2°4 to l 
Yi *2740 2°7 to 1 
1°2 *1324 6°6 to 1 
1°3 “0558 16°9 to 1 
1*4 ‘0211 46°4 to 1 
Beyond the limits 1-4—°6 0300 32°3 to 1 
= = 1°3—-7 "0954 9°5 to 1 
* - 1°2—°8 *2600 2°9 to l 
” a 1-1—'9 5713 ‘75 to 1 
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TABLE IX. 


Chances of Obtaining Certain Deviations in Index Values. 


Samp.es or 50 


Opsonic Index in 


| terms of the Fraction of the Total Odds against 
Mean Area bounded by the | the occurrence of | 
corresponding ordi- such a Deviation | 
nate (Total Area=1) or a Greater | 
| | 
7 ‘0112 88°3 to 1 
| 8 ‘0618 15°2 to 1 
°9 2283 3°4 to 1 
Ez ‘2160 3°6 to 1 
1:2 ‘O770 12 tol 
1°3 0248 39°3 to 1 
3eyond the limits 1°3—-7 “0360 26°8 to 1 
ef a 1:2—8 ‘1388 62 to 1 
“ - 1‘1—°9 "4443 1°3 to 1 


TABLE X. 


Chances of Obtaining Certain Deviations in Index Values. 


Sampies or 100 


Opsonic Index in 


terms of the Fraction of the Total | Odds against 
Mean Area bounded by the | the occurrence of 
corresponding ordi such a Deviation 
nate (Total Area = 1) or a Greater 
“” “0006 1666 to 1 
‘8 "0216 45°4 to 1 
9 ‘1661 50 to | 
rhe "1618 5°2 to 1 
1°2 0420 22°8 to 1 
1°3 “0109 90°8 to 1 
Beyond the limits 1°3—-7 *O115 86°6 to 1 
fe “i 1°2—°8 0636 14:7 to 1 
a = 1°1—°9 *3279 2°1 to 1 


If we adopt about a ten to one chance as the limiting value for evidence of 
differentiation, then the limits, for single determinations of the index, are in the 
ease of 25’s, roughly 13—7; for 50's, 1°25—75; for 100’s, 1.2—85. Everyone 





has his own standard of accuracy or reliability ; all we mean is that unless a given 
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index is beyond these limits, the chances are less than ten to one against its being 
drawn from a normal population. 


We now come to a most important qualification 
of the value of these tables. 


In the table as given, the measure is of the probable deviation of a sample 
from the real mean of the samples as determined from a long count. If the control 
mean, i.e. the denominator of the fraction which we call the opsonic index, has 
been determined from a large number of cells, say 1000, and the tested serum 
is measured on the basis of a small sample, then our results are directly applicable ; 
but this supposition does not represent the usual state of affairs. In ordinary 
practice, two small samples, each of not more than 100 cells, serve as control and 
test, there is therefore, so to speak, a double source of variation; odds based 
on the results obtained above would over-estimate the reliability of any given 
determination. This is a point of considerable importance which has, we think, 
escaped the attention of writers on the subject; its elucidation would involve 
a careful examination of the whole theory of index distributions, and would 


TABLE XI. 


Constants cf the Experimental Index Frequencies. 


z , 
| | | 
Samples Mean | Mode | o Ai | Bo ° 








Kg rT 
°5’s 800 1°0362 9432 | °2585 ‘7099 | 4°3914 “9619 — 
50’s* 400 1°0205 “9536 | *2204 "6952 4°8945 | °3635 | 11°2690 | 
100’s_ | 200 10120 | 9689 | °1553 3654 | 3°7138 | 90382 | — | 
| | | | 
Material m v a B2-3 |Skewness y Pp | P (Fit) 
25’s oe — — | 1:3914 3600 | 326°1439 27°4967 | °0853 
50’s 6°6345 | —8°5166  5°6343 | 1°8945 2912 .- — | °8953 
100’s — — — *7138 *2780 491°3142 48°7621 “9355 
Fe. 
Equations : y=antilog 43°6116 Fa xe ). 


A change of *1 in the index is taken as the statistical unit of grouping and index ‘25 as zero, 


5°6343 


, x2 —6°6345 85166 tan~! _* 
y =6°2092 (1 +-~ xe 
y 5°6343? 


A change of *1 in the index is taken as the statistical unit of grouping and index -45 as zero. 
The origin is at 1°4468. 
491°3142 
: eet sore = 
y=antilog 71°8335 (. eG oe i 


A change of ‘1 in the index is taken as the statistical unit of grouping and index °45 as zero. 


* A slight slip has been found in the values of 8; and f, for this series, they should be -6836 and 
4°9855. This error does not, however, visibly alter the form of the curve in Graph 7. 
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accordingly present features of interest quite apart from the particular case of 
the opsonic index. These remarks apply with special force to medical statistics 
which so frequently involve the use of ratios or indices. 


The importance of the problem tempted us to devote some little time to a 
theoretical consideration of the subject and we arrived at a few analytical results 
which might, we thought, be of some value. It is, however, unnecessary for 
us to publish these incomplete investigations, since Professor Karl Pearson has 
succeeded in obtaining a full solution of the problem*. We therefore confine 
ourselves to an empirical solution which is sufficient for our immediate object. 
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Opsonic Index. 


Grapu 6. Curve of 25’s secondary without replacement. 


* See the following paper. 


RETIRE 


ERE 

















523 


M. GREENWOOD AND J. D. C. Wurrs 


‘wapuy auosdQ 
vel St Gt bk 


*payjop—juemteovpder ynoyyTAs :snonulyu0o—juemteoe[der YIM 8G JO SAAINO 


*), Havay 














2 apow 


1 apow 





ee | 


- OL 


0% 


og 


- Ov 


og 


09 


-OL 


os 





*houanbaag 











524 On the Distribution of the Means of Samples 


The actual number of bacteria in each sample of 25, 50 or 100 cells was copied 
out on a separate slip of cardboard, so that we had 800, 400 and 200 cards 
corresponding to the samples of 25, 50 and 100 cells. The cards used were small 
pieces of correspondence cards which, being smooth, could be shuffled very easily. 


With these cards the following experiments were performed. Taking the case 
of the 400 cards relating to the fifties, the cards were thoroughly shuffled and one 
drawn out. The number on the card was written down, the card not being 
replaced. The remainder were once more shuffled, another drawing was made and 
the process was repeated until the whole 400 had been drawn. The cards were 
then once more shuffled and the process repeated, the second set of numbers beir z 
written down one by one under the first set of values. Each number of the first 
set was then divided by the corresponding number in the second set and the 400 


— 


Frequency. 


_ Mode 
Mean 


— — -— oe = =. lhe 2 wee | Mee anne 
Opsonic Index. 


Grappa 8. Curve of 100’s secondary without replacement. 
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indices thus obtained were analysed as usual. The same experiment was carried 


out in the case of 25’s and 100’s. Tables XI, XI a and Graphs 6, 7 and 8, contain 
the constants and data of the distributions and Table XII the chances of obtaining ij 
indices within assigned limits deduced from the areas. 











TABLE XTa. 
Experimental Determination of Index Distributions. 
SaMPLES oF 25 Sampues or 50 SamMpues or 50 Sampres or 100 
(Drawings not Replaced) | (Drawings not Replaced) | (Drawings Replaced) | (Drawings not Replaced) | 
Index | 
Se a aa ee |r aco eee 
| | 
| | Observed | Calculated Observed | Calculated | Observed Calculated | Observed | Calculated 
| | s— +4 |) (oe f ie iE onde 
‘4— 5 3 frase! asi me a 4 | 
| — 6 16°5 11°243 | (3-989 3 3:007 | Be! 
6— “7 45 40°468 13-990 | 15°5 15308 2 1 1°500 
T— °8 83°5 89°850 37°969 | 34:5 | 39°863 | 10% 11°697 
8— 9 | 133°5 | 121-939 66°729 69 65°497 | 35°5 36-009 
‘9—1°0 | 130 134-204 80-213 71 77°683 | 55°5 52-780 
1°0—1°1 | 105 122-028 71-759 78°5 72001 | 41°5 46°583 | 
1-1—1°2 83 96°635 51973 | 59°5 54817 | 34 29°003 | 
1°2—1°3 55°5 69°188 32°676 | 28 35°425 | 13 |. “Srey 4 
1°3—1°4 57°5 46°012 18°792 | 22 19°873 4 5°850 |, 
ly—15 37 28-990 10-249 | 13 9°330 | 2 2145 | 
1-5—1°6 19°5 17°560 5452 3 4°334 | 2 { 1:00 
£6—17 11°5 10°339 2-870 | 1°5 1°717 | 
17-18 | 10° | 5-967 52 | 15 814 | | 
1-38—1°9 4 3°397 803 ) — 
1°9—2°0 2 1918 *432 >1°735 _ 
2-0—2°1 2 | (20 6-5 | — | | 
2*1—2°2 _ , — 1 ( _- | 
+293 na , a on ; | 
Totals 800 803°265 400 399-908 400 400°169 200 203°744 





A comparison of these tables with nos, VIII—X indicates how far the limits of 


reliability are narrowed by a more complete study of actual conditions (see also 


Table XV infra). 


It will be noticed that these experiments do not assume any special distribution 
for the original variates. They merely afford a partial answer to this question—if 
we had two sets of 400—200, 800—samples exactly the same as the 400 (200, 800) 
actually counted and formed at random 400 (200, 800) indices from them, one set 
always furnishing the denominator and no sample being used more than once, how 
are the resulting indices likely to be distributed? In strictness, therefore, these 
results are only applicable to our particular case. If, however, we assume that our 
counts are numerous enough and accurate enough to give the distributions of 
sample means of the size used, and for the emulsion used, with fair exactitude, 
then we can perform another experiment. 
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Let us suppose that the distribution of samples of 50 that we obtained is 
actually the distribution of samples of 50 in general, then the following experiment 
may be considered. The 400 cards were shuffied and drawn as before except that 
after each drawing the card was replaced, the process being carried on until 400 
numbers had been written down, a second series of 400 numbers was then obtained 
in the same way and the indices were computed. Tables XIII, XIV and Graph 7 








TABLE XIII. 

Constants of the Experimental Frequency Distribution Obtained by the Method of 
Replacement. (Mean, Mode and o, Expressed in Terms of the Opsonic Indez.) 

Mean Mode o Bi Bo r b 


my mys 


1°0196 “9697 *2056 2077 3°2155 62°7551 44°6934 9°0844 51°6706 


ay ag kK Bo-3 Skewness r 
6°7337 37°9597 — 8536 °2155 "2428 *8988 
" fA av \9°0844 : Ps 51°6706 
ation : =80°0808 —— — oemrar : 
Equation y ( a oan) ( sraam) 


A change of ‘1 in the index is taken as the statistical unit of grouping and index °45 as zero. 


TABLE XIV. 


Chances of Obtaining Certain Deviations, on the Basis 
of an Experimental Curve from Samples of 50 
(Method of Replacement). 


| Fraction of Total | Gite wsiliek 


Opsonic Area bounded by | oti 5 iis 
Index corresponding ‘ey or a Gleuiae 
| Ordinate ; 
6 | 0075 | 132 «to 1 
5 | 0458 20°8 tol 
‘8 | "1454 5°88 to 1 
9 “3096 2°23 to 1 
rez “3165 2°16 to 1 
1°2 "1795 4°57 to l } 
1°3 ‘0910 999 tol | 
14 0413 23°2 tol } 
1°5 ‘0167 58°9 to l 
1°6 *0059 168 to 1 


| 


give the statistical constants and other information which resulted from an analysis 
of this experiment. It will be noticed that the table of chances does not differ 
appreciably from that due to the other experiment (XV). On this account, and 
because of the doubtfulness of the assumption in the case of 100’s and the smaller 
practical interest in the case of 25’s, the experiment was not repeated in these 
cases. 
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TABLE XV. 


Comparison of Odds against Certain Index Deviations as Calculated 
on Different Assumptions for Samples of 50. 


indies Using True Experimental Curve Experimental Curve 

, Mean without Replacement | with Replacement 
s 883 tol | 213 tol 20°8 tol 
8 162 tol | 616 to 1 5°88 to 1 
9 3°38 to 1 | 2°26 to 1 2°23 to 1 
ts 3°63 to 1 2°19 to 1 2°16 to 1 
1°2 12°0 tol 4°45 to 1 4°57 to 1 
1°3 39°3 to 1 8°81 to 1 9°99 to 1 


Without in the least wishing to underestimate the shortcomings of both our 
material and method, it appears to us that Table XII can fairly be used as a 
rough practical test of the accuracy of the opsonic method in the case of the 
Tubercle Bacillus and that it confirms the opinions expressed in our former paper. 
A theoretical discussion of the index problem in general is reserved for another 
paper and more suitable material. It is now proper to consider very briefly some 
practical points. 


In the first place, it has been asserted that the Tubercle Bacillus is an 
unfavourable organism to select for the purpose of studying the accuracy of the 
opsonic technic. While we cannot assert that no more favourable example could 
have been chosen, we are of opinion that the Tubercle Bacillus is certainly not 
the least satisfactory of the organisms commonly dealt with. For instance, Tubercle 
Bacilli act, or are acted on, almost entirely as single individuals, while Staphylococci 
or Streptococci occur in groups or chains respectively. These groups or chains are 
counted as so many individuals which, from the point of view of phagocytosis, may 
or may not be correct. The opsonic power required to attract a group of cocci to 
the leucocyte may not be greater than that associated with the ingestion of a 
single coccus. Under such conditions we have the difficulties as to the correct 
treatment of doublets and triplets—to which allusion has been made by Harvey 
and McKendrick—in an acute form. 


Next as to the practical accuracy of the method for diagnostic purposes or in 
the control of dosage. The results here communicated, which are in agreement 
with our previous observations and also, mutatis mutandis, with those of Harvey 
and McKendrick, seem to demonstrate that single observations falling within the 
limits of 20°/, on either side of the mean cannot, in the case of Tubercle Bacilli, 
be regarded as strongly suggestive of abnormality. 


If it is possible, as will often 
be the case when it is a question of diagnosis alone, to repeat the observation 
several times under comparable conditions, the opsonic method may be expected 
to yield results of a certain value. When, however, the opsonic index is used to 
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control injections of Tuberculin, so that no two determinations are made under 
the same conditions, we do not think the method has more than a very restricted 
field of usefulness. This leads us on to refer to a cognate matter. 


Glynn and Cox use language* which seems to imply that we have attacked the 
opsonic index in the sense of denying that differences in phagocytic power as 
between normal and diseased persons occur. We have not asserted this nor does 
any statement of ours bear such an interpretation. We have merely pointed 
out, a contention not seriously challenged, that the claims advanced by certain 
enthusiasts to the effect that all variations of 10°/, from the normal are probably 
or necessarily significant, must be entirely unjustifiable. 


At the risk of being tedious, we desire to emphasise once more the distinction 
between vaccine therapy on the one hand and opsonic estimations on the other. 
If our work tends, as it certainly may tend, to unsettle the opinion of some readers 
as to the relative accuracy of the tuberculo-opsonic index, we wish in no way to be 
ranked among detractors of the claim of vaccine therapy to be regarded as a new 
and powerful weapon in the struggle against disease. On the contrary we are 


convinced of its present importance and look to the future for its still greater 
development. 


In conclusion we would again remark that the results here published are in 


essential agreement with those communicated by Harvey and McKendrick and 
ourselves in recent numbers of Biometrika. 


Conclusions. 


1. Large phagocytic counts extending to 15,000 or 20,000 cells are markedly 
skew and exhibit signs of heterogeneity. 


2. The exclusion of clumped bacilli has not been the main cause of the poor 


fit obtained. 

3. The heterogeneity may be due to difficulties in counting. 

4. The distributions of the means of samples of 25, 50 and 100 cells drawn 
from these populations are also markedly skew and excellently represented by 
Pearson’s type curves. 
5. This skewness cannot be due entirely to homotyposis. 

6. The odds against a second sample of the same material giving an opsonic 
index, in terms of the true mean, outside of the limits 1:2 to ‘8 are, in the case of 
samples of 25, about 3 to 1, in the case of 50’s about 6 to 1 and for 100’s about 
15 to 1. 

7. When the index is estimated on the basis of two small samples and without 
determining the mean of the normal series from a long count, the odds are reduced 
to 1°5 to 1, 2 to 1 and 4°6 to 1 respectively, for samples of 25, 50 or 100. 


* Glynn and Cox: Biochemical Journal, Vol. tv. p. 300. 
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8. Even in the case of samples of 100, it is doubtful whether single determi- 
nations which give values even beyond the limits 13 to ‘7 can be regarded as 
satisfactory evidence of differentiation. 


9. The two previous conclusions are based on the analysis of four experimental 
frequencies. 


19. The Tubercle Bacillus is not necessarily an unfavourable test of the 
opsonic method. 


11. The opsonic index is less useful as a means of controlling treatment than 
for diagnostic purposes, since in the latter case it may be possible to make several 
comparable determinations. 


12. Differentiation in respect of phagocytic power (Tubercle Bacilli) can 
possibly be established by the opsonic index in some cases. 

















ON THE CONSTANTS OF INDEX-DISTRIBUTIONS AS 
DEDUCED FROM THE LIKE CONSTANTS FOR THE 
COMPONENTS OF THE RATIO, WITH SPECIAL REFER- 
ENCE TO THE OPSONIC INDEX. 


By KARL PEARSON, F.R.S. 


(1) GIvEN a and y two variables, the frequency constants of both of which 


| are known, we often require in statistics the frequency constants of their ratio or 
4 index: 1=a/y. 


If the coefficients of variation are small, we have with the usual notation for 
means, standard deviations, moment coefficients, etc.*: 


x . , Pu re 
t= (=) {1 + Ux? — TayV2Vy — ee an 4G . ey etc 











oe eS ee eres (1), 

a\?2 ( 4p,. 2, Zs 

M,=0/7= (=) i + 0,7 — 2g %y Try + £ a - = a 

3p — Pu’, Boy — 2p? | 2po'Pn =. 
+ 2 4+ Ht BPG P cocrcceccececseses (11), 
vy Y “LY 

\3 (iu. are. 3Di2 . DPeo — 6py? — Bue pe 

aty= (2) {es — Hs — BP Be Pn — Spat — Bat 
— 3(Pa— MePu) _ 9 (Pis— Ha Pu) , 3 (md — Ha") ete} -.-(iii), 

“ry xy y* 
= (Z\* | pe be 6p»  4ps 4p: 5 ) . 

M, = (5) {es + 7 + By = ay - ay + ee aute ke veh baeekeenunseaeseae (iv). 


These formulae go to the order of the 4th power in the coefficients of variation, 
but of course this is not to the same order of approximation in M,, M, and M,. 


(2) It will be seen at once that these approximate formulae would be prac- 


tically unworkable if 2 and y were correlated, as we should have to find 3rd and 4th 
order product moments. 


* A rule denotes a mean value, o a standard deviation, v, =c/mean, is a coefficient of variation ; 


M2, M3, My are the moment coefficients for x, se’, us’, wy’ for y, Mz, M3, My for the index i, and 
Puy=S (x - Z)“ (y - y)"/N, where N is the total number of pairs. Thus p39=y3, Poo=se’, etc. 
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If they be uncorrelated we have : 

kets f a He Me | * \bis 
i= (=) ! + Uz y + 7 € tc. adnan Se Waleed eee Te seine eneeten ay, 

z\2 Qu. Sn’ ~ Sy” 

Doe a Oa 2 So . ees Oe tes ie Ss ste 5) \bis 
M,=02= (=) {v= + 0,7 — F + 30,70,7? + 7 7 ete veccevenones (ii), 

/H\3 : Eg Cu. my 2 ae 2 CR 
M,= \; ) (2 —& 4 Se = is ss Ha’) + et Buse hid tia Salen eee (iii), 
y |e Pf wy y ) 


Gstudecerenan vada e eruiocvaakeesie eet 


Formulae (i) to (iv) are extensions of formulae given by me in a paper on 
spurious correlation*. Formulae (i) to (iv) are due to Dr M. Greenwood, Jun., 
who obtained them in dealing with the problem of the distribution of the opsonic 
index. They show at once two noteworthy but not yet noted points, namely, 
(a) if the distribution of both # and y be symmetrical, i.e. w; and p,’=0, M; will 
not be zero or the distribution of indices must be skew; () the mean of the ratio 
of two numbers picked out of the same series is certainly greater than unity if the 
series be symmetrical, and will probably be always greater than unity even if it be 
not, ie. 7>%/¥, which is unity for the same # and y seriest. Dr Greenwood 
found, however, that these formulae did not give with sufficient accuracy the 
constants of the index distribution. This was probably due to two causes: 
(a) clearly we ought only to keep to the square order in M, and the cubic order in 
M, if we retain only to the 4th order in M,; or if we keep to the higher terms 
in M, and M;, we must go further with M,; and (0) the values of v, or v, are not 
so small, that the convergency is sufficient when we take these lowest terms of the 
expansions. It seemed accordingly desirable to find some other way of attacking 
the problem, and Dr Greenwood asked me for suggestions. The problem he had 
in view was the distribution of the opsonic index when the blood of the same 
individual taken in the same manner at the same time was treated as test and as 
normal. If a wide range of values of the opsonic index could thus be obtained, it 
would cast some light on what deviations from unity must be looked upon as 
significant, when test and normal were different individuals. 


x\4 "6 ope 
M,= (2) {B+ Mt + Sea 4 ote.} 
y ay 


\e yx 


(3) The idea that suggested itself to me was a fairly simple one, namely to 
tabulate the y-frequencies to a variable z=1/y. The units of the z-frequency 
groups will not be equal, but they will all be sufficiently small for us to concentrate 
their frequencies at their mid-points. We can then calculate their moments easily. 
Let v,, v2, ¥3, vy, be the moments of # about the zero value of #, and v,’, 1’, vs’, vy be 
the moments of z about the zero value of z. Then 

t= 2X o, 


* R.S. Proc. Vol. tx. p. 492. 


+ Given two dice, it would be advantageous to bet that the ratio of the number of pips on the two at 
a cast will exceed unity. 
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and if m,, m., m3, m, be the moments of 7 about its zero value, we have: 
m=Py, m=Pe, m=Ps, m= Py, 
where P,, is the wth product moment of z and «# about axes through their zero 
values. In the special case in which z and w are uncorrelated, as in the opsonic 
index, 
M, =VyV;, Mz=VeVe, Ms=VzyVz, My =%M%, 


or we can obtain any moment about the zero of 7 by multiplying the corresponding 
moments of # and z. 


These moments are then transferred by the usual formulae 

M,=m,—m,- M,=m,—3m,m,+ 2m, M,= m,— 4m;m, + 6mm? — 3ms4 
to the mean as origin and the type of frequency calculated in the usual way from 
the corresponding £, and f,. 

In Greenwood and White’s data we have, for the three series discussed below, 
elementary subranges rising by 32, ‘24 and ‘16 of a bacillus per leucocyte for 
the distribution of the means of counts of 25, 50 and 100. I find that for such 
distributions, the value of the variate z will only be affected by about a unit in 
the third place of decimals in the worst cases, i.e. the lowest values of y in samples 
of 25, whether we use for z (i) the mean of the inverses of the start and finish 
of the subrange, (ii) the mean of the inverses of all the 32, 24, or 16 hundredths 
in the subrange, or (iii) the inverse of the mid-point of the subrange. I have 
accordingly adopted the last as the simplest value of z for practical purposes. 

(4) Illustration of the method. I. Greenwood and White's 200 samples of 
100 counts. 

The data are given in Table I. The moments of the frequency distributions 
for y and z as variates about the zero of those variates were then found by tables 
of powers of numbers and a calculating machine*. 


I. Distribution of 40,000 indices for 200 samples of 100 counts. 
For #: v, =3°67620, v», =13°67643, v, =51°50333, v, = 196°40357. 
For z:  v,{=0°275165, v./=0°076603, v. =0°021577, vv, = 0°0061494. 

These give: 
m, =101156, m,=1:04766, m,=1:11126, m,=120776, 


which transferred to the mean give for frequency constants of the 40,000 possible 
indices : 


by = 02440, o = ‘1562, « =°2651 
ps, = 00218, 8, = 3123, Mean = 1°01156, 
py = 00235, B.= 3°9472, Mode = ‘9774 


The distribution is accordingly of Type IV. 


* I have cordially to acknowledge help from Alice Lee, D.Se., Julia Bell, M.A., and Amy Barrington, 
who have each worked out nearly the whole of one distribution for me, and from H. Gertrude Jones, 
who has prepared the diagrams. 
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We find: r=16°5111, vy = — 9°9138 (since ps is positive), 
| = 52743, Yo = 8638°6, 


and the distribution of the 40,000 indices possible is given by: 


ro pelea sca ee 
a 35) —9°2556 9-9138 tan-! 52743 (i) 
‘27818 : 


Y =8638°6 (1 4 


where Y is the frequency at distance X from the origin of the curve which is at 
the opsonic index *6949, the unit of X being absolute opsonic index measurement. 
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It will be seen, in accordance with the results drawn from the approximate 
formulae, that the mean value of the opsonic index is greater than unity and that 
there is marked skewness, i.e. no approach even with samples of 100 to a Gaussian 
distribution. Clearly the most probable value of the index is less than unity, and 
is rather farther from unity than the mean. 


Greenwood and White have made an experimental determination by 200 
random drawings. In Diagram I the frequencies of the above curve are reduced 
to 1000 total and plotted against their experimental data increased to 1000. This 
of course exaggerates the apparent deviations, but enables our three diagrams to 
be compared among themselves. 


(5) Further Illustrations. The following frequencies are also deduced from 
Greenwood and White’s results for 400 samples of 50 and 800 samples of 25. 


With the same notation as before the constants of the distribution of the 
opsonic indices are as follows: 
Il. Distribution of 160,000 indices for 400 samples of 50 counts. 
For #: », =3°68140, v, =13°82272, v, =52°92733, v, = 206°69803. 
For z: » ='277213, / ='078453, 4’ = "022683, _v, = 006705, 


m, = 1:02053, m,=1°'08443, m,=1°'20055, m,=1°'38591, 


fo = "042948, 9 o ='2072, « =°5415, 
v3='006191, 2, ='4837, Mean = 1°0205, 
fy = "007567, 8,= 41022, Mode = ‘9611. 


The distribution is of Type IV: 


r = 20°86018, vp = — 2267044, 
a= '62536, Yo = 169740, 


and the equation to distribution of the 160,000 indices possible is given by: 


; X 
- tas X?_ \ —11-43009 22-67044 tan-! ——_ td 
Y = 16°9740 1+ ‘39107 é *62536 (ii). 


The origin of the curve is at the opsonic index *3409. 


Diagram II gives this theoretical distribution reduced to 1000 cases and 


compared with Greenwood and White’s experimental data for 400 increased to 
1000 also. 
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III. Distribution of 640,000 indices for 800 samples of 25 counts. 
For «: v, =368160, v, =14:00990, v, =55°03056, v, = 22290230, 
For z: v, ='281298, vy.’ ="082055, _—v,’ = 024861, _—v,/ = 007834, 

m, = 103563, m,=114959, m,;=1°36811, m,=1°74622, 


2='077063, o='2776, « = 1°2480, 
ps; = 017943, B,=°7060, Mean = 1:0356, 
Py = "025583, B= 43088, Mode = 9330. 
The distribution is of Type VI*: 
r = — 31:25271, qh = 51°67822, 92 = 18°42551, 
a = 152075, log y, = 26°6813462, 


and the equation to the distribution of the 640,000 indices possible is given by: 
Y =antilog 266813462 x (X — 1:52075)##51 X-nviem 
The origin is at the opsonic index — 1:4304. 


Diagram III gives this distribution reduced to a total of 1000 indices and set 
against Greenwood and White’s experimental curve for 800 increased to 1000. 
This exaggerates the apparent deviations, but enables us to compare with the 
samples of 50 and 100. Considering the relative paucity of Greenwood and White’s 
drawings, and the fact that we have fitted our curves to their non-replaced material, 
whereas a random sample would probably be better represented by the replaced 
material, the fits seem fairly close. The actual goodness of fit test was not applied, 
as at the time of writing the paper means of drawing the curves on a large scale 
and mechanically integrating them were not accessible. 


(6) It is as well to look at a combined table of the frequency constants of 
the three distributions given above. 


No. of Counts | No. in Distribution 


average based on of Means | Mean Mode | ’ By | By 
a — . —— — ian 
25 | 800 1°0356 | ‘9330 "2776 7060 4°3088 
50 400 1°0205 | ‘9611 *2072 "4837 4°1022 
100 200 1°0116 ‘9774 "1562 *3123 3°9472 


| 


It will be seen that these values form a very consistent relatively smoothly 
altering system. But the approach to normality is very slow. Even with counts 
of 100 the distribution of the opsonic index is markedly skew and platykurtic, and 
it would not be safe to treat the distribution as a normal curvet+. In all cases the 

* Type VI is adjacent to Type IV and in this case the curve is almost on the boundary line— 
i.e. Type V: see Rhind’s Diagram, Biometrika, Vol. vi. p. 389. 

+ The reduction of 8; and 82 to the Gaussian values 0 and 3 is, of course, not at the same rate as if 
we averaged 25, 50 and 100 opsonic indices. What we are doing here is to average the number of 


bacilli in 25, 50 and 100 cells on which the two factors of the index are based—a very different 
process. 
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modal or most probable value of the opsonic index of an individual tested by 
himself is less than unity and his average value greater than unity. 


If we actually suppose the distributions normal and varying round the means 
with the standard deviations given above we have the following results: 





aa , Number of Counts 
Significance laid on an event 


| which will not happen more |- 
frequently than 


| 25 50 100 

eas paws 

| Once inten trials ... | S8—1-48 | -68—1:36 | -75—1-27 
Once in eight trials... 61—1°46 *70—1°34 ‘77—1°25 
Once in six trials aes *65—1°42 ‘73—1°31 *80—1°23 


This table is to be read in the following sense: If an opsonic index were based 
on a count of 50 cells, then once in eight trials an individual tested against himself 
would have an index lying outside the limits ‘70 to 1:34. Or, again, with 100 
counts once in six trials an individual tested against himself would have an index 
lying outside the limits *80 to 1:23. It may be added that for most of the purposes 
of practical life or of exact science we should not consider an “improbability ” 
which could happen once in ten trials as marking a significant differentiation. 
Much greater degrees of improbability would be required. In the case of medicine, 
however, much less certainty may be demanded of a judgment, and probably no 
weight would be given to an isolated opsonic determination in diagnosis. Still 
‘the matter gives ground for pause, the opsonic index of the same material tested 
by itself has a very wide range round unity. 


The following table will show how closely Greenwood and White’s experimental 
determinations agree with our theoretical evaluation of the constants : 


Samples of 25 Samples of 50 Samples of 100 
| Constant of | 
Distribution feria (eae eo ee 
| | G.+W. P. | G.+W. P. G.+ W.* | P. 
Mean Ue a | 10362 | 1:0356 | 10205 10205 10120 10116 
Mode nile ons | 9432 | ‘9330 9536 ‘9611 “9689 ‘9774 
| Standard Deviation "2585 ‘2776 2204 2072 1553 1562 
By ‘7099 | ‘7060 “6952 4837 3654 3132 
By 4°3914 4°3088 4°8945 4°1022 3°7138 3°9472 


| 
| 

The only substantial divergences are in Greenwood and White’s values for f, 
and 8, in the case of samples of 50. I have been through my results again and 
can find no error. Dr Greenwood has been through his figures and finds slight 


* They have used the boundary curve Type V in the case of samples of 100. 














KARL PEARSON 541 


slips (see footnote, p. 521) in the values of 8, and 8,. These are not sufficient to 
account for the divergence of their results from mine, and from the trend of their 
samples of 25 and 100. I am-inclined to think the divergence is due to the 
presence of the outlying index 2°16 in their chance drawings—a result my curve 
shows to be exceedingly improbable not once in 3000 drawings, and the like of 
which does not occur in their samples of 50 (drawings replaced) (see p. 525). It 
would clearly be possible from their experimental method of drawing samples 
of 800, 400 and 200 indices to closely approach the corresponding theoretical 
distributions of 640,000, 160,000 and 40,000 represented in our curves. Our 
results then are in every way confirmatory of theirs, but place on a rather more 
satisfactory theoretical instead of experimental footing their deduction of index 
distributions. The general conclusion seems to be that except in the case of 
an extremely low or extremely high value of the opsonic index, little if any 
weight whatever ought to be placed on a single determination of this index. 
Hence the method would not be valid when applied to cases in which, owing 
to the evolution of a morbid process or the action of some drug, very few observa- 
tions can be made under the same conditions, i.e. it must be of doubtful application 
in treatment. Further the concentration obtained by basing the index on a count 
of 100, rather than one of 50 or even 25, while sensible is not very rapid. It 
would require very large numbers—much beyond every-day practice—to reduce 
in a marked manner this variation of the opsonic index from unity, when an 
individual even is tested against himself. Generally the diagrams indicate that 
an extreme variation in excess is more likely to occur than an extreme variation in 
defect, but that the most probable index, when the individual is tested against 
himself, will be one somewhat less than unity. 
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MISCELLANEA. 


I. On the Value of the Teachers’ Opinion of the General Intelligence 
of School Children. 


Compiled from the Tables and Reductions of 


H. GERTRUDE JONES, University College, London. 


A great deal of random criticism has recently been made of school-teachers’ estimates of 
general intelligence. We have been told that those estimates are absolutely without value, that 
any correlations obtained from them are idle and that the personal equation involved is so great 
that no use can legitimately be made of them. Good illustrations of this type of criticism 
will be found in the notices of Heron’s memoir on The Influence of Defective Physique and 
Unfavourable Home Environment on the Intelligence of School Children, which Mr G. U. Yule 
has considered it desirable to publish in two separate journals*, As it does not appear that 
those who criticise the teacher’s estimate of general intelligence in this way have made any 
experimental inquiry into the matter themselves, it may be of value to publish an investiga- 
tion of a preliminary character made some time ago as an attempt to ascertain whether such 
classifications are really idle. More ample researches will shortly be published, but it may be 
said in advance here that they serve to confirm the present result. 


The data are taken from schedules filled in by over 20 Aberdeen teacherst, who were not 
specially prepared for the task. 249 boys were taken from 4 different schools ; their ages ran from 6 
to 14 and they belonged to classes which were termed “ Infants,” I, II, III, 1V, V, and VI. They 
were judged not only by different teachers, but in different schools, and their mental capacity 
was appreciated in four groups: (i) Excellent=E., (ii) Good=G., (iii) Moderate = M., (iv) Dull=D. 
This is a pure appreciation of general intelligence by a variety of teachers in a variety of schools. 
Next we have taken the individual boy’s place in class as shewn by examination results and 
divided by the number of boys in the class. This may be adopted as a measure of the boy’s 
examinational intelligence. If the teacher’s estimate of general intelligence be of small value we 
should expect that when allowance is made for difference of class and age there would be little 
relation between examinational test and general intelligence. 


The following characters were taken out and tabled, ‘‘ Mental Capacity,” Age, Position in 
Class (measured as above stated) and Class or Standard, and the following six tables obtained. 


* Journal of the Statistical Society, Vol. txxi11. p. 547, and School Hygiene, Vol. 1. p. 473. 
+ For the Royal Commission on Physical Training (Scotland). We have gratefully to acknowledge 
the loan of these schedules by Professor Matthew Hay. 
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TABLE I. 


Mental Capacity. 




















Age | G@ | M | D | Totals 
6 4 uf|ul 5 31 
v 2 15 | 10 6 33 
8 5 19 6 4 34 
9 3 18 12 — 33 

10 2 20 11 1 34 

1l 2 9 12 5 28 

12 1 11 10 3 25 | 

13 — 9 13 | 2 24 | 

14 — | 4 3);— " 

Totals 19 116 | 88 26 249 
TABLE II. 
Mental Capacity. 

Standards E and G | M and D | Totals | 
Infants, I, II, III 77 61 138 
IV, V, VI me 58 53 lll | 

as... sn 135 | 114 249 | 
| 











TABLE III. 


Class or Standard. 














Age Infants| I IT [II IV V VI | Totals 

6 30 1 = Sy oe Bd 31 

? — 30 2 1 oo = ra 33 

8 _ 12 21 1 se > ee 34 

9 me 3 3 21 6 ae he 33 

10 _ i i 10 22 l 34 

11 — <a 1 — 4 22 1 28 

12 _ x =e 2 9 13 25 

| sas — — — — 7 17 24 
| ms _ ‘ont -|— - on 7 7 
| Totals 30 | 46 27 35 33 | 39 39 | 249 
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TABLE IV. 
Mental Capacity. 
Placein | = | a | M | D | Totals 
3 coe 2 a 18 
‘08 3 § a 11 
13 4 10 1 ts 15 
18 3 12 | — 17 
23 1 12 sh 17 
28 s 11 1 | — 12 
33 ot 10 an 13 
38 = 7 oa ll 
“43 aa 5 5 | 10 
48 i 11 ST 20 
53 = 8 a a 14 
58 1 5 “oa 15 
63 — 2 7 | — 9 
68 = 4 6 gi 10 
13 = 2 4 3 9 
"7 ai dine 8 l 9 
3 = des 11 es ll 
“88 ws ‘ate 5 5 10 
“93 = ~_ 2 6 8 
‘98 <a on 2 s 10 
Totals 19 | 116 88 26 249 
TABLE V. 
Age. 
lace i | | 
— 6 7 | 8 | 9 | 10 11 | 12 | 13 | 14 | Totals 
i u 
‘03 2 si—]| 5 3 2 3|/—|— 18 
08 3 1 “a 2 1 1 1 =e 11 
13 S43 se 1 ne 1 2 1 15 
18 4| 1 ee Sa 1 — i}; 3} — 17 
23 i 3 6 1] 1 2 a oS 17 
28 as | 2 3 1 4 = 1 | ae 12 
“33 —j—} 4 3 4 a SE a < 13 
38 “ae oe eo ae 3 1 il - 2 11 
£43 as a ane 1 1 2 2 ee 10 
48 3 a 1 1 4 4 2 | 2 20 
53 1 a 3 us 3 = = | 1 14 
‘58 1 — | 3 2 1 4 2 Si = 15 
63 l 2} 1 l 4 en ee ee 9 
68 — | l | 2 2 3 — ae 2\;— 10 
"3 |} 3 2 iat io 2 , ae 9 
7 2 Se 2 Pj—_|—]-at— 9 
“83 ie oa 2 2 8 1 Fs 3 11 
“88 —| 1 1 1 1 | 32 1 | 2) 1 10 
93 2] 1 1 a | an ft <2 2 . ha 8 
‘98 OSA Cae foot ae a oe oP 10 
| | | | 
Totals 31 | 33 | 34 28 25 2 | 7 249 
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: TABLE VI. 


Class or Standard. 








l l 
Place in | Infants| I 1 | ur | iv | Vv | VI {| Totals 
Class | 
| rad 
03 2 | 3 — | 5 5 |} 2 1 18 | 
‘08 3 | 1] — 1 2 2 2 ns 4 
13 “sa ie & 2 1 2 2 15 | 
18 ey .34 % 2 2 2 3 17 
23 i. he 4 3 — | 3 3 17 
28 a 3 2 1 5 | 1 12 
3: aes 1 4 5 1 1 l 13 
38 £5 we 2 1 2 | 1 3 11 
“43 ae 1 l 1 3 2 10 
48 .s 2 1 | 5 | 6 | 20 
53 4 ee 3 1 2 1 14 
58 eos 2 1 is 2 15 
63 1 1 2 2 ae le ee “ir 
68 a 2 1 2 3 l 1 10 
13 1 Th eee Bae 1 et 2 Boa 
x8 1 3 l 2 Li = 1 9 
83 1 “pe 1 a 2 4 11 
88 ies 1 2 l l 3 2 10 
93 “Ts ae 1 1 ae l 2 8 
98 4 “Ae ee 2 2 10 
| | | | | 
| | | | | 
Totals 30 | 46 | 27 | 35 | 33 | 39 | 39 | 249 
j | | | 








Let us look at the results of these Tables in succession : 
Table I. Mental Capacity and Age. The mean ages are as follows: 
E: 887 | 
V. a- 
G: 9 89 | All classes: 9°91. 
M: 10°33 
D: 9°38 


There is thus a decrease of intelligence with age until we come to the worst class and here 
there is a drop in age. There are very certainly cross-currents at work, e.g. clever boys go 
young to school, and leave early, and extreme cases of dullness go to “special” schools at 
9 to 11, and school lessens at any rate the outward appearance of stupidity. Worked out by the 
correlation ratio method the correlation is 7=‘1830. This clearly should be considered as 
negative, but while the regression is not linear and therefore the correlation ratio sensibly 
greater than the correlation coefficient, we will take r= —*1830 because we believe the “ dullness” 
of the younger children to some extent wears off with continued school life or the dull are drafted 
elsewhere*. 


* The reader should examine the last column of Table I. We think several of the younger children 
would ultimately pass into ‘‘ special schools” or be sharpened by school life. The same bi-modal 
character is obvious in column M, if to a lesser extent. 


4 
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Table II, Mental Capacity and Standard. It is interesting to look at the data for this as a 
contingency table. This is given below as Table II bis : 
TABLE II bis. 


Class or Standard. 














. | 

bey Infants | I II Ey kv iow | VI {| Totals 
‘3 

&.| E a; 1) aa is ea 8 19 
35 | G i + 15 | 18 | 21 | 10 22 116 
ae 10 l4 7 | 15 | 10 | 18 | 14 88 
8\ D 5 2 a 8 2 26 
s Totals] 30 6 | 27 | 35 | 33 | 39 | 39 | 249 








The greatest contributions to the contingency arise from (i) the defect of Dulls in 
Standard III and (ii) the excess of Dulls of Standard V, accompanied as it is by a defect of 
Good mental capacities. These may be personal to the teachers of these standards, or to their 
selection of children, which was made for the purposes of physical measurement. But the 
corrected contingency coefficient which should be much emphasised by personal equation is only 
*2123*, Now we shall see that the correlation between age and standard is very high, i.e. ‘9366; 
hence the correlation between mental capacity and age for constant standard = — ‘0461, or is 
practically negligible. In other words, allowing for the personal equation of the master with 
regard to intelligence, age counts for very little in his estimate of general intelligence. If we 
even out the irregularities of the above table by forming out of it a mere four-fold table as 
in Table II we obtain r= — ‘0558. We expect that this is really the more correct measure of the 
true relationship between standard and general intelligence. 

Table IIT. Age and Standard. The correlation here is clearly very high ; determined as a 
correlation ratio it is ‘9366. Practically children are moved up a class or standard a year. 

Table IV. Mental Capacity and Place in Class. The correlation as found by correlation 
ratio is *7246. The actual mean places of each category are : 

E: -122 | 

G: as + All categories: °456. 
M: ‘500 . 

D: ‘851 

This list itself shews that the teacher’s estimate of general intelligence represents very closely 
the examinational place which a student will take and for many inquiries may quite legitimately 
replace that estimate of ability. 

Table V. Age and Place in Class. Found by the ordinary product moment method the 
correlation coefficient is ‘0484, shewing that with increased age a slightly worse class position is 
taken, i.e. any advantage from age is counterbalanced by the clever elder boys leaving school 
earlier. The value, considering the numbers dealt with, is, however, hardly significant. 

Table VI. Place in Class and Standard. The mean places are as follows : 

Infants *457 
I "494 
II "445 | 
Itl 393 » All classes: °456. 
EV "388 
b ‘493 | 
VI “493 


* This must be taken negative if treated as a correlation coefficient. 


| 
| 
| 
| 
| 
| 
| 
{ 
| 
| 











| 
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This series is very irregular. The » for the series is *1553, but this is rather a measure of the 
irregularity than of anything else*. If we take the correlation coefficient for a four-fold table 
as follows 


> 


Infants 


I, Ul, Il i V, V1} Totals 








05 to *45 71 53 124 
‘50 to 1°00 67 58 125 
Totals 138 = (| ol 249 








we find r=0580, a number of scarcely any practical significance 


and shewing that really place 
in class is, as it should be, independent of standard. 


We have then to sum up the following results : 


Intelligence and Age : ‘ . 1 3= — ‘1830 
Intelligence and Position in Class . r= *7246 
Intelligence and Standard : .  yg=— 0558 (— 2123) 


Position in Class and Age 'og== — 0484 
Position in Class and Standard - Ta= ‘0580 (*1553) 
Age and Standard 134 9366 

In the case of two of these correlations for which some doubt might exist, we have placed in 
brackets the higher values found by different processes. The results for both values of the 
partial correlation coefficient have been worked out. 

Correlation of “General Intelligence” and Examinational Place, for constant age and 
constant standard : 


With lower values of ryy and 74: py ="70+°02. 


With higher values of 744 and ry: pig ='99+°01. 


We may safely draw from these results the conclusion that the teacher’s estimate of general 
intelligence is not a purely idle character, wholly valueless owing to the personal equation of 
the teacher. Whatever else it may be, it is highly correlated with the place which the student 
will take in an examination test and accordingly has at least one quite definite significance. 


The first of these results shews that with the lesser and more probable values of the correla- 
tions we obtain, allowing for standard and age, a high correlation (‘70) between the teacher’s 
estimate of general intelligence and actual examination measure of capacity. There are, 
however, cross-currents at work in elementary schools ; the one is the selection of the notably 
dull children about 10 who are drafted into “special schools.” This causes rather a defect of 
dull children in the II, III and IV standards ; the next is the selection of the more intelligent 
children in the highest standards (this is most obvious in V) to leave the elementary schools. 
The total effect is to make a somewhat high correlation ratio for standard and intelligence. If 
we took all the children in a class there ought to be sensibly no correlation between position 
in class and standard, but we actually find for classes II, III and IV excess of good places 
(averages *445, 393, °388 instead of -456) and for the high classes V and VI excess of bad places 
(averages *493, -493). This seems to be the result of the same cross-currents, the selection of the 

* It must be remembered that not every child, but only 30 to 40 in each standard were physically 
examined, and the number so dealt with is not by any means the number in the class. This selection of 
children accounts for the mean place in class being 456 and not -500. 
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dull to go and of the dull to stay at different ages. If we now give the higher values to the 
correlations produced by these tendencies in elementary schools, we find on correction not that 
they have tended to disguise the real value of the “general intelligence” character by giving it 
a spurious high intensity, but that they have tended to /ower its true importance. Those who 
complain of the different percentages of the divers grades of intelligence of the children in the 
different standards of the same elementary school often seem unaware that the same sample does 
not pass through all standards, at one stage there is a rejection of the very dull and at another 
the retention of the dull class. The main facts brought out are the substantial correlation 
between teacher’s estimate of general capacity and examination test, and the fact that extreme 
allowance for age and standard tends to emphasise this relationship rather than to shew that 
the teacher’s estimate is of little value. These observations were not made ad hoc, but they 
have been confirmed by other observations made on far larger numbers with additional safeguards 
for accuracy. The results of these will shortly be published, but it seems worth while indicating 
at a time when hasty criticisms are being made of the value of “general intelligence” estimates 
that the teacher’s appreciation of mental capacity does mean something, and has a very direct 
and practical value, when properly registered and handled. 


EB. P. 


II. Note on the Separate Inheritance of Quantity and Quality 
in Cows’ Milk. 


3y KARL PEARSON, F.RS. 


Under the above title Professor James Wilson, of the Royal College of Science, Dublin, read 
a paper in May, 1910, which is published in Vol. x11. pp. 470—479, of the Scientifie Proceedings 
of the Royal Dublin Society. 


He states that: “It is a very general opinion that the milk of high-yielding cows is usually 
poorer and that of low-yielding cows richer in quality” (p. 470), and again he writes : 

“If we group together all the low-yielding cows, and find their milk invariably high in quality, we 
may infer that low yield and high quality are of the nature of concomitant variations. If we group the 
high-yielding cows together, and find their milk invariably of low quality, we may infer that high yield 
and low quality run together. But if we take these groups and any other groups we can form, and find 
that the quality varies the same way in them all—that is that there are low qualities, high qualities, and 
medium qualities in every one of them—then we are justified in inferring that the quantity and quality 
of the milk are independent of each other. And this is what we do find” (p. 471). 

The italics are mine. It would be difficult to find a paragraph containing more fallacies in as 
many words. Two qualities may be associated together, even closely correlated, and yet one will 
not be invariably accompanied by the other. High, low and medium values of one character 
may occur with any given value of a second and yet there be high correlation between them. 
The whole problem turns on the extent of the variation of tie first character for a given value of 
the second, and not the invariable appearance of one with the other. The author may have 
some conception of correlation other than absolute association, but neither this paragraph, nor 
his general treatment of the subject, shews any signs of it. 


Professor Wilson’s material is involved in the table below. 


Now there can be no doubt of the value of such returns and of tabling them, although we 
regret the clubbing together of the tail frequencies, that common error of the non-statistical 
trained mind, which renders determination of the true correlation so difficult. 
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How now does Professor Wilson extract from these data the independence of quality and 
quantity of milk on which he insists, and which is chosen as the very title of his paper? He 
has prepared six elaborate diagrams, and he leaves the reader from mere inspection of these to 
come to the same conclusion as himself, i.e. that the variation is the same in all his groups and 
the quality independent of the quantity. Not being able myself to judge from such diagrams 
whether there is complete independence between quality and quantity, I thought it worth while 
to ascertain the mean percentages of fat and the variations in fat of his four groups. The results 
are as follows, the probable errors being added : 





Quantity of Milk Means Standard Deviations 


Under 500 gallons oie 3°710+ ‘010 372 + 007 
From 500 to 600 gallons ... 3°673 + ‘007 *319 + °005 
From 600 to 700 gallons ... 3°651 + ‘008 *321+ °005 
Over 700 gallons ... vax 3°640 + ‘008 ‘275+ 006 


The means and variabilities are in percentage of butter-fat. 

For the whole population of 2866 cows the mean is 3°672 and the variability is measured by 
3230. It will accordingly be clear that there is (i) a slight but perfectly definite decrease in 
quality with increase of quantity, the low-yielding cows have a higher, the high-yielding cows a 
lesser fat percentage than the average Ayrshire, although the increase and decrease only amount 
to about 1°/, of the average yield, and (ii) the variation in quality is not the same for cows of 
each quantity category. Cows which give high quantity have far less variation in their yield of 
butter-fat than cows of low yield. The number of cows is so large relatively that the probable 
errors are small and there is no doubt about the results being significant. Owing to the manner 
in which the tails of the above frequency have been clubbed together it is not possible to work 
out the proper correlation by the ordinary product moment method. The means, however, 
suggest that the correlation cannot be very far from linear, and the correlation ratio » may be 
fairly taken to measure r, This gives us r= —°075 with a probable error of +°013, a result, 
if small, yet definitely significant. 


It seemed worth while putting these numerical results on record, as a further protest against 


somewhat dogmatic inferences drawn from mere inspection of elaborate but unnecessary graphs. 
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